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Abstract  of  Dissertation  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  Doctor  of  Philosophy 

NONLINEAR,  UNBONDED,  POST-TENSIONED  TENDON  ELEMENT 
FOR  CONCRETE  STRUCTURES 

By 

Yuh-Yi  Lin 

Chairman;  Dr.  Marc  I.  Hoit 
Major  Department;  Civil  Engineering 

A general  tendon  element  is  developed  to  model  the 
prestressing  reinforcement  for  the  analysis  of  either 
pretensioned  or  post-tensioned,  prestressed  concrete 
structures.  Two  versions  of  the  element  are  provided.  One 
is  the  tendon  element  without  friction,  which  can  be  used 
for  bonded  prestressed  concrete  structures.  The  other  is 
the  tendon  element  with  friction,  which  is  mainly  developed 
to  simulate  the  external  prestressing  tendon  used  in 
post-tensioned,  prestressed  concrete  structures. 

The  formulation  of  the  tendon  element • without  friction 
is  based  on  the  continuity  along  the  tendon.  According  to 
this  criterion,  the  tendon  stiffness  matrix  better  models 
real  structures.  The  initial  force  matrix  is  formed  by 
using  force  eqilibrium  at  each  node.  The  prestress  loss  due 
to  the  deformation  of  the  structure  is  calculated  from  the 
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nodal  displacements.  After  the  loss  is  calculated,  the 
final  tendon  force  is  then  evaluated. 

The  immediate  losses,  due  to  the  friction  and 
anchorage  slip,  are  included  in  the  tendon  element  with 
friction.  After  accounting  for  the  losses,  the  initial 
force  matrix  is  then  formed.  Because  of  the  existence  of 
friction  at  the  harped  points,  the  tendon  element  will  be  a 
combination  of  truss  and  continuous  elements.  The 
structural  model  is  determined  from  the  sliding  condition 
at  each  node.  Then,  the  formulation  of  the  stiffness  matrix 
can  be  determined  from  the  structural  model  of  the  tendon 
element.  A force  distribution  method  is  proposed  to  solve 
for  the  tendon  force  including  the  effects  of  friction. 

This  method  uses  the  unbalanced  force  matrix  and  the  slip 
stiffness  matrix  in  which  the  friction  effects  are 
considered. 

Four  examples  are  presented  to  verify  the  validity  and 
applicability  of  the  elements.  The  effects  of  friction  on 
the  strength  of  the  structure  are  shown  by  a parametric 
analysis . 
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CHAPTER  1 
INTRODUCTION 

1 . 1 Background  of  External  Prestressina 

The  development  of  segmental  bridge  construction, 
mostly  pioneered  in  Europe,  has  been  spreading  into  the 
United  States  and  elsewhere  in  the  world.  The  major 
advantages  of  this  type  of  bridges  are  a reduced  total  cost 
and  duration  of  construction,  especially  in  medium  and  long 
span  bridges  [1,2,3]. 

Most  of  the  segmental  bridges  use  post-tensioned 
prestressing  cables  buried  inside  the  concrete  sections. 
They  may  either  be  bonded  or  unbonded.  In  such  systems,  the 
ducts  of  adjacent  segments  have  to  be  carefully  positioned 
in  order  to  match  precisely  at  all  joints.  Extra  effort  has 
to  be  made  to  place  the  ducts  properly.  In  addition,  the 
ducts  interfere  with  casting  and  consolidation  of  concrete 
[4,5].  For  economical  reasons  and  quickened  paced 
construction,  external  tendons  have  been  adopted.  In  recent 
years,  several  bridges  have  been  built  with  prestressing 
tendons  placed  outside  the  concrete.  Due  to  the  many 
advantages,  an  increased  usage  of  external  prestressing  may 
be  expected. 

The  concept  of  external  prestressing  was  introduced  by 
Dischinger  [4].  He  used  them  in  the  first  two  prestressed 
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concrete  bridges  ever  built  between  1928  and  1936.  In 
France  and  the  United  Kingdom,  external  prestressing  has 
been  used  for  the  past  several  years.  Jean  Muller  [4]  was 
the  first  engineer  in  the  United  States  to  apply  the 
concept  of  external  prestressing  to  a series  of  bridges, 
such  as  Long  Key,  Seven  Mile,  Channel  Five  and  Niles 
Channel  bridges  of  the  Florida  Keys. 

Not  only  the  erection  of  new  bridges  but  also  the 
repair  of  existing  bridges  have  contributed  to  the 
development  of  external  prestressing  techniques  [4].  When 
existing  bridges  need  additional  prestressing  and  no  empty 
, pre-buried  duct  is  available,  external  prestressing  is 
the  only  way  to  solve  the  problem.  For  new  bridges,  one  of 
the  advantages  of  external  prestressing  is  the  flexbility 
to  replace  the  cables. 

The  design  of  externally  prestressed  bridges  can  be 
divided  into  three  types  [4].  The  first  is  wholly  external 
prestressing  tendons  which  are  continuous  over  several 
spans.  The  Long  Key  Bridge  and  the  Seven  Mile  Bridge  in 
Florida  are  good  examples  of  this.  They  were  constructed  by 
using  the  span-by-span  erection  method.  In  this  form  of 
construction,  the  precast  segments  are  assembled  using 
post-tensioned  cables.  Prestressing  tendons  run  though  the 
inside  of  the  box  girder  and  are  anchored  in  the  pier 
segments.  Some  other  bridges,  such  as  Vallon  des  Fleurs,  La 
Banquiere  in  France  and  Bubiyan  in  Kuwait,  were  built  in  a 
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similar  way  [4].  The  profile  of  external  tendons  used  in 
Seven  Mile  bridge  is  shown  in  Fig.  1.1. 

The  second  type  is  wholly  external,  non-continuous 
prestressing.  An  example  is  the  Sermenaz  Viaduct,  built  by 
Bouygues  using  the  cantilever  method.  The  cable  pattern 
includes  the  classical  cantilever  and  closure  tendons  which 
are  placed  outside  the  concrete  and  anchored  in  vertical 
ribs  provided  in  each  segment.  According  to  Virlogeux  [4], 
this  type  needs  more  anchorages  and  is  less  attractive  than 
the  other  two  groups. 

The  third  one  is  mixed  prestressing.  The  method  was 
developed  for  the  La  Fleche  bridge.  Internal  tendons, 
confined  in  the  upper  slab,  were  designed  to  resist  the 
deadload  only.  After  deck  completion,  external  cables  were 
used  as  continuous  prestressing  tendons  and  designed  to 
balance  all  other  loads  [4]. 

Advantages  of  the  external  prestressing  system  have 
been  disscussed  by  Virlogeux  [4]  Sowlat  and  Rabbat  [5]. 

They  can  be  explained  as  follows:  First,  duration  of 
construction,  a major  economic  concern  for  all  construction 
companies,  is  reduced  by  eliminating  duct  laying  and 
adjustment  operations.  Second,  the  layout  of  the  tendon  can 
be  made  more  streamlined  and  the  efficiency  of  the 
prestressing  is  thus  improved.  This  is  so  because  wobble 
friction  loss  does  not  exist  and  the  total  friction  losses 
are  reduced.  Third,  placing  and  consolidation  of  concrete 
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Fig.  1.1  External  Prestressing  in  Seven  Mile  Bridge 


5 

are  simplified  and  web  thickness  and  the  dead  weight  of  the 
structure  are  reduced.  Fourth,  the  injection  procedure  is 
simplied  and  the  problem  of  concrete  intrusion  in  the  ducts 
during  casting  is  eliminated.  Therefore,  it  is  easier  to 
check  the  injection  process  and  to  access  and  repair 
inadequately  grouted  areas.  Finally,  replacement  of  broken 
and  corroded  tendons  are  possible  with  proper  design. 

1.2  Research  Needs 

Although  external  tendons  have  come  into  use 
gradually,  very  little  information  about  testing  and 
modeling  is  avaliable.  Sowlat  and  Rabbat  [5]  tested  the 
segmental  concrete  girders  with  external  tendons  and  show 
that  the  ultimate  strength  is  adequate  under  1983  AASHTO 
Specifications  [6]  or  ACI  318-83  Code  [7]. 

For  the  computer  analysis  and  design,  there  is  a need 
to  develop  the  model  of  the  external  tendon  according  to 
it's  specific  behavior.  According  to  the  AASHTO 
Specifications  [6]  or  ACI  318-83  Code  [7],  external  tendons 
belong  to  the  category  of  unbonded  tendons  and  the  strength 
of  the  member  is  calculated  based  on  these  provisions.  In 
reality,  however,  external  tendons  are  not  totally  un- 
bonded. They  act  semi-bonded  at  the  deviation 
blocks (harping  points)  [8].  Although  Soongswang  [9] 
developed  a model  for  external  tendons  assuming  complete 
slip  at  each  deviation  block,  further  research  about 
simulating  the  semi-bonded  behavior  of  external  tendons  is 
still  needed. 
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For  unbonded  tendons,  the  displacements  along  the 
tendon  are  not  continuous  because  of  lack  of  contact 
between  the  tendon  and  concrete.  Also  the  stiffness  of 
tendons  are  excluded  in  the  analysis.  But  for  external 
tendons,  the  tendon  does  come  into  contact  with  the 
concrete  at  the  harping  points,  and  then  the  deformation 
field  is  continuous.  Therefore,  the  stiffness  should  be 
included. 

1.3  Objectives  of  the  Study  and  Scope  of  the  Work 

The  objectives  of  developing  the  tendon  element  include 
the  following: 

First,  provide  a more  reasonable  and  precise  analysis 
for  external  prestressing.  In  general,  the  element  can  be 
applied  to  pre-tensioned  bonded  steels  as  well  as  post- 
tensioned  unbonded  steels. 

Second,  simulate  the  semi-bonded  behavior  at  the 
deviation  blocks  and  find  out  it*s  effects  on  the  whole 
structure.  The  friction,  induced  by  the  deformation  of  the 
structure,  is  used  to  model  the  semi-bonded  effects. 

Third,  provide  the  complete  procedures  for  the 
analysis.  This  contains  the  tensioning  process,  with 
friction  loss  and  anchorage  slip,  application  of  prestress 
load  and  external  loads. 

Fourth,  study  the  friction  effects  on  the  ultimate 
strength  of  the  structure.  When  the  heavy  load  is  applied 
on  the  structure,  the  friction  at  the  deviation  block 
becomes  critical  to  the  ultimate  strength  of  the  structure. 
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Fifth,  study  the  effects  of  the  tendon  profile  on  the 
behavior  of  the  structure.  The  tendon  profile  dominates  the 
deformation  and  strength  of  the  structure,  and  effects  the 
magnitude  of  friction  at  the  harping  point.  This  friction 
can  change  the  distribution  of  the  tendon  force  and  thus 
affects  the  behavior  of  the  structure. 

1.4  Arrangement  of  the  Following  Chapters 

Chapter  2 presents  the  literature  review  about  the 
computer  modeling  of  the  prestressing  steel.  The  tendon 
element  without  friction,  assuming  complete  slip  at  the 
harping  point,  is  discussed  in  Chapter  3.  Chapter  4 
explains  the  tendon  element  with  friction  and  presents  a 
method  to  calculate  the  tendon  force.  Tensioning  process, 
with  friction  loss  and  anchorage  slip,  are  included  in 
Chapter  5.  Examples  and  numerical  results  are  stated  in 
Chapter  6.  Conclusions  and  recommendations  are  included  in 
Chapter  7. 

Appendix  A is  the  input  guide  for  the  tendon  element. 
Appendices  B and  C provide  the  source  program  of  the  tendon 
element  without  and  with  friction,  respectively. 


CHAPTER  2 
LITERATURE  REVIEW 

2 . 1 Introduction 

The  rapid  development  of  external  tendons  has  given 
rise  to  the  need  of  precise  analysis  for  segmental  bridges. 
Computer  models  of  prestressing  steels  have  been  given  in 
previous  works  [9,10,11,12,13,14].  Truss  elements  or  frame 
elements  were  used  to  model  the  prestressing  steel  in  the 
analysis  of  prestressed  concrete  structures.  The  prestress 
losses  induced  by  friction  between  prestressing  steel  and 
duct,  elastic  shortening,  anchorage  slip,  shrinkage  and 
creep  of  concrete  and  relaxation  of  the  prestressing  steel 
were  considered  in  the  analyses.  The  different  behaviors 
between  bonded  and  unbonded  tendons  were  also  included. 

From  the  past  reseach,  there  is  no  specific  element 
available  to  model  external  tendons.  Although  Soongswang  [9] 
proposed  a 2-D  tendon  element  to  simulate  external 
presrtressing,  the  assumptions  of  the  element  can  not 
really  match  the  behavior  of  external  tendon. 

This  chapter  deals  with  the  past  research  of 
prestressing  steel,  including  application  of  prestress  and 
the  formation  of  stiffness.  The  following  sections  will 
discuss  models  of  tendons,  calculation  of  prestress  losses, 
stiffnesses  of  tendons  and  internal  tendon  force. 
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2 . 2 Models  of  Tendons 
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Kang  and  Scordelis  [10]  divided  a steel  tendon  into  a 
series  of  straight  segments  for  use  in  prestressed  concrete 
frame  analysis.  These  segments  were  modeled  by  truss 
elements  or  frame  elements.  The  discrete  representation  of 
the  prestressing  steel  or  reinforcement,  using  one- 
dimensional elements,  has  been  widely  used  [11].  The  layout 
of  tendon  in  the  element  is  shown  in  Fig.  2.1.  Greunen  and 
Scordelis  [12]  and  Millanes  and  Calzon  [13]  used  a similar 
method  to  model  prestressing  steel  in  the  analysis  of 
prestressed  concrete  slabs  and  composite  bridges, 
respectively.  Soongswang  [9]  presented  a tendon  element, 
assuming  complete  slip  at  the  harped  point,  to  model 
external  tendons  in  the  analysis  of  segmental  bridges.  The 
profile  and  the  degrees-of-freedom  of  the  element  is  shown 
in  Fig.  2.2. 

The  major  difference  between  truss  elements  and  tendon 
elements  is  in  the  formation  of  the  stiffness  for  the 
tendons.  The  basic  term  of  stiffness  of  a tendon  can  be 
expressed 

k=EpAs/L  (2.1) 

in  which 

Ep  = Young's  modulus  of  the  steel  tendon. 

Ag  = Area  of  the  tendon. 

L = Length  of  the  truss  element  or  total  length  of 
the  tendon. 
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Fig.  2.1  Prestressing  Steel  Segment 
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Profile  of  the  Tendon  Element 
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The  terms  of  the  stiffness  matrix  of  the  tendon  element  are 
coupled  to  each  other  due  to  the  continuity  along  the  whole 
length.  The  movement  of  each  degree-of-freedom  of  the 
tendon  element  will  contribute  the  elongation  of  the 
tendon.  But  the  movement  of  each  node  in  the  truss 
element's  model  only  affects  the  elements  connected  to  the 
node. 

For  the  tendon  with  a smooth  profile,  such  as 
parabola,  the  model  of  truss  elements  can  be  satisfied  if 
enough  elements  are  used.  This  is  because  the  angle  change 
of  two  neighbor  elements  is  very  small  and  can  be  ignored. 
But  for  a harped  tendon,  the  model  of  truss  elements  mat 
not  be  reasonable.  Because  prestressing  steel  is  basically 
a continuous  member,  the  tension  along  the  tendon  should  be 
the  same  excluding  the  friction  effects.  At  the  harped 
point  the  tension  of  each  side  will  be  the  same,  and  the 
tendon  is  continuous.  For  the  truss  element  model,  the 
harped  point  is  discontinuous.  Finally,  the  different 
models  will  yield  different  structural  behaviors, 
espacially  if  the  stiffness  of  the  tendon  is  taken  into 
account  in  the  analysis. 

2 . 3 Calculation  of  Prestress  Losses 

In  traditional  and  simplified  analysis,  only  the  key 
points  of  the  span  are  concerned  and  the  initial  force 
along  the  tendon  is  assumed  to  be  constant.  In  more  precise 
analysis,  the  varied  force  along  the  tendon  should  be 
considered. 
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For  pretensioning  bonded  tendons,  losses  of  prestress 
are  due  to  the  instaneous  and  long-term  deformations  of  the 
concrete.  In  this  case,  the  initial  force  along  the  tendon 
is  constant  and  equal  to  the  jacking  force.  The  losses  are 
obtained  internally  as  long  as  the  change  of  elongation  of 
the  tendon  has  been  evaluated  [13]. 

For  post-tensioned  structures,  losses  of  prestress  are 
due  to  friction,  anchorage  slip,  relaxation  and 
deformations  of  the  concrete  [15,16].  These  effects  are 
described  as  follows: 

2.3.1  Friction 

When  the  tendon  slides  through  the  duct,  friction 
resistance  is  developed.  As  a result,  the  tension  will 
decrease  as  the  distance  or  the  angle  change  from  the 
measured  point  to  the  jacking  end  increases  [15,16].  The 
following  equation  is  wellknown  to  calculate  the  friction 
loss 

Px  = Po  e"(kX+Ma)  (2.2) 

where 

Px  = Force  at  a distance  x from  the  jack. 

Pq  = Force  at  the  jack. 

k = Wobble  friction  coefficient. 

= Curvature  friction  coefficient. 

a = Angle  change  in  tendon  profile  over  x. 

The  prestressing  force  at  each  of  the  steel  segment  ends 
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can  be  calculated  by  using  the  above  equation  successively 
from  the  tensioning  end.  Kang  and  Scordelis  [10]  used  the 
average  of  the  values  at  the  two  ends  as  the  segment  force. 
2.3.2  Anchorage  Slip 

This  occurs  when  the  tendon  is  anchored  after  the 
force  is  transfered  to  the  concrete.  Typically,  the 
magnitude  of  slip  is  in  the  range  from  1/8  in.  to  1/4.  in. 
depending  on  the  prestressing  system  used  [15]. 

To  simplify  the  analysis,  the  slip  is  assiimed  to  be 
uniformly  distributed  over  the  length  of  the  tendon.  Then 
the  anchorage  loss  can  be  calculated  from  the  expression 
[15] 

^^anc  “ (aL/L) E (2.3) 

where 

AL  = amount  of  slip. 

L = tendon  length. 

Ep  = elastic  modulus  of  the  prestressing  steel. 

Actually,  the  loss  is  mostly  confined  to  a region 
close  to  the  jacking  anchorage.  For  more  elaborate 
analysis,  an  iterative  procedure,  accounting  for  friction, 
was  used  by  Greunen  and  Scordelis  [12].  The  length  affected 
by  the  slip  loss,  shown  in  Fig.  2.3,  can  be  expressed 
by  [17] 

^anc  “ 2AliAsEp/Po  ( 1 “ e”2^^anc)  (2.4) 

in  which 


q = (/ia+kx) /x 


(2.5) 
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Fig.  2.3  Prestress  Variation  Before  and  After  Anchorage 
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The  above  euation  is  based  upon  the  assvimption  that  the 
anchorage  loss  is  develped  over  length  La^c  average 

prestress  along  Lane  equal  to  the  average  of  Pq  and  P2. 
After  Lane  has  been  determined,  the  prestress  at  any  point 
along  the  length  Lane  he  evaluated  by  using  Eg.  (2.4). 

2.3.3  Relaxation  of  Steel 

When  the  prestressing  steel  has  been  held  for  some 
time,  there  is  a gradual  reduction  of  stress  in  the  steel 
resulting  from  relaxation.  The  amount  of  relaxation 
includes  the  intensity  of  steel  stress  as  well  as  time.  The 
ratio  of  reduced  stress  fp  to  initial  stress  fp^  can  be 
estimated  by  [18] 

fp/fpi  = 1 - log  t /10(fpi/fpy  - 0.55)  (2.6) 

where 

fpy  = 0.1%  offset  yield  stress, 
t = the  time  in  hours  after  stressing, 
and  fpi/ffpy  is  not  less  than  0.55. 

2.3.4  Deformation  of  the  Concrete 

This  deformation  is  caused  by  the  stressing  of 
tendons,  application  of  external  loads,  the  creep  and 
shrinkage  of  concrete  with  time.  The  losses  due  to  this 
deformation  is  automatically  covered  in  the  solution 
process  [13]. 

2 . 4  Stiffnesses  of  Tendons 
For  pretensioned  members,  the  stiffnesses  of  the 
tendons  are  used  throughout  the  whole  analysis.  This  is 
because  the  steel  and  the  concrete  are  already  bonded 
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together  when  the  prestress  is  released.  This  effect  has 
been  considered  in  past  work  [10,12,13,14]. 

For  post-tensioned  members,  the  prestress  is 
transfered  to  the  concrete  gradually  by  a tensioning 
process.  At  force  transfer,  the  stiffness  of  the 
prestressing  steel  is  not  included  when  the  prestress  and 
dead  load  are  applied.  After  force  transfer,  the  stiffness 
due  to  the  prestressing  steel  is  considered  if  the  tendon 
is  bonded.  For  unbonded  tendons,  there  is  no  composite 
action  between  the  steel  and  the  concrete  and  the  stiffness 
of  the  tendon  is  not  included  in  the  analysis.  In  this 
case,  the  unbonded  tendons  exert  a stay-like  action  that 
contributes  to  the  overall  stiffness  but  they  are  not 
included  in  the  stiffness  of  sections.  Millanes  and  Calzon 
[13]  used  the  following  eqation  to  obtain  the  stiffness 
matrix 

k = EpAp(Yp)T(Yp)/Lo  (2.7) 

Matrix  (Yp)  has  the  form: 

(^p)  “(^If^i/Ci, f^n'^n'^n)p  (2.8) 

Directional  coefficients  a^,  b^  and  c^  equal  zero  for  all 
nodes  not  crossed  by  cable  p and  have  the  values  for  the 
remaining  nodes. 

Soongswang  [9]  formed  the  stiffness  of  the  tendon 
element,  based  upon  the  continuous  behavior  along  the  whole 
length.  But  he  did  not  include  the  stiffness  through  the 
whole  analysis  due  to  the  unbonded  behavior. 
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The  above  algorithm  follows  the  classical  method  used 
in  the  prestressed  concrete  analysis  [15].  With  or  without 
stiffness  of  the  tendon  corresponds  to  the  net  or 
transformed  section  used  in  the  classical  analysis. 

2 . 5 Internal  Tendon  Force 

Kang  and  Scordelis  [10]  used  an  iterative  method  to 
find  the  difference  of  the  tendon  length  to  calculate  the 
tendon  force.  For  bonded  tendons,  the  strain  increment  is 
evaluated  by 

= (Lc“^)/^o  (2.9) 

in  which 

Lc/  %)  = "the  length  of  the  prestressing  steel  segment 
for  the  current  and  previous  iteration,  respectively. 

Lq  = the  original  length  based  on  the  original 
locations  of  the  two  end  points  in  the  frame  element. 

For  post-tensioned  unbonded  structures,  the  following  eqation 
is  used: 

A€  = (SLc  -SLp  )/2Lo  (2.10) 

where  the  summation  is  made  for  all  the  steel  segments  in 
the  tendon.  The  total  strain  in  current  stage  e is  obtained 
by  summing  up  Ae . And  the  corresponding  stress  can  be 
obtained  from  the  stress-strain  curve. 

Usually  the  deformation  field  of  the  bonded  tendons  is 
continuous  and  the  same  as  that  of  the  concrete  and  the 
strain  and  tendon  force  can  be  calculated  directly.  But  for 
the  unbonded  tendons,  the  deformation  field  is 
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discontinuous  everywhere  except  at  the  anchors  and  the 
tendon  force  can  not  be  calculated  directly. 

Van  Zyl  and  Scordelis  [19]  proposed  the  use  of  an 
average  extension  factor  which  makes  it  possible  to 
directly  calculate  the  change  in  force  resulting  from  a 
change  in  tendon  length  for  unbonded  tendons.  If  a unit 
force  is  applied  at  both  anchors,  the  tendon  force  in 
segment  i,  f^,  can  be  calculated  by  using  Eg.  (2.2)  when 
the  effect  of  friction  is  considered.  The  elongation  in 
segment  i due  to  the  unit  force  is  expressed  by 


ALi  - 1^  fi/ (Ep  Ap)  (2.11) 

The  total  elongation,  aL,  of  the  tendon  of  n segments  is 

aL  — E AL^ , 1—1 ,2,3,...,n  (2* 12 ) 

The  average  extension  factor  is  defined; 

■^ef  ~ (^p  Ap)/E(li  fi)  (2.13) 

Then  the  force  change  at  the  anchors  is  obtained: 

APq  — AL  A0f  (2.14) 


And  the  force  change  in  each  segment  can  be  calculated  by 
mutiplying  aPq  fi. 

Greunen  and  Scordelis  [12]  adopted  the  criterion  for 
unbonded  tendons  and  used  the  extension  ratio  in  the 
analysis  of  prestressed  concrete  slab.  The  ratio  is  defined 
as  the  extension  of  a frictionless  tendon  over  that  of  a 
tendon,  including  the  effect  of  friction  due  to  a unit 
force  applied  at  both  anchors,  and  is  given  by 
Re  = L/Z(li  fi)  > 1 ; i=l,2, ,n 


(2.15) 
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After  calculating  the  total  elongation,  aL,  and  using  the 
extension  ratio,  Rg,  the  change  in  anchor  force,  aPq/  is 
calcualted  as 

APq  = (Ap  Ega  Rq)AL/L  (2.16) 

where  Eg^  is  the  average  modulus  of  elasticity  of  the 
prestressing  steel. 


CHAPTER  3 

TENDON  ELEMENT  WITHOUT  FRICTION 
3 . 1 Introduction 

The  purpose  of  tendon  element  is  to  analyze  the 
behavior  of  prestressed  concrete  structures  which  contain 
prestressing  steel.  Because  of  the  different  properties  and 
behavior  between  concrete  and  steel  tendons,  it  is  more 
plausible  to  use  a tendon  element  combined  with  concrete 
elements  (shell  elements)  than  to  use  only  concrete  elements 
in  the  analysis  of  prestressed  concrete  structures. 
Basically,  the  element  can  be  used  in  pretensioning  or 
post-tensioning  tendons. 

The  tendon  element  introduced  here  has  4 or  6 
nodes (Fig.  3.1)  which  connect  to  a shell  or  membrane  element 
used  to  model  concrete  beams. 

To  be  capable  of  handling  the  3-D  prestressing  steel, 
the  tendon  element  has  twelve  degrees  of  freedom,  with  two 
at  node  II  and  III  and  three  at  node  JJJ  and  JJ  and  only 
one  at  node  I and  J,  in  the  local  coordinate  system.  In  the 
global  coordinate  system,  three  degrees  of  freedom  per  node 
are  required  because  of  the  3-D  capabilities. 

The  tendon  element  discussed  here  does  not  consider 
the  force  variation  at  the  harping  points,  that  is,  sliding 
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Fig.  3.1  Tendon  Element  in  the  Local  Axes 
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will  happen  at  all  intermediate  nodes.  Therefore,  the 
tendon  element  can  be  modeled  as  a continuous  tension 
element.  The  tension  at  any  point  along  the  tendon  is 
assumed  to  be  the  same.  For  external  prestressing,  more 
elaborate  modeling,  due  to  the  force  variation,  will  be 
discussed  in  the  next  chapter. 

3 . 2 Formulation 

3.2.1  Stiffness  Matrix  Without  Friction 

In  general,  the  sag  of  the  tendon  compared  to  the  span 
is  very  small,  therefore  the  change  of  the  angles  is 
negligible.  Small  strain  theory  is  applicable  in  this 
derivation.  The  derivation  of  the  stiffness  matrix  without 
friction  is  developed  in  a local  coordinate  system  q 
(Fig.  3.2). 

The  stiffness  is  derived  by  definition; 

Kij=the  force  at  DOF  i due  to  the  unit  displacement 
at  DOF  j.  When  qi=l  and  the  other  displacements  are  zero, 
then  the  deformation  of  the  tendon  is  a1=-1  and  the  force 
of  the  tendon  will  be  equal  to  -EA/L.  From  the  equilibrium 
of  each  node  (Fig.  3.2),  we  can  easily  find  the  forces  at 
each  degree  of  freedom  as 
At  node  I 

Kii=  EA/L 
At  node  II 

K21=  -EA/L(l-COSai) 

K3i=  -EA/L(cos/3i) 


ig.  3.2  Derivation  of  Stiffness  Matrix  (q=l) 


At  node  III 
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K4i=  -K21 
Ksi=  -K31 
At  node  JJJ 

Kei=  -EA/L(l-COSa2) 

K?i=  EA/L(cos/32) 

Ksi=  EA/L(coS72) 

At  node  JJ 

Kgi=  -Kgl 
Kl0/1=  “K71 
Kll/1=  -Ksi 
At  node  J 

^12  1 1=  -EA/L 

where 

E = modulus  of  the  tendon. 

A = area  of  the  tendon. 

L = total  length  of  the  tendon. 

“1/^1  = the  angles  between  segment  II-III  and  X 

and  Y axis,  respectively. 

“2/  Ph  72  = the  angles  between  segment  JJJ-JJ  and  X 

and  Y and  Z axis,  respectively. 
Similarly,  the  other  stiffnesses  terms  can  be  derived.  The 
stiffness  matrix  can  be  stated  as  follows; 
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3.2.2  Transformation  Between  Local  and  Global  Stiffnesses 
Let  [T]  represent  the  transformation  matrix  between 
local  and  global  displacements.  We  have 

{q}=[T]{r}  (3.2) 


where 


{q}  is  the  matrix  of  local  displacements, 
{r}  is  the  matrix  of  global  displacements. 
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In  the  same  way,  we  can  find 

{Q}=[T]{R}  (3.3) 

where 

(Q)  is  the  force  matrix  in  local  coordinate. 

{R}  is  the  force  matrix  in  global  coordinate. 
The  force-displacement  relationship  in  local  coordinate  can 
be  written 

{Q)=[K]L{q)  (3.4) 

From  Eg. (3.2) 

{Q)=[K]LCT]{r}  (3.5) 

and, from  Eg. (3.3), 

[T]{R}=[K]L[T]{r}  (3.6) 

Using  the  property  of  orthogonality, 

{R}=[T]T[K]L[T](r}  (3.7) 

Thus 

{R}=[K](r)  (3.8) 

where 

[K]=[T]T[K]l[T]  (3.9) 

This  is  the  equation  most  generally  used  for  finding  the 
member  stiffness  matrix  in  global  coordinates  [20]. 

3.2.3  Initial  Force 

For  pretensioned  members,  the  initial  force,  T,  is 
the  jacking  force  transfered  to  the  concrete.  For  post- 
tensioned  members,  the  initial  force,  T,  is  the  force 


28 


deducting  all  the  instaneous  losses  from  the  jacking  force. 
In  this  version,  the  force  variation  along  the  tendon  is 
not  considered  and  the  initial  force  should  be  estimated 
before  being  used  in  the  program. 

For  the  harped  tendon,  the  initial  force  is  the 
balanced  force  on  the  nodes  (Fig.  3.3).  The  direction  of  the 
initial  force  defined  in  ANSR  [21]  is  the  direction  that 
the  tendon  is  stressed.  For  example,  when  a straight  tendon 
is  being  stressed  in  the  simple  supported  beam,  the 
direction  of  the  initial  force  is  shown  in  Fig.  3.4. 

When  the  prestressing  force,  T,  is  applied  to  the 
tendon,  the  initial  forces  on  each  degree  of  freedom  can  be 
evaluated  by  the  equilibrium  equations  at  each 
node  (Fig.  3.3).  The  local  force  matrix  is  expressed  by 


1 X12  i 

[F]l  = T ! -Y2  i (3.1 
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i -X12 
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1 

T is  the  prestressing  force.  The  other  notations  are 
the  same  as  in  the  stiffness  matrix.  The  force  matrix  in 
global  coordinates,  from  Eg.  (3.3),  can  be  written 
[F]=[T]T[F]l 
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Initial  Force  ot  Tendon  Element 
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Direction  of  the  Initial  Force 
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3.2.4  Resisting  Internal  Force 

The  resisting  internal  force  is  defined  as  the  tendon 
force,  including  the  losses  due  to  the  displacements,  in 
the  current  step.  Therefore,  the  magnitudes  of  resisting 
internal  forces  are  the  initial  forces  minus  the  forces 
induced  by  the  nodal  displacements.  This  can  be  calculated 
by  sustituting  T-aT  into  T in  Eq.  (3.10),  where  aT  is  the 
force  due  to  the  nodal  displacements.  To  evaluate  AT,  we 
can  use  Eq.  (3.8)  to  get  the  force  matrix  R in  the  global 
system.  Because  the  tendon  is  continuous  along  the  whole 
length,  the  force  change  will  be  constant.  Only  the  force 
components  at  one  node  are  needed  for  calculation.  If  node 
I exists,  the  components  of  the  force  change,  AT,  are  R^, 

R2  and  R3.  Otherwise,  the  components  are  R4,  R5  and  Re.  The 
sign  convention  used  is  that  the  resisting  forces  are  the 
same  sign  as  the  initial  forces. 

3 . 3 Computer  Implementation 
ANSR  [21]  is  used  in  this  study.  Eleven  interface 
subroutines  are  called  from  the  base  program  to  perform 
element  computations.  Only  six  subroutines  are  utilized  in 
tendon  element  for  linear  analysis  without  frictional 
forces.  The  functions  of  each  subroutine  are  as  follows: 
3.3.1  Subroutine  INELll 

This  subroutine  is  called  from  the  base  program  once 
for  the  tendon  element.  The  purposes  of  the  subroutine 
include  the  following: 
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First,  input  and  initialize  the  element  group  data. 

The  labeled  common  block  /INFGR/  for  the  group  of  tendon 
elements  is  defined: 

COMMON  /INFGR  / NGR, NELS , MFST, IGRHED ( 10 ) , NINFC , LSTAT, 

* LSTF , LSTC , NDOF , DKO , DKT , ERROR (1000) 

in  which 

NGR  = Element  group  number. 

= 11  for  tendon  element. 

NELS  = Number  of  elements  within  the  group. 

MFST  = Element  number  of  the  first  element  in 

this  group. 

DKO  = Initial  stiffness  damping  factor. 

= 0 for  tendon  element. 

DKT  = Current  tangent  stiffness  damping  factor. 

= 0 for  tendon  element. 

IGRHED (10)=  Element  group  heading. 

The  above  are  read  from  the  element  input,  and  the 
others  stated  below  are  set  up  in  this  subroutine. 

LSTAT  = Length  of  words  of  state  information 
variables. 

= 108  for  tendon  element. 

LSTF  = Length  of  words  of  stiffness  information 

variables. 

=648  for  tendon  element. 

LSTC  = Length  of  words  of  stiffness  control 

variables. 
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= Length  of  words  of  variables  in  /INFEL/. 
= 756  for  tendon  element. 

NDOF 

= Number  of  element  degrees  of  freedom. 
= 18  for  tendon  element. 

The  element  input  data  is  read  next  and  stored  in  the 
common  block  /INFEL/.  The  labled  common  block  /INFEL/  for 
tendon  element  includes  the  following; 

COMMON  /INFEL/  IMEM, KST, LM ( 18 ) , NODE (6) , XX (6) , YY ( 6) , 

ZZ (6) ,PROP(l) , AREA, XL, B (7) ,SEL,VTOT, 
STOT,TT(3,3) ,T,SKP(324) 

in  which  IMEM,  NODE (6),  AREA  and  T are  read  from 
input,  and  the  others  are  calculated  in  the  program.  The 
definitions  of  these  variables  are  stated  as  follows: 


IMEM 

= Element  number. 

KST 

= Stiffness  update  code. 
= 1 (initialization) . 

= 0 (for  linear  case) . 

LM(18) 

= Location  matrix. 

NODE (6)  = Node  numbers  from  I to  J. 

XX  (6) 


YY(6) 

= Local  coordinates  of  the  six  tendon 

ZZ  (6) 

nodes . 

PROP 

= Young's  modulus  of  the  tendon. 

AREA 

= Cross  sectional  area  of  the  tendon. 

XL 

= Element  length. 

B(7) 

= Direction  parameters,  see  Eg.  (3.13). 
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SEL 


Current  tendon  force 


VTOT 


= Current  tendon  deformation 


STOT 


Current  tendon  stress 


TT(3,3)  = Transformation  matrix 


T 


Initial  force  of  the  tendon 


SKP(324)  = Current  stiffness  matrix  (=18*18) 


The  above,  except  SKP(324) , are  state  information 
variables.  SKP(324)  is  the  stiffness  information  variable. 
The  lengths  of  these  variables  are  calculated  and  saved  in 
COMMON  /INFGR/. 

The  third  step  is  to  initialize  the  global  degree  of 
freedom  numbers.  This  can  be  achieved  by  calling  subroutine 
NCODLM  eighteen  times.  Then  the  values  are  stored  in  LM 
array. 

The  fourth  step  is  to  transform  the  global  coordinates 
of  the  nodes  of  the  tendon  into  the  local  coordinates  by 
calling  subroutine  TRANl.  In  TRANl  the  transformation 
matrix  TT  is  calculated  by  calling  subroutine  VECTOR  and 
CROSS,  and  the  local  coordinates  can  be  easily  found  by 
using  the  following  formula: 

XX(N)=TT(1,1) *X(N)+TT(1,2) *Y(N)+TT(1,3)*Z(N) 
YY(N)=TT(2, 1) *X(N)+TT(2,2) *Y(N)+TT(2,3) *Z (N)  (3.12) 

ZZ (N)=TT(3,1) *X(N)+TT(3,2) *Y(N)+TT(3,3) *Z (N) 

The  subroutine  saves  the  direction  parameters  in  an 
array  B.  This  is  to  avoid  unnecessary  recalculation. 
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Due  to  the  geometry  of  the  tendon  element,  it  is 
easier  to  use  direction  parameters  in  order  to  find  the 
components  of  the  tendon  force  in  the  global  axes.  The 
direction  parameters  are  defined: 

B(l)  = cos/9i 
B(2)  = cos/92 
B(3)  = cosqi 

B(4)  = COS02  (3.13) 

B(5)  = COS72 

B(6)  = 1.  - COSai 

B(7)  = 1.  - COS02 


Finally,  the  stiffness  matrix  profile  is  calculated 
by  calling  the  subroutine  BAND  and  the  element  data  is 
transferred  to  tape  by  calling  the  subroutine  COMPAC. 

3.3.2  Subroutine  STIFll 

This  subroutine  performs  the  calculation  of  the 
element  stiffness.  Because  of  excluding  friction  force  in 
this  stage,  the  subroutine  is  only  called  when  the  total 
element  stiffness  matrix  is  to  be  formed.  The  purposes  of 
this  subroutine  are  stated  as  follows: 

First,  compute  the  local  stiffness  matrix  S (including 
the  six-node  tendon,  five-node  tendon  and  four-node 
tendon) . 

The  stiffness  matrix,  derived  in  section  3.2.1,  is 
calculated  in  this  subroutine.  For  the  general  purposes, 
the  stiffnesses  of  the  four-node  tendon  and  the  five-node 
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tendon  are  included  here.  These  tendons  are  distinguished 
by  the  coordinates  of  the  nodes,  and  their  stiffnesses  are 
modified  from  the  six-node  tendon. 

Next,  transform  the  local  stiffness  matrix  S,  twelve 
degrees  of  freedom,  into  the  global  stiffness  matrix  SK, 
eighteen  degrees  of  freedom.  The  transformation,  based  upon 
Eq.  (3.6),  is  calculated  by  calling  subroutine  TRAN2.  The 
size  of  the  transformation  matrix  is  12  x 18. 

Finally,  store  the  global  stiffness  matrix  with  the 
standard  form  (SKP)  in  the  element  common  block  INFEL. 

3.3.3  Subroutine  INITll 

This  subroutine  is  called  at  the  beginning  of  each 
load  step.  In  original  ANSR  the  prestressing  force  was  not 
included  in  the  options  of  initial  forces,  but  changes  have 
been  made  to  adapt  for  this  purpose.  Then  the  subroutine 
initll  will  be  called  by  the  base  program.  The  purposes  of 
the  subroutine  are  as  follows: 

First,  compute  the  initial  forces,  FI,  using  the 
prestressing  force  and  the  direction  cosines.  The  local 
initial  forces  of  the  tendon,  derived  in  section  3.2.2,  are 
calculated  by  calling  subroutine  FORVER. 

Finally,  transform  the  local  initial  forces  into  the 
global  forces  RF.  According  to  Eq.  (3.11),  the 
transformation,  from  twelve  degrees  of  freedom  to  eighteen 
degrees  of  freedom,  is  done  by  calling  subroutine  TRAN3. 


3.3.4  Subroutine  RINTll 
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The  subroutine  performs  the  element  force  recovery.  It 
is  called  in  the  beginning  of  the  analysis  and  after  each 
state  determination.  It*s  purposes  include  the  following: 

First,  it  evaluates  the  incremental  tendon  force,  DSL, 
due  to  the  displacements  of  the  nodes.  Because  of 
neglecting  frictional  forces,  the  change  of  tendon  force 
will  be  the  same  along  the  tendon.  Therefore,  only  the 
forces  at  the  end  of  the  tendon  are  concerned.  And  the 
force  components,  SA,  SB  and  SC,  of  the  end  node  are 
calculated  by  using  the  global  stiffness  matrix,  SKP,  and 
global  displacement  matrix,  q. 

Second,  it  calculates  the  step  loading  factor,  fac. 

Next,  the  routine  updates  the  total  tendon  force  SEL 
and  total  stress  STOT.  The  tendon  force  SEL  can  be  easily 
found  by  substracting  the  tendon  losses  DSL  from  the 
previous  tendon  force  which  is  the  initial  force  in  the 
linear  case. 

The  routine  then  computes  the  resisting  internal 
forces  FT  with  respect  to  the  local  coordinates  using  the 
tendon  force  SEL.  This  is  done  by  calling  subroutine 
FORVER. 

Finally,  the  local  resisting  internal  forces,  FT,  are 
then  transformed  to  global  resisting  internal  forces,  FE, 
by  calling  the  subroutine  TRAN 3 . There  is  no  damping  force 
in  this  element,  so  the  force  matrix  FD  is  set  to  zero. 


3.3.5  Subroutine  STATU 
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Subroutine  STATU  is  the  state  determination 
subroutine.  Since  the  element  is  linear,  there  is  no  state 
determination.  Therefore,  it  sets  the  stiffness  update  code 
KST=0  for  linear  analysis  without  frictional  forces.  This 
is  because  the  stiffness  matrix  is  not  changed  during  the 
analysis. 

3.3.6  Subroutine  OUTS 11 

The  subroutine  performs  the  output  of  the  tendon 
element.  It  outputs  element  current  state.  If  the  print 
indicator  KPR=0,  no  output  is  processed,  otherwise,  the 
output  will  include  tendon  force,  stress  and  strain  in  the 
current  state. 


CHAPTER  4 

TENDON  ELEMENT  WITH  FRICTION 
4 . 1 Introduction 

The  prestressing  tendons  used  in  this  study,  running 
inside  the  box  girder  (but  not  through  the  section) , are 
completely  unbonded  from  the  concrete  (Fig.  4.1).  The 
locations  of  the  tendons  are  controlled  by  the  deviation 
blocks  cast  monolithically  with  segments  [8,9].  Because  of 
the  high  tension  of  the  tendons,  a high  vertical  force  will 
be  produced  in  the  deviation  block  and  consequently 
develop  a large  friction  resisting  force  against  the 
movement  between  the  tendon  and  concrete.  This  friction  and 
the  bonding  stress  at  the  block  yield  the  semi-bonded 
behavior.  But  the  bonding  stress  is  small  compared  to  the 
friction  and  can  be  neglected  or  can  be  incorporated  in  the 
friction.  When  the  bridge  is  deformed,  the  displacements  of 
the  deviation  block  will  be  accompanied  by  the  deformation 
of  the  tendon.  This  deformation  will  induce  changes  in  the 
tendon  force  and  friction  force  and  finally  might  change 
the  structural  behavior  of  the  tendon. 

Usually,  friction  losses  in  this  type  of  unbonded 
tendons  can  be  divided  into  two  categories.  One  is  the 
initial  friction  between  tendon  and  duct  when  the  tendons 


39 


40 


cn 

fa 


Profile  of  the  External  Tendon 


41 


are  being  tensioned.  The  other  is  the  friction  between 
deviation  block  and  tendon  caused  by  the  deformation  of  the 
stucture.  The  former  is  always  taken  into  account  in 
analyzing  post-tensioned  prestressed  concrete  structures. 
However,  the  latter  is  often  ignored  by  most  engineers.  For 
bonded  tendons  this  is  reasonable,  because  sliding  between 
tendon  and  concrete  will  not  happen.  But  in  external 
tendons,  sliding  between  tendon  and  concrete  block  might 
happen.  In  this  case,  friction  plays  an  important  role  and 
cannot  be  neglected.  This  chapter  presents  a rational 
method  to  determine  the  tendon  force  and  find  the  effects 
of  the  friction  on  the  behavior  of  the  bridge. 

4.2  Structural  Model 

There  are  twelve  degrees  of  freedom  for  one  tendon 
element  in  the  local  coordinate  system  (Fig.  3.1).  Different 
from  the  frictionless  tendon,  the  force  along  the  tendon 
will  not  be  constant  but  varied  with  the  tendon  segment.  As 
in  the  finite  element  approach,  the  tendon  is  broken  into 
five  truss  elements  at  the  beginning  of  the  solution 
process.  After  the  displacements  have  been  solved,  the 
force  of  each  segment  can  be  found.  When  the  difference  of 
the  ajacent  segment  forces  is  larger  than  the  friction  at 
the  node,  sliding  will  happen  and  the  element  will  be 
treated  as  continuous.  On  the  other  hand,  if  friction  is 
not  overcome,  the  element  will  be  discontinuous  at  the 
nodes  and  should  be  regarded  as  a truss  element  (Fig.  4.2). 


Assume  T1  > T2 
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Therefore,  the  sliding  conditions  of  the  nodes  strongly 
affect  the  structural  behavior  of  the  tendon. 

After  each  node  has  been  checked,  the  tendon  will  be  a 
conbination  of  tendon  elements  and  truss  elements.  There 
are  many  combinations  of  the  tendon  as  shown  in  Fig.  4.3- 
4.5.  When  the  areas  of  the  tendon  become  large,  the 
differences  of  structural  behavior  between  tendon  element 
and  truss  element  will  be  obvious.  In  pratice,  the  external 
tendons  are  always  used  in  heavily  loaded  bridges  and  the 
areas  of  tendons  are  not  negligible.  The  stiffnesses  of  the 
tendons  also  form  an  important  part  of  the  structure.  For 
these  reasons  the  tendon  element  is  used  in  analyzing  the 
prestressed  concrete  bridge. 

4 . 3 Formulation 

4.3.1  Stiffness  Matrix 

The  stiffness  matrix  of  the  tendon  element  is 
dependent  upon  the  structural  model.  As  stated  earlier,  the 
element  matrix  is  composed  of  a combination  of  tendon  and 
truss  stiffnesses.  There  are  sixteen  combinations  of  the 
structural  model  as  shown  in  Fig.  4. 2-4. 4.  To  show  the 
derivation  of  the  stiffness  matrix,  consider  case  7 
composed  of  two  truss  elements  and  one  tendon  element. 

The  stiffness  is  derived  by  definition: 

Kij=the  force  at  DOF  i due  to  the  unit  displacement 
at  DOF  j.  When  q2=l  and  the  other  displacements  are  zero, 
then  truss  1 will  be  in  tension  and  the  tendon  will  be  in 
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Fig.  4.3  Types  of  Structural  Models  ( Part 
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Fig.  4.4  Types  of  Structural  Models  ( Part 
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Fig.  4.5  Types  of  Structural  Models  ( Part 
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compression.  The  deformations  of  truss  1 and  the  tendon  and 
truss  2 are  1 and  -cosai  and  0,  respectively.  The  force  of 
truss  1 equals  to  EA/L^.  The  force  of  the  tendon  will  be 
uniform  along  the  length  and  equals  to  -EA/L,  where 
L=L2+L3+L4.  From  the  equilibrium  of  each  node  (Fig.  4.6),  we 
can  easily  find  the  forces  at  each  degree  of  freedom  as 
At  node  I 
Ki2=  -EA/Li 
At  node  II 

K22“  EA ( l/L^+COSa^^/Lj 
K32=  -EA*COSai*COSj8i/L 
At  node  III 

K42=  EA*COSai*(l  - COSai)/L 
Ks2=  -K32 
At  node  JJJ 

Kg2=  -EA*COSai*(l  ~COSq2)/L 
Ky2=  EA*COSqi*COS/92/L 
Ks2=  EA*COSai*COS72/L 
At  node  JJ 

Kg2=  -EA*COSai*COSa2/L 
^10/ 2=  “K72 
Ki1/2=  -Ks2 

At  node  J 

Ki2/2=  0 

in  which 

“1/  /9l/  71  = the  angles  between  segment  II-III  and  X 

and  Y axis,  respectively. 


Truss  2 
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Fig.  4.6  Derivation  of  Stiffness  Matrix  (q2  =1 ) 
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= the  angles  between  segment  JJJ-JJ  and  X, 
Y and  Z axis,  respectively. 

Similarly,  the  other  stiffnesses  terms  can  be  derived.  The 
stiffness  matrix  can  be  stated  as  follows: 


’1  I2  <I3 

1a 

I5  16 

I7 

18 

I9 

I10  111 

I12 

L/Lil  -L/Li  1 0 1 

1 1 1 

0 

0 1 0 
1 

0 

1 0 

1 0 

0 1 0 
1 

1 0 

1 1 

X1*X11 

1 

Xl*Yl|-Xl*X12 

1 

X1*Y2 

|X1»Z2 

1 -X1»X2 

1 

-X1*Y2  |-X1*Z2 

1 

1 0 

1 1 
1 

-Y1*X11 

1 

-Y1  |n*X12 

1 

-Y1*Y2 

|-Y1*Z2 

1 Y1*X2 

1 

Y1*Y2  |Y1*Z2 
1 

0 

I 

1 

Xll^ 

I 

X11*Y1|-X11*X12 

1 

X11*Y2 

1 1 
1 X11*Z2|-X11*X2 
1 

1 

-X11*Y2|-X11*Z2 

0 

1 

Yl^  |-Y1*XX2 

1 

Y1*Y2 

1 Y1*Z2 

1 -Y1*X2 

-Y1*Y2  |-Y1*Z2 

0 

1 

1X122 

-X12*Y2 

|-X12*Z2 

1 X12*X2 

1 

X12*Y2|X12*Z2 

0 

CM 

CM 

1 Y2*Z2 

1 -Y2*X2 

1 

-Y22  I -Y2*Z2 

0 

Syrametric 

1 ZZ2 

1 -Z2*X2 

1 

-Z2*Y2  I -Z22  I 

1 t 

0 

IX22+L/L2 

1 1 

X2*Y2  |X2*Z2  1 

1 1 

-L/L2 

1 1 
Y2^  1 Y2*Z2  1 

0 

Z2^  I 0 


I L/L2 

in  which,  the  notations  are  defined  in  Chapter  3. 

After  the  structure  model  has  been  determined,  the 
number  of  sliding  nodes  can  be  evaluated  and  then  the 
tendon  stiffness  can  be  formed  and  put  in  the  corresponding 
block  of  the  12  x 12  element  stiffness  matrix.  The  truss 
stiffness  then  can  be  formed  and  added  to  the  corresponding 
degrees  of  freedom  in  the  element  stiffness  matrix. 

Whenever  the  structure  model  is  changed,  the  element 
stiffness  should  be  reformed. 


4.3.2  Initial  Force  Matrix 
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The  post-prestressing  force  is  applied  to  the  bridge 
in  the  form  of  equivalent  loading.  The  derivation  of  the 
matrix  was  stated  in  Chapter  3 and  omitted  here.  Fig.  3.3 
shows  the  direction  and  magnitude  of  the  initial  force.  The 
variation  of  the  initial  force  along  the  tendon  is  not 
shown  here  and  will  be  discussed  in  Chapter  5. 

4.3.3  Resisting  Internal  Force  Matrix 
4. 3. 3.1  Determination  of  truss  forces.  After  the 
displacements  are  solved,  the  deformation  of  each  truss 
element  can  be  calculated  by  using  the  formula 

Ai= [ (dxi+i-dxi) 2+ (dyi+i-dyi) 2+ (dzi+i-dzi) 2 ] 0 . 5 (4.2) 

where 

dxi=  X displacement  at  node  i. 
dyi=  y displacement  at  node  i. 
dzi=  z displacement  at  node  i. 

Then,  the  force  of  each  truss  element  can  be  found: 

Si=TIi+EAAi/Li  (4.3) 

where 

Si=  the  force  of  truss  element  i. 

TIi=  the  initial  force  of  the  segment  i. 

Ai=  displacement  of  truss  element  i. 

Li=  length  of  truss  element  i. 

4. 3. 3. 2 Determination  of  the  nodes*  sliding  conditions.  If 
the  force  of  each  segment  is  known,  then  the  friction  at 
each  node  can  be  calculated.  The  friction  model  used 
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assumes  that  the  friction  force  is  equal  to  the  normal  force 
times  the  coefficient  of  friction.  At  node  i,  if  the 
driving  force  DF  is  larger  than  the  static  friction  force 
SFF,  then  the  node  is  sliding.  Otherwise,  the  node  will  not 
slide.  Because  the  tendon  profile  in  the  deviation  block  is 
not  a strictly  sharp  angle  but  is  a smooth  curve.  However, 
the  transition  lenth  is  very  small,  a circular  curve  can  be 
assumed.  Therefore,  the  direction  of  the  normal  force  N is 
assumed  to  bisect  the  angle  S and  the  direction  of  the 
driving  force  DF  is  vertical  to  that  of  N as  shown  in 
Fig.  4.7.  The  magnitudes  of  the  forces  DF  and  SFF  can  be 
expressed  as  follows; 

DF=(Si+i-Si)cos(90-(?i/2)  (Si+i  > Si) 

N=(Si+i+Si)sin(90-<?i/2)  (4.4) 

SFF=/igN 

where 

^i+l/  Si  = the  forces  of  the  segments  connected  the 
node  i. 

/ig=  static  friction  coefficient  betwwen  tendon  and 
concrete. 

4 *3. 3. 3 Determination  of  tendon  forces,  if  both  end  nodes 
of  the  truss  element  do  not  slide  according  to  the  sliding 
criterion,  the  element  will  be  a real  truss  element  and  the 
force  Si  will  be  the  member  force  and  does  not  need  to  be 
revised.  Otherwise,  the  tendon  force  will  be  determined  by 
using  force  distribution  based  upon  the  stiffness  of  each 


DRIVING  FORCE  DF  :=(Sl-S2)«a)S(g(}-0/2) 
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Determination  of  Node's  Sliding — One  Node  Case 
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truss.  For  expressing  the  concept,  consider  the  case  with 
one  node's  sliding  which  is  shown  in  Fig.  4.7. 

4. 3. 3. 4 One  node's  sliding.  Let  Si  and  Si+i  designate  the 
force  and  of  truss  i and  i+1,  respectively.  Then 

Si=  TIi  + (AEAi)/Li  (4.5) 

Si+i=  TIi  (A®^i+l) /^i+1 

where  the  notations  were  defined  before.  Assume  Si  is  larger 
than  Si+i  and  the  driving  force  DF  is  larger  than  friction 
force  SFF. 

DFi=(Si  - Si+i)  cos(90  - ^i/2) 

Ni=(Si  + Si+i)  sin(90  -&i/2) 

SFFi=Ni/ig  (4.6) 

DFFi=Ni/ijj 

UBFi=DFi  - DFFi 
where 

^i=  the  angle  formed  by  member  i and  i+1. 

Md=  dynamic  friction  coefficient  between  tendon  and 
concrete. 

DFFi  = the  dynamic  friction  force  at  the  node  i. 

UBFi  = the  unbalanced  force  at  the  node  i. 

Then,  distribute  the  unbanlanced  force  UBFi  to  Si  and 
Si+i  by  using  the  distribution  factors  KRi  and  KRi+i.  The 
stiffness  of  each  member  is 


Ki=AE/Li 

Ki+i=AE/Li+i 


(4.7) 
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Since  the  friction  will  affect  the  force  distribution,  so 
the  contribution  of  the  friction  is  considered 

Kf=^d (Ki+i“Ki) sin ( 90-5  i/2 ) (4.8) 

Sum  up  the  stiffnesses 

ZK=Ki+Ki+l+Kf  (4.9) 

Then,  the  distribution  factors  can  be  defined: 

KRi=Kj^/ZK  (4.10) 

KRi+i=Ki+l/SK 

Because  S^  is  larger  than  Sj^+i,  the  friction  force  will 
reduce  S^  and  increase  Sf+i.  The  final  force  of  each 
segment  can  be  calculated  by 

Si=Si  - (UBFi*KRi)/COS(90-5i/2)  (4.11) 

Si+l=Si+i  + (UBFi*KRi+i)/cos(90-5i/2) 

For  the  tendon  without  friction,  the  tendon  force  can 
be  easily  found  by 

Si=Si+i=  TIi  + AEA/L=  T +AEA/L  (4.12) 

where  a— Aj_+Aj_+Q_,  L=L]_+L2+L3+L4+L5 

By  using  Eq.  (4.5-11)  and  pluging  Ms=Md=0f  we  can  get  the 
same  result  as  the  continuous  tendon  without  friction.  This 
shows  that  the  force  distribution  can  be  applied  for  the 
tendon  without  friction  as  well. 

4 *3. 3. 5 Mutiple  nodes'  sliding.  If  more  than  one  node  is 
sliding,  there  will  be  more  than  one  unbalanced  force. 
Because  of  the  existence  of  friction,  there  is  difficulty 
in  distributing  the  unbalanced  forces  with  regard  to 
direction.  The  unbalanced  forces  may  be  distributed  to  the 
nodes  in  opposite  directions,  thus  it  is  hard  to  decide  the 


final  direction  of  the  friction  force.  To  solve  the 
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problem,  a slip  force-displacement  relationship  is  formed. 
4. 3. 3. 6 Slip  force-displacement  relationship.  The  tendon  is 
assumed  as  the  truss  structure  at  the  beginning  of  the 
calculation  of  the  tendon  forces.  If  the  nodes  are  sliding, 
then  the  nodes  will  be  released.  Suppose  the  tendon  is 
modeled  by  the  structure,  with  a roller  support  at  each 
node,  shown  in  Fig.  4.8.  Since  the  element  is  continuous, 
the  degree  of  freedom  at  each  node  is  along  the  curve, 
connected  with  two  adjacent  segments.  In  this  model,  there 
is  no  external  force  but  there  is  a fixed  end  force  in  each 
segment.  The  force-displacement  relationship  [20]  can  be 
stated  as 

{R}  = [Kn]  {rl>  + {SF}  (4.13) 

where 

{R}  = the  external  force  matrix. 

= {0} 

[Kn]  = the  stiffness  matrix  of  the  truss  structure. 

{rl}  = the  slip  displacement  matrix. 

{SF}  = the  fixed  end  force  matrix. 

If  friction  is  considered,  then  the  above  is  modified 
to 

(0)  = [Kn]  {rl}  + {SF}  + {F}  (4.14) 

where 

{F}  = the  friction  force  matrix  between  tendon  and 
concrete.  The  friction  force  F includes  the  initial 
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Fig.  4.8  Slip  Degrees  of  Freedom  in  Tendon  Element 


57 


friction  FI,  due  to  the  fixed  end  force,  and  the 
incremental  friction  FD,  due  to  the  slip  displacement  (see 
Fig.  4.9).  Break  F into  FI  and  FD  and  then  rearrange 
Eg.  (4.14),  we  can  get 

(UBF)  = [Kn]  {rl}  + (FD)  (4.15) 

where 

(UBF)  = the  unbalanecd  force  matrix. 

= -{SF}-{FI}  (4.16) 

Introduce  the  slip  stiffness  matrix  [Kl]  and  Eg.  (4.14) 
can  be  written 

(UBF)  = [Kl]  (rl)  (4.17) 

If  SF^  is  smaller  than  the  maxmium  initial  friction 
Fli  at  node  i,  then  UBF^  will  be  egual  to  zero  and  the  node 
won't  slide.  Therefore,  the  degrees  of  freedom  of  the  above 
matrices  are  dependent  upon  the  number  of  sliding  nodes. 

At  most  there  are  four  degrees  of  freedom  in  one  tendon. 

The  unbalanced  force  on  degree  of  freedom  i can  be 
calculated  by  using  Eg.  (4.4-6). 

4. 3. 3. 7 Derivation  of  slip  stiffness  matrix.  The  slip 
stiffness  matrix  can  be  divided  into  two  parts;  One  is 
based  on  the  truss  stiffness  [Kn]  and  the  other  is  the 
friction  terms  [Kf]  induced  by  the  distributive  forces. 
Define 

Klij  = the  net  force  at  DOF  i due  to  the  unit 
displacement  at  DOF  j . 

Knij  = the  driving  force  at  DOF  i due  to  the  unit 
displacement  at  DOF  j . 


JOMII 
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Fig.  4.9  Free  Body  Diagram  of  Joint  II  Before  Force  Distribution 
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Kfij  = the  friction  force  at  DOF  i due  to  the  unit 
displacement  at  DOF  j . 

Nij  = the  normal  force  at  DOF  i due  to  the  unit 
displacement  at  DOF  j . 

Ki  = the  stiffness  of  each  truss  segment. 

= EA/Li 

When  rli=l  and  the  other  displacements  are  zero,  then 
the  force  of  each  segment  will  be 
Si=Kl 

S2=-K2  (4.18) 

53=84=35=0 

From  the  equlibrium  of  each  node  in  tangential  direction, 
we  can  find 

Knii=  Kl  + K2 
Kn2i=-K2 

Kn31=Kn4i=  0 (4.19) 

The  nomal  forces  induced  by  the  unit  displacement  are  found 
from  the  equlibrium  in  normal  direction 
Nii=  (K1-K2)  sin  a 

N2i=  K2  sin  a (4.20) 

N3i=  N4i=  0 

where  a=  90  - Bi/2 

Next,  the  friction  force  at  each  node  is  determined  by 

Kfii=  ± Nil 

= ± Md  (K1-K2)  sin  a 
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Kf2i=  T ^ld.  N21  (4.21) 

= T K2  sin  a 

Kf3i=Kf4i=  0 

In  matrix  form,  the  above  can  be  restated  by 

[Kf]=  Md  [N]  (4.22) 

The  sign  of  Kfn  will  be  dependent  upon  the  sign  of  UBFi. 
This  is  because  the  direction  of  the  incremental  friction 
will  be  different  from  that  of  the  initial  friction  when 
the  unbalanced  force  is  not  compatible  with  the  slip 
degree-of-freedom.  If  this  happens  the  incremental 
displacement  will  cause  the  total  normal  force  and  the 
total  friction  to  be  reduced  and  finally  the  stiffness  will 
be  decreased.  Similarly,  the  sign  of  Kf2i  will  be  opposite 
to  that  of  UBF2.  When  rli=l,  the  segment  II-III  will  be  in 
compression  and  will  decrease  the  total  normal  force  and 
frictional  force  at  node  III.  Therefore,  the  direction  of 
the  incremental  friction  will  be  opposite  to  that  of  the 
initial  friction  at  node  III.  The  derivation  of  slip 
stiffness  matrix  for  rli=l  is  shown  in  Fig.  4.10. 

When  rl2=l  and  the  other  displacements  are  zero,  the 
force  of  each  segment  is  given  by 

Si=  0 

S2=  K2 
S3=-K3 
84=85=0 


(4.23) 
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Fig.  4.10  Derivation  of  Slip  Stiffness  Matrix  (rll 
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From  the  equlibrium  of  each  node  in  both  tangential  and 
normal  directions,  the  terms  of  {Kn)2  and  {N}2  are  given  by 
Kni2=-K2 
Kn22=  K2+K3 

Kn32=-K3  (4.24) 

Kn42=  0 

and 

Ni2=  sin  a 

N22=  (K2-K3)  sin  a (4.25) 

N32=  K3  sin  b 
N42=  0 

The  friction  terms  can  be  evaluated 
Kfi2=±  Md  K2  sin  a 
Kf22=±  Md  (K2-K3)  sin  a 

Kf32=T  Md  K3  sin  b (4.26) 

Kf42=  0 

where  b=  90  - $2/2. 

Again,  the  sign  of  Kf22  is  the  same  as  that  of  UBF2 
and  the  sign  of  Kfi2  is  compatible  to  that  of  UBFi.  But  the 
sign  of  Kf32  is  opposite  to  that  of  UBF3  (Fig.  4.11).  The 
reasons  are  as  stated  earlier.  Similarly,  the  other 
stiffnesses  can  be  derived.  The  slip  stiffness  matrix  is 
expressed  as  follows: 

[Kl]  = [Kn]  + [Kf] 


(4.27) 
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Fig.  4.11  Derivation  of  Slip  Stiffness  Matrix  (rl2  =1) 
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(4.28) 


0 

0 


(4.29) 


! 0 tK3  sinb  ±(K3-K4)sinb  ±K4  sinb 

0 0 TK4  sinb  ±(K4-K5)sinb 

We  have  to  notice  that  the  [Kl]  matrix  is  not 
symmetrical  if  the  profile  of  the  tendon  is  not  symmetrical 
or  the  loading  is  not  symmetrical. 

4. 3. 3. 8 Determination  of  internal  forces.  After  the 
unbalanced  force  matrix  and  the  slip  stiffness  matrix  are 
formed,  the  slip  displacement  (rl)  can  be  solved  by  using 


(rl}=  [Kl]"l  (UBF) 


(4.30) 


The  internal  force  of  each  segment  is  dependent  upon 
the  relative  slip  displacement  in  each  segment.  Therefore, 
the  force  of  each  segment  can  be  expressed  by 
Ql=  Si  + Kl  rli 


Q2=  ^2  + K2(rl2  -rli) 

Q3=  S3  + K3(rl3  -rl2)  (4.31) 

Q4=  S4  + K4(r4  -r3) 

QS=  S5  - K5  rl4 


4.4  Solution  Stratecfv 
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The  determination  of  the  tendon  forces  will  be 
strongly  related  to  the  sliding  conditions  of  the  nodes. 

Not  only  one  node  can  slide  but  at  most  four  nodes  can 
slide.  A comprehensive  method  is  proposed  to  solve  the 
problem.  The  procedures  can  be  stated  as  follows: 

1.  Assume  the  tendon  is  a combination  of  five  truss 
elements  and  apply  the  initial  forces  to  the  structure. 

2 . Find  the  displacements  at  the  nodes  and  the  force  of 
each  truss  element. 

3.  Find  the  sliding  condition  of  each  node  by  using  Eg.  4.4. 

4.  Decide  the  structural  model  of  the  tendon. 

5.  Solve  the  tendon  forces  by  using  the  method  stated 
earlier  and  calculate  the  resisting  internal  forces. 

6.  Reform  the  stiffness  matrix. 

7.  Resolve  the  structure. 

8.  Repeat  procedure  2-7  until  the  forces  are  balanced. 

4 . 5 Computer  Implementation 
4.5.1  Introduction 

The  tendon  element  with  friction  is  installed  in  ANSR. 
six  interface  subroutines  are  utilized  for  this  element. 

The  element  is  modified  from  the  frictionless  tendon 
element  which  is  included  in  Chapter  3 . The  programming 
maintains  the  fundamental  frame  of  the  original  version 
except  the  necessary  changes  for  the  specific  purpose.  In 
subroutine  INEL,  INIT  and  OUTS  only  a few  portions  are 
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changed  or  added.  But  in  subroutine  STIF,  STAT  and  RINT 
nearly  all  the  programming  is  rewritten  because  the 
approach  is  totally  different.  The  source  code  for  these 
subroutines  are  contained  in  Appendix  C.  The  functions  of 
each  subroutine  are  stated  as  follows: 

4.5.2  Subroutine  INELll 

The  subroutine  is  called  first  when  the  elment  is 
being  used.  The  purposes  of  the  subroutine  include  the 
following: 

First,  input  and  initialize  the  element  group  data  and 
store  the  informations  in  labeled  common  block  /INFGR/.  The 
variables  are  defined: 

COMMON  /INFGR/  NGR,NELS,MFST, IGRHED(IO) ,NINFC, LSTAT, 

LSTF , LSTC , NDOF , DKO , DKT , EPROP (4,250) 
in  which,  the  following  variables  are  changed: 


LSTAT 

= 170 

LSTF 

= 648 

NINFC 

= 818 

The  others  are  the  same  as  those  described  in  Chapter  3 . 

Next,  read  the  element  data  and  store  the  informations 
in  the  common  block  /INFEL/.  The  labeled  common  block 
/INFEL/  for  this  element  includes  the  following: 

COMMON  /INFEL/  IMEM, KST, LM { 18 ) , NODE ( 6 ) , INS ( 6 ) , INSO ( 6) , XX ( 6 ) , 

YY(6) ,ZZ(6) ,PROP(3) ,AREA,XL(5) ,B(7) ,SEL(5) , 
STOT(5) ,TT(3,3) ,T,QL(12) ,SKP(324) 
in  which  IMEM,  NODE (6) , AREA  and  T are  read  from  input,  and 
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the  others  are  calculated  in  the  program.  The  definitions  of 
new  variables  are  stated  as  follows; 


INS (6) 

= Code  of  node's  sliding. 
= 0 for  no  sliding. 

= 1 for  sliding. 

= 0 for  initialization. 

INS0(6) 

= Previous  code  of  node's  sliding. 
= 0 for  initialization. 

PROP(l) 

= Young's  modulus  of  the  tendon. 

PROP (2) 

= Static  friction  coefficient  between  tendon  and 
concrete. 

PROP (3) 

= Dynamic  friction  coefficient  between  tendon  and 
concrete. 

XL(5) 

= Length  of  each  segment  in  the  tendon. 

SEL(5) 

= Current  tendon  force  of  each  segment. 

ST0T(5) 

= Current  tendon  stress  of  each  segment. 

QL(12) 

= Current  displacements  in  local  coordinate  system 

The  third  step  is  to  initialize  the  global  degree  of 
freedom  numbers.  This  is  done  by  calling  subroutine  NCODLM 
eighteen  times.  Then  the  values  are  stored  in  LM  array. 

The  fourth  step  is  to  transform  the  global  coordinates 
of  the  nodes  of  the  tendon  into  the  local  coordinates  by 
calling  subroutine  TRANl.  And  then  store  the  local 
coordinates  in  array  XX(6) , YY(6)  and  ZZ(6).  The 
transformation  is  discussed  in  Chapter  3. 

The  fifth  step  is  to  evaluate  the  direction  parameters 
and  store  in  an  array  B. 


68 


The  next  step  is  to  initialize  array  INS(6),  INS0(6) 
and  QL(12) . From  the  solution  strategy,  the  element  is 
assumed  to  be  a combination  of  five  truss  elements  at  the 
beginning.  Therefore,  the  array  INS  and  INSO  are  set  to 
zero  in  this  stage.  Finally,  the  stiffness  matrix  profile 
is  determined  by  calling  the  subroutine  BAND  and  the 
element  data  is  transferred  to  tape  by  calling  the 
subroutine  COMPAC. 

4.5.3  Subroutine  STIFll 

This  subroutine  is  used  to  form  the  element  stiffness. 
Because  of  friction,  the  stiffness  of  the  tendon  will  not 
be  unique,  however,  will  be  strongly  dependent  upon  the 
nodes'  sliding  conditions  during  the  iterative  analysis.  The 
nodes'  sliding  conditions  dominate  the  structural  model,  and 
furthermore  affect  the  element  stiffness.  Therefore,  the 
structural  model  should  be  determined  before  the  stiffness 
matrix  is  formed.  The  functions  of  this  subroutine  include 
the  following: 

First,  initialize  the  local  stiffness  matrix 
ST(12xl2) . 

Second,  compute  the  matrix  ST  according  to  the  nodes' 
sliding  conditions.  There  are  sixteen  combinations  for  the 
structural  model.  The  flow  chart  for  deciding  the  structural 
model  and  the  stiffness  are  shown  in  Fig.  4.12.  The  element 
stiffness  is  actually  a combination  of  tendon  stiffness  and 
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Fig.  4.12  Flow  Chart  for  Determination  of  Structural  Model 


truss  stiffness.  The  formation  of  tendon  stiffness  was 


70 


provided  in  Chapter  3 and  will  not  be  stated  here.  Also, 
the  derivation  of  combination  of  tendon  and  truss  stiffness 
was  stated  earlier.  There  are  several  subroutines  set  up 
seperately  in  this  program  in  order  to  provide  the 
different  element  stiffness  for  the  specific  structural 
model.  These  include  subroutine  TRU46,  TEN46,  TT36,  TEN36 
and  TEN16.  The  purpose  of  subroutine  TRU46  is  to  form  the 
stiffness  for  the  truss  structure  from  node  JJJ  to  node  J. 
Subroutine  TEN46  is  provided  to  form  the  stiffness  for  the 
tendon  element  from  node  JJJ  to  node  J.  Subroutine  TT36  is 
used  for  the  combination  of  the  tendon  element,  from  node 
III  to  node  JJ,  and  the  truss  element,  from  node  JJ  to  node 
J . Subroutine  TEN36  is  used  for  the  tendon  element  from 
node  III  to  node  J.  And  subroutine  TEN16  is  used  for  the 
tendon  elewments,  from  node  I to  node  J,  or  from  node  I to 
node  JJ,  or  from  node  II  to  node  JJ,  or  from  node  II  to  J. 

Next,  transform  the  local  stiffness  matrix  ST,  twelve 
degrees  of  freedom,  into  the  global  stiffness  matrix  SK, 
eighteen  degrees  of  freedom.  The  transformation,  identical 
to  the  one  dicussed  before,  is  done  by  calling  subroutine 
TRAN2.  The  12  X 18  transformation  matrix  is  used  in  this 
calculation. 

Finally,  store  the  global  stiffness  matrix  with  the 
standard  form  (SKP)  in  the  element  common  block  INFEL. 


4.5.4  Subroutine  INITll 
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The  main  purpose  of  this  subroutine  is  to  apply  the 
prestressing  forces  to  the  structure.  It  is  called  at  the 
beginning  of  each  load  step.  The  functions  of  this 
subroutine  are  expressed  as  follows; 

The  first  step  is  to  calculate  the  loading  factor  fac 
which  is  the  fraction  of  the  initial  load  to  be  included  in 
the  current  load  increment.  Next,  multiply  the  factor  fac 
to  initial  force  T to  get  current  initial  force  matrix  P 
and  determine  the  equvalent  loading  by  calling  subroutine 
FORVER.  Subroutine  FORVER  has  been  changed  to  adapt  the 
friction  version.  Because  the  tendon  force  might  not  be  the 
same  along  the  tendon  due  to  friction.  Therefore,  P matrix 
was  used  instead  of  a constant  value.  The  local  force 
matrix  is  expressed  by 


r ^ 

I -P(l) 

I P(l)  - P(2)*B(3) 
P(2)*B(1) 
P(2)*B(3)  - P(3) 
-P(2)*B(1) 

P(3)  - P(4)*B(4) 
[F]l  = -P(4)*B(2) 

-P(4)*B(5) 

: P(4)*B(4)  - P(5) 

I P(4)*B(2) 
P(4)*B(5) 

P(5) 


(4.32) 


in  which  B matrix  was  defined  earlier. 

The  final  step  is  to  transfoirm  the  local  force  matrix 
into  the  global  force  matrix  by  calling  subroutine  TRAN3.  The 
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transformation  matrix  was  stated  earlier  and  the  formulation 
was  based  upon  Eq.  (3.11)  in  Chapter  3. 

4.5.5  Subroutine  STATU 

The  subroutine  provides  current  state  information 
which  dominates  the  change  of  the  element  stiffness.  This 
can  be  done  by  comparing  the  current  nodes'  sliding 
conditions  with  the  previous  nodes'  sliding  conditions.  If 
there  are  no  changes  in  sliding  conditions  then  the  state 
code  KST  will  be  set  to  zero,  otherwise  KST  will  be  set  to 
one.  But  before  the  comparisons,  the  new  silding  conditions 
should  be  determined  first.  The  procedures  for 
determination  of  sliding  conditions  basically  follows  the 
solution  strategy  discussed  earlier.  It's  purposes  include 
the  following: 

First,  it  sets  the  previous  sliding  code  INSO  to  the 
current  code  INS  before  recalculates  the  current  code. 

Second,  it  transforms  the  global  nodal  displacements  Q 
into  local  nodal  displacements  and  updates  the  current 
local  nodal  displacements  QL.  The  transformation,  based 
upon  Eq.  (3.2),  is  done  by  calling  subroutine  TRAN4. 

Third,  it  calcualtes  the  additional  deformation  TD 
(not  including  the  deformation  due  to  initial  force)  of 
each  truss  element  using  the  relative  nodal  displacements. 

Fourth,  it  determines  the  loading  factor  fac  used  in 
the  calculation  of  the  total  resisting  internal  forces.  Fac 
is  the  total  fraction  of  the  initial  load  that  has  been 
applied. 
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Fifth,  it  initializes  the  force  of  each  tendon  segment 
SEL  which  is  set  to  the  corresponding  truss  force  TRUSF. 

The  truss  force  is  calculated  by  adding  the  initial  force  T 
to  the  nodal  force  induced  by  the  additional  deformation 
TD.  The  step  matches  the  assumption  that  the  tendon  is 
supposed  to  be  five  truss  elements  at  the  beginning  of  the 
solution  process. 

Sixth,  it  initializes  the  current  sliding  index  array 

INS. 

Seventh,  it  determines  the  sliding  condition  of  each 
node  by  comparing  the  driving  force  DF  with  the  maximum 
friction  force  FF.  If  the  node  slides,  calculates  the 
unbalnced  force.  Otherwise,  sets  the  unbalanced  force  to 
zero.  Then  the  structural  model  can  be  decided  according  to 
the  nodes'  sliding  conditions. 

Eighth,  it  forms  the  slip  stiffness  matrix  SF  and 
solves  the  slip  displacements  R.  The  slip  stiffness  matrix 
SF  includes  the  non-frictional  part  AK,  independent  of  the 
direction  of  the  unbalanced  force,  and  the  frictional  part 
BK,  dependent  upon  the  direction  of  the  unbalanced  force. 
The  slip  displacements  are  solved  by  calling  subroutine 
SOL.  The  major  task  in  subroutine  SOL  is  to  determine  the 
slip  degrees  of  freedom  of  the  tendon  first  and  then 
retrieve  the  corresponding  slip  stiffness  matrix  and  the 
unbalanced  force  matrix  and  finally  call  equation  solver 
subroutine  GROUT  to  get  displacements. 
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Ninth,  it  calculates  the  tendon  force  according  to 
Eq.  (4.31).  Next,  recheck  the  unsliding  nodes  by  comparing 
the  updated  sliding  conditions  INS  and  the  sliding 
conditions  of  the  last  iteration  IN.  If  the  two  sliding 
indices  are  different,  then  repeat  the  procedures  from 
seventh  to  ninth  until  they  are  the  same. 

Finally,  set  state  code  KST  by  comparing  the  sliding 
indices  INS  and  INO. 

4.5.6  Subroutine  RINTll 

The  subroutine  evaluates  the  element  force  recovery. 

It  is  called  in  the  beginning  of  the  analysis  and  after 
each  state  determination.  Because  the  force  of  each  segment 
SEL  is  already  determined  from  subroutine  STATU,  the  force 
does  not  need  to  be  calculated  here.  The  purposes  of  this 
subroutine  are 

The  first  step  is  to  calculate  stress  STOT  of  each 
segment  of  the  element. 

Next  step  is  to  form  the  resisting  internal  force 
matrix  FT  by  calling  subroutine  FORVER.  And  then,  transform 
FT  in  local  axes  into  FE  in  global  axes  by  calling 
subroutine  TRAN3 . 

The  final  step  is  to  set  damping  force  matrix  FD.  The 
damping  force  is  not  considered  in  the  tendon  element,  so 
the  force  matrix  is  set  to  zero. 


CHAPTER  5 

TENSIONING  PROCESS 
5. 1 Introduction 

For  the  pretensioned  member,  the  initial  force  along 
the  whole  length  can  be  regarded  as  a constant.  But  for  the 
post-tensioned  member,  the  initial  force  will  not  be 
constant  because  of  friction  and  anchorage  slip.  The 
instaneous  losses  should  be  accounted  for  in  the  initial 
force. 

Friction  loss  is  developed  between  the  steel  and  the 
duct  when  the  steel  slides  through  the  duct.  The  total 
friction  loss  is  the  sum  of  the  wobble  friction  due  to  the 
unintentional  misalignment,  and  the  curvature  friction  due 
to  the  intentional  curvature  of  the  tendon  [15].  For 
external  prestressing,  the  wobble  friction  does  not  exist 
because  the  tendon  is  outside  of  the  concrete  and  no 
reaction  is  provided  for  the  friction.  Therefore,  only  the 
curvature  friction  is  developed  in  the  external  tendons. 
Tensioning  at  one  end  is  most  commonly  used  for  external 
prestressing  but  tensioned  at  both  two  ends  is  also 
included  in  this  chapter. 

Anchorage  slip  will  happen  in  pretensioned  members  as 
well  as  post-tensioned  members.  But  slip  loss  is 
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insignificant  in  beams,  pretensioned  by  long-line  method, 
due  to  the  great  length  of  the  tendon  and  often  is  ignored. 
The  magnitude  depends  upon  the  particular  prestressing 
system  or  hardware  used  and  is  available  from  the 
manufacturer  or  from  laboratory  tests.  Usually,  the  loss 
is  more  important  in  short  tendons. 

The  following  sections  will  provide  the  method  to 
calculate  the  immediate  losses. 

5.2  Friction  Loss  During  Tensioning 
5.2.1  Tensioning  at  One  End 

If  the  tendon  is  stretched  at  one  end  and  anchored  at 
the  other  end,  the  tendon  force  will  gradually  decrease 
with  the  distance  from  the  point  to  the  tensioned  end.  The 
curvature  friction  loss  can  be  predicted  by  using  Eg.  (2.2). 
However,  the  direction  of  the  friction,  described  in 
Chapter  4,  is  assumed  to  bisect  the  angle  between  two 
segments.  According  to  this  assumption,  the  equation  can  be 
simplified. 

Suppose  the  tensioned  end  is  at  node  I,  the  tendon 
force  of  each  segment,  shown  in  Fig.  5.1,  is  expressed  by 
Tfif  Tf2,  Tf3,  Tf4  and  Tfs.  Let  f^,  f2,  f3  and  f4  designate 
the  friction  at  node  II,  III,  JJJ  and  JJ,  respectively.  If 
the  jacking  force  is  T,  then 


Tfi  = T 

Tf2  = T - fi 

Tf3  = T - fi  - f2 


(5.1) 


77 


Fig.  5.1  Friction  Loss  During  Tensioning  at  One  End 
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Tf4  = T - fi  - f2  - f3 

Tf5  = T - fi  - f2  - f3  - f4 

The  friction  loss  at  each  node  is  derived  by  the  following: 
At  node  II,  the  normal  force,  N,  is  (shown  in  Fig.  5.2) 

N = (Tfi+Tf2)sin  a 

= (2T  - fi)sin  a (5.2) 

where  a = (90-^i/2) 

From  the  classical  Coulomb's  friction  law,  the  friction, 
fl,  eguals  to  the  normal  force,  N,  times  the  curvature 


friction  coefficient, 

fl  = ^^di  N 

(5.3) 

Substitute  N,  from  Eq.  (5.2),  into  Eq.  (5.3) 

fl  = A*di  (2T  - fi)sin  a 

(5.4) 

Rearrange  the  above  equation,  put  f^  on  the 

same  side 

We  get 

fl(l+Mdisin  a)  = 2MdiTsin  a 

(5.5) 

Solve  fl 

fl  = 2^diTsin  a/A 

(5.6) 

where  A = 1 + /j^isin  a 

Similarly,  the  friction  at  node  III,  jjj  and 

JJ  can  be 

expressed  by 

f2  = 2Mdi(T-fi) sin  a/A 

(5.7) 

f3  = 2^di (T-fi“f2) sin  b/B 

(5.8) 

f4  = 2/idi  (T-fi-f2“f3)  sin  b/B 

(5.9) 

where  b = 90  - ^2/2  and  B = l + ^disin  b 

After  determiming  the  friction  loss  at  each  node,  the 

tendon  force  can  be  calculated  by  using  Eq.  (5.1). 
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Fig.  5.2  Derivation  of  Friction  Loss  at  Node  II 


5.2.2  Tensioning  at  Both  Ends 
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If  the  tendon  profile  is  syimetrical  and  stretched  at 
both  ends,  the  tendon  force  Tfs  and  tf4  will  equal  to  Tfi 
and  Tf2,  respectively.  And  Eq.  (5.1)  can  be  adjusted  to  the 
following; 

Tfl  = Tfs  = T 

Tf2  = Tf4  = T - ffi  (5.10) 

Tf3  = T -ffi  -f2 

The  magnitude  of  f^  and  f2  are  described  in  Eq.  (5.6)  and 
Eq.  (5.7).  Fig.  5.3  shows  the  tendon  force  with  tensioning  at 
both  ends. 

5 . 3 Anchorage  Slip 
5.3.1  Tensioning  at  One  End 

For  external  prestressing,  Eq.  (2.3-4)  can  not  be 
applied  because  of  the  high  friction  at  the  harping  point 
and  no  wobble  friction.  The  slip  loss  will  be  resisted  by 
the  reverse  friction  at  each  node.  The  length  affected  by 
the  loss  will  depend  upon  the  reverse  friction  at  each 
node.  If  the  friction  is  high  enough  at  the  particular 
node,  the  loss  will  be  stopped  at  that  node.  Otherwise,  the 
effect  will  extend  to  the  next  node  and  may  affect  the 
whole  tendon. 

Let  Tsi,  Ts2/  Ts3,  Tg4  and  Tg5  be  the  tendon  force 
after  the  anchorage  loss.  Given  the  length  of  slip  loss  aL, 
the  loss  of  the  tendon  force  at  segment  I-II,  Afanc/ 
expressed  by 


can  be 
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Fig.  5.3  Friction  Loss  During  Tensioning  at  Both  Ends 
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(5.11) 


'^^anc  “ ^p/^1 

where  Li  is  length  of  the  segment  I-II. 

The  tendon  force  of  segment  I-II  is  then 

Tsl  = Tfi  - Afanc  (5.12) 

The  procedures  to  determine  the  intial  force  after 
slip  loss  are  explained  as  follows  (shown  in  Fig.  5.4); 

5. 3. 1.1  Determine  the  sliding  condition  at  node  II. 
According  to  the  sliding  criterion  described  in  Chapter  4, 
compare  Tgi  and  Tf2  and  determine  the  sliding  condition.  If 
the  node  II  does  not  slide,  the  tendon  force  will  remain 
the  same  and  no  adjustment  is  necessary.  On  the  contrary, 
if  node  II  slides,  then  the  tendon  forces  of  segment  I-II 
and  II-III  should  be  updated. 

The  procedures  are  similar  to  those  described  in 
Chapter  4 except  for  the  slip  stiffness  matrix.  The  truss 
stiffness  [Kn]  stays  the  same  but  the  signs  in  the  friction 
matrix  [Kf]  will  be  somewhat  different,  because  the  effect 
of  the  slip  loss  will  decrease  with  the  distance  from  the 
tensioning  end.  Therefore,  the  direction  of  the  unbalanced 
force  in  this  case  will  be  fixed,  the  sign  in  the  matrix 
[Kf]  will  be  constant.  To  distinguish  the  difference,  the 
friction  matrix,  [Kf]i, 


is  restated 


Tsl  = Tfl  - fane 
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Fig.  5.4  Anchorage  Slip  Loss  (Tensioning  at  One  End) 
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[Kf]i  = Mdi 


(Kl-K2)sln  a K2  sin  a 
-K2  sin  a (K2-K3)sin  a 
0 -K3  sin  b 

0 0 


1 

0 0 ! 

! 

K3  sin  a 0 ! 

CK3-K4)sin  b K4  sin  b 
-K4  sin  b (K4-K5)sin  b 

i 


(5.13) 


Treat  Tgi  and  Tf2  as  fixed  forces  at  segments  I-li  and 
II-III,  respectively.  Using  Eg.  (4.3-5)  and  replacing  /ig  and 
Md  with  and  l^i,  calculate  the  unbalanced  force,  UBFij. 
And  then  solve  rli  by  using  Eg.  (4.30)  and  replacing  [Kf] 
with  [KfJi.  Notice  that  there  is  only  one  degree-of-freedom 
in  this  case.  From  Eg.  (4.31),  the  tendon  force  Tgi  and  Tg2 
are  updated  by 


Tsl  = Tgi  + K1  rli 

Ts2  = Tf2  - K2  rli  (5.14) 

5 *3. 1.2  Determine  the  sliding  condition  at  node  ttt.  Compare 
Tg2  and  Tf3.  If  node  III  does  not  slide,  no  force 
adjustment  is  necessary.  Otherwise,  this  is  a two  degrees- 
of-freedom  distribution.  Only  node  III  has  the  unbalanced 
force,  UBFiij.  The  unbalanced  force  at  node  II,  UBFji,  is 
zero.  After  forming  the  unbalanced  force  matrix,  solve  rli 
and  rl2  by  using  Eg.  (4.30).  The  updated  tendon  forces  are 
Tsl  = Tgi  + K1  rli 

'^S2  = Tg2  + K2  (rl2  - rli)  (5.15) 

Ts3  = Tf3  - K3  rl2 

5 >2. 1.2  Determine  the  sliding  condition  at  node  JJJ.  Compare 
Tg3  and  Tf4.  Again,  if  no  sliding  happens  at  node  JJJ,  then 
no  refinement  is  made.  Otherwise,  calculate  the  unbalanced 
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force  at  node  JJJ.  The  unbalanced  forces  at  node  II  and  III 
are  zero.  After  forming  the  unbalanced  force  matrix  and 
slip  stifness  matrix,  solve  rl^,  rl2  and  rl3  by  using 
Eq.  (4.30).  The  updated  tendon  forces  are 
Tsl  = Tgi  + K1  rli 
Ts2  = Ts2  + K2  (rl2  - rli) 

Ts3  = Ts3  + K3  (rl3  - rl2)  (5.16) 

Ts4  = Tf4  - K4  rl3 

5 . 3 . 1 . 4 Determine  the  sliding  condition  at  node  JJ.  Compare 
the  updated  force  Tg4  and  Tf5.  If  no  sliding  happens,  the 
effect  of  the  slip  loss  reaches  to  node  JJJ  only  and  no 
refinement  is  necessary.  If  node  JJ  slides,  then  this  is  a 
four  degrees-of-freedom  case.  Again,  only  node  JJ  has  the 
unbalanced  force.  The  slip  displacement  matrix,  (rl),  can 
be  solved  by  using  Eq.  (4.30).  The  initial  forces  after  slip 
loss  are 

Tsl  = Tgi  + K1  rli 

Ts2  = Tg2  + K2  (rl2  - rli) 

Ts3  = Tg3  + K3  (rl3  - rl2)  (5.17) 

Ts4  = Tg4  + K4  (rl4  - rl3) 

■^sS  = Tf5  - K5  rl4 

5.3.2  Tensioning  at  Both  Ends 

Suppose  the  tendon  is  tensioned  at  both  ends,  then  the 
tendon  force  Tgs  and  Tg4  will  equal  to  Tgi  and  Tg2, 
respectively.  The  same  criterion  is  used  here.  But  at  most 
two  nodes  need  to  be  checked  (shown  in  Fig.  5.5). 


Tsl  = Tfl  - fane 
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Fig.  5.5  Anchorage  Slip  Loss  (Tensioning  at  Both  Ends) 
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5.4  Computer  Implementation 
The  tensioning  process  is  added  in  subroutine  INEL. 
After  the  instaneous  losses  are  calculated,  the  initial 
force  will  be  stored  and  retrieved  in  subroutine  INIT.  For 
the  particular  applications,  including  the  tendon  stiffness 
in  the  analysis  will  be  the  user's  option.  The  procedures 
are  done  in  subroutines  STAT  and  STIF.  The  changes  to  the 
previous  element  routines  are  stated  as  follows: 

5.4.1  Subroutine  INELll 

The  labeled  common  block  /INFEL/  has  been  expanded; 


C0MM0N/INFEL/IMEM,KST,LM(18) ,N0DE(6) ,INS(6) ,INS0(6) ,XX(6) , 

* YY(6) ,ZZ(6) ,PR0P(5) , AREA, XL ( 5) , B (7 ) ,SEL(5) , 

* ST0T(5) ,TT(3,3) ,T,TI(5) ,QL(12) ,SHA1,SHA2, 

* BK(4,4) ,SK(5) ,IP,IDUM,ANSL,SKP(324) 

where  PR0P(5),  TI(5),  SHAl,  SHA2 , BK(4,4),  SK(5) , IP,  IDUM 


and  ANSL  are  added.  The  definitions  are  stated  as  follows: 

PROP(l)  = Young's  modulus  of  the  tendon. 

PROP (2)  = Static  friction  coefficient  between  the  tendon 

and  the  concrete  during  tensioning. 

PR0P(3)  = Dynamic  friction  coefficient  between  the  tendon 

and  the  concrete  during  tensioning. 

PROP (4)  = Static  friction  coefficient  between  the  tendon 

and  the  concrete. 

PROP (5)  = Dynamic  friction  coefficient  between  the  tendon 

and  the  concrete. 

TI(5)  = The  initial  force  of  each  segment  after 

friction  loss  and  anchorage  slip  loss. 


SHAl  = Sin(90-^i/2) 

SHA2  = Sin(90-(?2/2) 

BK(4,4)  = The  truss  stiffness  matrix.  One  part  of  the 

slip  stiffness  matrix. 

SK(5)  = The  truss  stiffness  of  each  segment. 

SHAl,  SHA2 , BK(4,4)  and  SK(5)  were  calculated  in  subroutine 
STAT  and  are  shifted  here  to  avoid  repeat  calculations. 

IP  = The  loading  segment  number  to  include  the 

tendon  stiffness  in  the  analysis. 

This  means  the  analysis  will  start  including  the  tendon 
stiffness  from  the  ip  th  loading  segment.  This  number  can 
be  greater  than  the  total  loading  segments  and  no  tendon 
stiffness  will  be  included  in  the  analysis. 

IDUM  = Dummy  interger  variable. 

ANSL  = The  amount  of  the  anchorage  slip. 

The  variables  for  the  tensioning  process  are  input 
after  the  material  properties.  These  include  NT,  IP  and 
ANSL,  where  NT  is  the  index  for  the  tensioning  at  one  end 
or  both  two  ends.  The  losses  will  be  calculated  after  array 
INS(6),  INSO(6)  and  QL(12)  are  initialized.  The  friction 
loss  is  calculated  first,  the  procedures  are 
(1)  Tell  the  number  of  ends  at  which  the  tendon  is 
stretched  during  tensioning.  If  the  tendon  is 
tensioning  at  both  ends,  check  the  tendon  profile's 
symmetry. 
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(2)  Distinguish  the  tendon's  profile,  4-node  or  5-node  or 
6-node. 

(3)  Apply  Eq.  (5.6  -10)  to  calculate  the  tendon  forces  and 
store  them  in  TI(5).  If  the  tendon  is  not  6-node,  the 
force  adjustments  are  made. 

After  the  friction  loss  has  been  calculated,  the 
anchorage  slip  loss  is  then  evaluated.  The  step-by-step 
method  has  been  explained  in  Chapter  4 and  section  5.3.  The 
procedures  used  in  the  programming  are 

(1)  Calculate  the  stiffness  of  each  truss  segment,  SK(5) , 
and  initialize  the  loss,  TLOSS,  due  to  the  slip. 

(2)  Form  the  slip  stiffness  matrix,  SF. 

(3)  The  loss,  TLOSS,  changes  the  tendon  force  of  the 
segment  I-II,  TI(1),  and  induce  the  force 
redistribution.  First,  check  the  sliding  condition  of 
node  II.  If  it  slides,  then  check  node  III.  The  check 
will  continue  until  the  non-sliding  node  or  the  last 
node  JJ  is  reached. 

(4)  Each  sliding  condition  will  induce  the  redistribution 
of  the  tendon  force.  The  tendon  forces  will  be 
determined  based  upon  Eq.  (5.14-17). 

5.4.2  Subroutine  INiTii 

The  common  block  /INFEL/  is  changed  to  the  one 
discussed  in  the  above  section.  The  variable,  NSEG, 
indicating  the  number  of  the  loading  segment,  is  added  to 
the  common  block  /STLDPT/. 
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The  other  change  in  this  subroutine  is  to  replace  the 
the  force,  T,  with  the  initial  force  after  losses,  TI,  in 
calculating  the  current  force  matrix,  P.  Next,  the 
equivalent  loading  is  determined  by  calling  subroutine 
FORVER.  The  initial  force  matrix  is  shown  in  Eg.  (4.25). 

The  magnitude  and  the  direction  of  the  loading  is  shown  in 
Fig.  5.6. 

5.4.3  Subroutine  STIFll 

The  tendon  stiffness  is  excluded  when  the  index,  IP, 
is  smaller  than  the  current  loading  segment,  NSEG.  Even  if 
the  sliding  code,  INS(5),  has  been  changed,  and  the 
subroutine  is  called  to  reform  the  stiffness,  the  tendon 
stiffness  will  stay  zero  under  this  situation.  These 
statements  are  done  at  the  beginning  of  the  subroutine. 

5.4.4  Subroutine  STATll 

Except  that  the  common  blocks  /INFEL/  and  /STLDPT/  are 
changed,  the  common  block  /PRITC/  is  added  into  this 
subroutine.  The  only  variable  used  in  this  common  block  is 
the  iteration  number,  NIT. 

In  the  force  solving  process,  the  force,  T,  is 
replaced  by  the  corresponding  initial  force,  TIi,  in  the 
ith  segment. 

When  the  index,  IP,  equals  to  the  current  loading 
segment,  the  tendon  stiffness  should  be  added  to  the  total 
stiffness.  At  the  beginning  of  the  loading  segment,  i.e., 
the  first  iteration,  the  state  code,  KST,  is  automatically 
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set  to  one.  These  statements  are  put  in  the  end  of  the 
subroutine. 
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Fig.  5.6  Initial  Force  After  Immediate  Losses 


CHAPTER  6 

EXAMPLES  AND  PARAMETRIC  STUDY 
6 . 1 Introduction 

Four  examples  are  presented  in  this  chapter  to 
demonstrate  the  validity  and  applicability  of  the  tendon 
element  developed  here.  In  order  to  find  out  the  influence 
of  the  friction  on  the  structural  behavior,  a parametric 
study  regarding  the  cofficient  of  friction,  the  tendon 
force  and  external  force  is  presented. 

The  first  example  is  used  to  verify  the  tendon  element 
without  friction.  The  vertical  displacement  and  the  change 
of  the  tendon  force  are  compared  to  those  derived  from 
beam  theory.  This  example  can  be  applied  for  the 
pretensioned  prestressed  concrete  structure  because  no  slip 
between  the  tendon  and  the  concrete  is  allowed.  The  second 
example  shows  the  friction  loss  and  the  seating  loss  during 
the  tensioning  process  and  the  friction  at  the  deviation 
block  during  the  application  of  external  loading.  This 
example  shows  the  use  of  external  prestressing  for  the 
post-tensioned  prestressed  concrete  structure.  The  third 
example  shows  the  application  for  the  continuous  beam  with 
external  prestressing.  The  results  are  then  compared  to 
those  from  conventional  analysis.  The  final  example 
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introduces  the  link  element  [22,23]  to  model  the  interface 
of  the  concrete  segments  in  the  segmental  bridge  with 
external  prestressing. 

The  parametric  study  presents  the  effects  of  the 
different  parameters  to  the  whole  structural  behavior.  The 
vertical  displacements,  the  tendon  forces  and  the  gap 
widths  of  the  segmental  bridge  due  to  the  variations  of  the 
parameters  are  presented. 

6.2  Example  1 — For  the  Pretensioned  Bonded  Structure 

A simple-supported  beam  with  the  harped  tendon  is 
adopted  in  this  example.  The  configuration  of  the  tendon 
and  the  dimensions  of  the  beam  are  shown  in  Fig.  6.1.  The 
material  properties  used  for  the  analysis  are  shown  in 
Table  6.1.  The  prestressed  concrete  beam  is  modeled  by 
forty-eight  shell  elements  [22,23]  and  one  tendon  element 
in  the  xy  plane.  To  verify  the  finite  element  mesh  and  to 
account  for  the  contribution  of  each  loading,  three  loading 
segments  are  performed.  The  first  loading  segment  is  due  to 
the  prestressing  force,  the  second  one  is  due  to  the  self- 
weight of  the  beam  and  the  third  one  is  due  to  one  pair  of 
concentrated  loads  at  one-third  of  the  span  from  each 
support . 

Two  kinds  of  the  boundary  conditions  are  used  in  this 
example,  the  first  is  used  to  match  the  assumption  of  beam 
theory  and  the  end  supports  are  put  along  the  neutral  axis 
of  the  cross  section,  the  other  model's  end  supports  are 


Fig.  6.1  Pretensioned  Prestressed  Concrete  Beam — Example 
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put  at  the  bottom  of  the  cross  section.  In  the  first 
model,  two  cases  are  performed;  one  with  the  tendon 
stiffness  and  the  other  without  the  tendon  stiffness,  and 
the  poisson's  ratio  is  set  to  zero.  In  the  second  model, 
the  tendon  stiffness  is  considered  and  the  poisson's  ratio 
is  set  to  0.2. 

6.2.1  Verification  of  Displacements 

The  theoretical  value  of  the  displacement  at  midspan 
due  to  each  loading  case  is  calculated  from  beam  theory 
and  can  be  stated  as  follows; 

6. 2. 1.1  Camber  due  to  the  prestressing  force.  For  the 
harped  tendon  used  here,  the  camber  is  [15,16] 

Ai  = 23/216 (TeL2/EI)  (6.1) 

in  which 

T = prestressing  force. 

e = eccentricity  of  the  tendon  to  the  neutral  axis. 

L = span  length. 

E = Young's  modulus  of  concrete. 

I = moment  inertial  of  the  net  cross  section 
(w/o  tendon  stiffness) . 

= moment  inertial  of  the  transformed  section 
(w  tendon  stiffness) . 

6. 2. 1.2  Displacement  due  to  the  self-weiaht  of  concrete. 

The  displacement  due  to  the  distributed  self-weight  is 

Aw  = 5/384  (wLVeI)  (6.2) 

in  which  w is  the  concrete  weight  per  unit  length. 


6 . 2 . 1 . 3 Displacement  due  to  the  concentrated  loads . The 
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displacement  due  to  the  point  loads  (shovm  in  Fig.  6.1)  is 
expressed  by 

Ap  = 23/648 (PL3/EI)  (6.3) 

6.2.2  Verification  of  the  Change  of  the  Tendon  Force 

(With  Tendon  Stiffness) 

If  the  prestressed  concrete  beam  is  regarded  as  the 
composite  structure  of  concrete  and  steel  tendon,  then  the 
increment  of  the  tendon  force  can  be  derived  from  the 
virtual  work  [24] 

^le  + Tine  ^11  = 0 (6.4) 

in  which 

Ale  = the  displacement  along  the  tendon  due  to  the 
external  loads. 

Tine  “ the  increment  of  the  tendon  force. 

All  = the  displacement  along  the  tendon  due  to  the  unit 
tendon  force. 

The  equivalent  loading  of  the  tendon  can  be  divided  into 

two  parts;  one  is  flextural  loads  and  the  other  is  axial 

loads  (shown  in  Fig.  6.2).  Then 

All  = 2/EI[  (5/162)  sin2crL3 
+ (1-cosa) 2e2L/6 
+sina(l-cosa) eL2/9] 

+L/EA(2/3*COSa+l/3) +L/EgAs ( l/3+COs2a*2/3 ) (6.5) 

where  a is  the  sharp  angle  at  the  harped  point. 

For  the  self-weight  of  the  beam,  the  displacement  Ai^  is 
Aiw  = WL^/EI [-11/486  sina  L - 13/324 (1-COSa) e]  (6.6) 
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Fig.  6.2  Equivalent  Loading  of  Unit  Tendon  Force — Example 


For  the  point  loads,  the  displacement  a^p  is 

Alp  = (PL2/EI) *[-5/81  sina  L - (l-COSa)e/9]  (6.7) 

Substitute  Eq.  (6.6-8)  into  Eq.  (6.5),  the  change  of  the 
tendon  force  at  each  stage  can  be  evaluated. 

The  results  from  the  output  and  beam  theory  for  each 
case  are  listed  in  Table  6.2.  Compare  model  1 without 
tendon  stiffness  to  the  ideal  beam,  the  differences  of  the 
results  are  within  0.5%.  For  model  1 with  tendon  stiffness, 
the  increment  of  the  tendon  force  and  the  displacement, 
excluding  the  camber,  for  each  loading  are  very  close  to 
those  from  the  theory.  The  tendon  stiffness  has  more 
influence  in  the  camber  due  to  the  prestress,  because  the 
high  axial  loads  are  applied  at  the  ends  where  the  tendon 
stiffness  is  comparatively  high.  To  demonstrate  this,  the 
equivalent  loadings  without  the  axial  loads,  matched  to  the 
ideal  beam,  are  applied  to  the  same  beam  and  results  in  a 
higher  camber.  The  tendon  stiffness  reduces  the  vertical 
displacement,  except  for  the  camber,  about  1%.  These 
results  show  the  validity  of  this  element.  Model  2 
considers  the  effects  of  the  tendon  stiffness  and  the 
poisson's  ratio  and  the  results  are  satisfactory  compared 
to  the  ideal  beam. 

6 • 3 Example  2— For  External  Prestressina 
For  simplicity,  the  simple-supported  beam  with  a 
rectangular  cross  section,  shown  in  Fig.  6.3,  is  used  in 
this  example.  The  structure  is  a post-tensioned  prestressed 
concrete  beam  and  the  harped  tendon  is  unbonded  with  the 


Table  6.1  Material  Properties  of  Example  1 
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f'c(ksi)  Ec(ksi)  Wc(k/ft-^)  v Es(ksi)  As(in^) 

Fi(k) 

5.  4833  150  0 30000  2.5 

(0.2) 

300 

* v=0.2  for  model  2 

Table  6.2 

Results 

of  Example 

1 

Load* 

Model  1 

Ideal 

Model  1 

Composite 

Model  2 

Case 

(w/o  stiff) 

Beam 

(w  stiff) 

Beam 

1 

0.5263 

0.5236 

0.4818 

0.5093 

0.5016 

2 

-0.2350 

-0.23373 

-0.2317 

-0.22735 

-0.2335 

Y Dis. 

1+2 

0.2913 

0.28987 

0.2501 

0.28195 

0.2681 

3 

-0.5475 

-0.54455 

-0.5392 

-0.52968 

-0.5434 

1+2+3 

-0.2562 

-0.25468 

-0.2891 

-0.24773 

-0.2753 

Initial  force  = 

300. kips 

1 

272.8 

274.6 

274.5 

Tendon 

51 

27.2 

- 

25.4 

25 . 5 

Force 

1+2 

274.4 

- 

276.5 

— 

276.6 

52 

1.6 

- 

1.9 

1.96 

2 . 1 

1+2+3 

279.3 

- 

281.0 

— 

281.5 

53 

4.9 

— 

4.5 

4.57 

4.9 

P = 175  K 1 P = 175  K 

t t Anchorage  End 
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Fig.  6.3  External  Prestressing  Beam — Example 
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concrete  except  at  the  deviation  blocks.  The  material 
properties  are  shown  in  Table  6.3.  Different  from  the 
pretensioned  beam,  there  are  friction  and  anchorage  slip 
losses  during  the  tensioning.  Friction  at  the  deviation 
blocks  will  happen  after  the  structure  is  deformed.  The 
tendon  stiffness  is  included  in  the  analysis,  because  the 
tendon  contacts  with  concrete  at  the  deviation  blocks. 

Three  cases  are  performed  in  this  example.  The  first 
case  is  to  show  the  validity  of  the  instaneous  prestress 
losses  due  to  friction  and  anchorage  slip.  For  post- 
tensioned  prestressed  concrete  structures,  jacking  force  at 
the  tensioning  end  is  input  instead  of  the  initial  force. 
After  the  instaneous  losses  happened,  the  net  tendon  force 
is  then  applied  to  the  structure.  The  friction,  due  to  the 
deformation  of  the  deviation  block,  is  ignored  by  setting 
coefficient  of  friction  to  zero.  Therefore,  the  tendon  is 
continuous  and  the  final  tendon  force  is  uniform  along  the 
length.  The  output  is  shown  in  Table  6.4.  The  verification 
of  the  friction  loss  can  be  performed  by  using  Eg.  (2.2). 
The  tendon  forces  after  the  immediate  losses  are  shown  in 
Fig.  6.4. 

The  second  case  shows  the  application  of  the 
friction  induced  by  the  deformation  of  the  structure.  The 
initial  losses  are  neglected,  i.e.,  the  initial  force  is 
assumed  constant  along  the  length.  Two  loading  segments  are 
used,  the  first  one  is  the  prestress  and  the  self-weight  of 


Table  6.3  Material  Properties  of  Example  2 
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Case  1 

Case  2 

■kifk 

Case  3 *** 

fc(ksi) 

6. 

Ec(ksi) 

4415 

— 

Wc(lb/ft3)  150 

- 

— 

V 

0.2 

— 

As(in2) 

23.256 

- 

Eg(ksi) 

29000 

- 

Fi(kip) 

4200. 

N/A 

N/A 

Fj (kip) 

5000. 

4200. 

5000. 

Mi* 

0.2 

0.0 

0.2 

Mf** 

0.0 

0.3 

0.3 

*anc ( 

0.3 

N/A 

0.3 

- tne  same  as  case  1 

* Mi  - 

friction 

coefficient 

during  tensioning. 

**  Mf  - 

friction 

coefficient 

after 

force  transfer. 

***  the 

tendon  is  tensioned 

at  one 

end  only. 

Table  6.4  Output 

for  Case  1 & Case  3 of  Example  2 

Jacking  force  = 5000. 

k Mi  -0.2  5anc 

=0.3  in. 

Initial  Tendon 
II-III  III-JJJ 

Force 

JJJ-JJ 

After  friction  loss 
After  seating  loss 

5000.  4909. 

4677.  4764. 

4819. 

4819. 

Dis. 

Final  Tendon  Force 

Load  1+2  .1256 

Load  1+2+3  -.2853 

4504.  4627. 

4542.  4670. 

4638. 

4660. 
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Fig.  6.4  Immediate  Prestress  Losses — Example 
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the  beam  and  the  second  one  is  the  external  loads,  shown  in 
Fig.  6.3.  The  output,  shown  in  Table  6.5,  is  compared  to 
that  with  no  friction. 

The  third  case  shows  the  combined  effects  of 
instaneous  prestress  losses  and  the  friction  loss  due  to 
the  deformation  of  the  structure.  The  results  are  shown  in 
Table  6.4. 

6 . 4 Example  3 — For  the  Continuous  Girder 

This  example  shows  the  application  of  the  five-node 
tendon,  mostly  used  as  external  prestressing  in  continuous 
girders.  The  structure  has  two  equal  spans  and  is  post- 
tensioned  by  external  prestressing.  The  dimensions  of  the 
structure  are  shown  in  Fig.  6.5,  and  the  material  properties 
are  shown  in  Table  6.6.  Three  tendons  are  used  in  each  span 
and  overlapped  in  the  middle  pier  segment  to  assure  the 
continuity  of  the  structure.  The  configurations  of  the 
tendons  in  each  span  are  arranged  to  be  symmetrical  to  the 
middle  support.  Each  tendon  is  tensioned  at  one  end  only 
and  has  the  same  amount  of  anchorage  slip.  The  direction  of 
the  tensioning  is  also  shown  in  Fig.  6.5. 

Two  loading  segments  are  used  here.  The  first  is  the 
prestress  and  the  self-weight  of  the  girder,  the  second  is 
the  external  point  loads,  shown  in  Fig.  6.5.  The  results 

listed  in  Table  6.7.  For  the  purpose  of  comparison,  the 
same  structure  is  analyzed  again  by  the  conventional  method 
in  which  the  program,  STAN  [25],  is  used.  In  the 


Table  6.5  Output — Case  2 of  Example  2 
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/if=0.3  /if=0 

Tendon  Force  Tendon 


Load 

Dis. 

III-JJJ 

JJJ-JJ 

Dis. 

force 

1 

.8820 

4003. 

3972. 

4003. 

.8840 

3992. 

1+2 

.0083 

4049. 

4094. 

4049. 

.005625 

4064. 

1+2+3 

-.4092 

4070. 

4152. 

4070. 

-.4065 

4098. 

*load  1 — Prestress 
load  2 — Beam  Weight 
load  3 — External  loads 


Table  6.6  Materials  Properties  of  Example  3 
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fc' 

Wc 

Ec 

As  Eg 

Fj 

Mi 

Mf 

'^anc 

6. 

150. 

4415. 

5.814  29000 

1250 

0.2 

in 

• 

o 

0.3 

Table  6.7  Results  of  Example  3 


Load  Case  Displacement  (in.) 

ANSR  STAN 

1.  Prestress  . 5903  . 60554 

2.  Self-weight  -.49913  -.47749 

3.  External  Loads  -.1527  -.14845 

4.  Total  -.06153  -.02040 


*Displacement  was  measured  at  the  midspan  of  the  first  span. 
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Continuous  Beam — Example 
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conventional  model,  shown  in  Fig.  6.6,  the  beam  elements  are 
used  and  the  tendon  elements  are  replaced  by  the  equivalent 
loadings.  The  friction  loss  and  the  anchorage  slip  loss 
during  the  tensioning  are  calculated  due  to  the  assumption 
that  the  amount  of  the  anchorage  slip  is  distributed  along 
the  tendon.  No  friction  effects,  due  to  the  deformation  of 
the  structure,  and  no  tendon  stiffnesses  are  considered. 

The  vertical  displacements  from  the  conventional  analysis 
are  included  in  Table  6.7. 

6.5  Example  4 — For  the  Secrmental  Bridge 
This  example  has  the  same  dimensions  with  the  example 
used  in  section  6.3,  but  more  tendons,  with  different 
configurations,  are  used  to  control  the  gaps  between  the 
segments.  Link  elements  [22,23]  are  added  along  the 
interface  between  the  segments  in  order  to  simulate  the 
behavior  of  gap  openings.  The  dimensions  and  the  tendon 
PJ^ofiles  are  shown  in  Fig.  6.7.  The  material  properties  are 
listed  in  Table  6.8. 

Since  the  tendons  are  tensioned  at  both  ends,  the 
initial  forces  are  symmetrical.  Friction  at  the  interface 
of  the  segments,  modeled  by  link  elements,  is  included.  If 
the  link  elements  are  open,  the  structural  stiffness  will  be 
reduced.  The  reduction  of  the  stiffness  depends  upon  the 
number  of  open  link  elements.  To  study  the  effects  of  the 
friction  at  the  deviation  blocks,  two  cases  are 
investigated  in  this  example.  The  first  includes  the 
fi^iction  while  the  second  does  not. 
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Fig.  6.6  Structural  Model  Used  for  Compa r ison--Example 


See  Detail  A 
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Fig.  6.7  Example  for  Segmental  Bridge 
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The  output,  including  the  vertical  displacement,  the 
tendon  force  and  the  gap  width,  is  shown  in  Table  6.9.  The 
results  indicate  that  the  friction  reduces  not  only  the 
maxmium  gap  width  but  also  the  vertical  displacement  at  the 
midspan.  Therefore,  the  friction  can  reduce  the  possibility 
of  gap  openings  and  then  increase  the  ultimate  strength  of 
the  bridge. 

6 . 6 Effects  of  the  Friction 

Due  to  the  lack  of  any  test  data  about  the  friction 
coefficient  at  the  deviation  block  of  external 
prestressing,  the  study  uses  the  different  values  of  the 
friction  coefficient,  ranged  from  0.  to  0.5,  to  find  out 
the  effects  of  friction  on  structural  behavior.  The 
structures  of  example  2 and  4,  designated  regular  beam  and 
segmental  beam  respectively,  are  used  in  this  study. 

For  the  structure  of  example  2,  the  following 
parameters  are  used: 

1.  Initial  force  4200  kips  (immediate  losses  are 
ignored) . 

2.  External  load  from  0 to  500  kips. 

3.  Friction  coefficient  from  0.  to  0.3. 

Two  loading  segments  are  performed.  Fig.  6.8  shows  the 
load-deflection  curves  due  to  the  various  friction 
coefficients.  Fig.  6.9  shows  the  plots  of  the  tendon  forces 
at  midspan  as  functions  of  external  loads  due  to  the 


Table  6.8  Material  Properties  of  Example  4 
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Case  1 Case  2 


fc(ksi) 

6. 

- 

Ec(ksi) 

4415 

— 

Wc(lb/ft3) 

150 

- 

V 

0.2 

— 

As ( in2 ) 

7.752 

- 

Es(ksi) 

29000 

- 

Fj (kip) 

1650. 

- 

Mi* 

0.25 

- 

Mf** 

0.6 

0, 

Ms*** 

0.3 

- 

^anc(iri) 

0.3 

- 

- the  same  as  case  1 

* Mi  ” friction  coefficient  during  tensioning. 

**  Mf  - friction  coefficient  after  force  transfer. 

***Ms  - friction  coefficient  for  link  elements  (shear  transfer). 
****  the  tendon  is  tensioned  at  one  end  only. 


Table  6.9  Results  of  Example  4 


Mf Displacement  Max.  Gap  width  Tendon  Force 


0.  -0.5284  .001013 

0.6  -0.5232  .0005 


12  3 

1534  1550  1568 

1581  1571  1579 

(1519  1522  1546)* 


*the  tendon  force  at  the  tensioning  end. 


500 


(sdi>f)  pDo~i 


Displacement  at  Midspan  (in.) 


500 
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Fig.  6.9  Tendon  Force  at  Midspan  v.s.  External  Load  for 
Various  Friction  Coefficients  (Example  2) 


different  friction  coefficients.  The  plots  of  the  tendon^^^ 
forces  at  the  tensioning  end  versus  the  external  loads  due 
to  the  various  friction  coefficients  are  shown  in  Fig. 

6.10. 

the  structure  of  example  4,  the  following 
parameters  are  used: 

1.  Jacking  force  1250  kips  for  each  tendon. 

2.  The  tendons  are  tensioned  at  two  ends. 

3.  External  load  from  180  to  228  kips. 

4.  Friction  coefficient  from  0 to  0.5. 

Fig.  6.11  shows  the  load-deflection  curves  for  the 
various  friction  coefficients.  The  maximum  gap  width  verse 
the  external  load  for  the  various  friction  coefficients  is 
shown  in  Fig.  6.12. 

It  can  be  seen  that  the  displacement  decreases 
slightly  with  the  increase  of  the  friction  coefficient. 

This  can  be  explained  by  the  difference  between  the  tendon 
stiffness  and  the  truss  stiffness.  Furthermore,  the  open  of 
gaps  can  reduce  the  total  stiffnesses  of  the  structure,  so 
the  difference  is  more  obvious  in  example  4 than  in  example 
2.  When  the  structure  is  under  the  specific  exrternal  load 
and  the  friction  coefficient  is  big  enough,  then  all  the 
nodes  of  the  tendon  will  not  slide  and  the  tendon  behaves 
as  the  truss  elements.  At  this  time,  the  behavior  of  the 
tendon  remains  the  same  despite  the  increase  of  friction. 
Further,  the  increase  of  the  external  loads  will  reduce  the 
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Fig.  6.12  Maximum  Gap  Width  v.s.  External  Load  for 
Various  Friction  Coefficients 
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resistant  forces  at  the  deviation  blocks  and  cause  the 
nodes  to  slide,  and  the  tendon  behaves  as  a continuous 
element.  So  the  slope  of  the  load-deflection  curve  will  be 
slightly  decreased  due  to  the  increase  of  the  external 
loads.  This  tendency  is  more  significant  in  example  4 than 
in  example  2 due  to  the  gap  openings.  Under  the  same 
external  loads,  the  tendon  force  at  midspan  increases  with 
the  friction  coefficient  but  at  the  same  time  the  tendon 
force  at  the  tensioning  end  decreases  with  the  friction 
coefficient.  This  indicates  that  the  friction  increases  the 
strength  at  the  midspan  but  reduces  the  strength  at  the 
harped  point.  The  sudden  drop  of  the  tendon  force  should  be 
noticed  especially  the  friction  is  big  enough.  From  Fig. 
6.12,  it  can  be  deary  seen  that  the  friction  helps  reduce 
the  maximum  gap  width  by  about  100%  of  that  with  no 
friction. 

The  results  point  to  the  importance  of  the  friction  at 
the  deviation  blocks  for  the  prestressed  non-segmental  and 
segmental  bridges. 
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CHAPTER  7 

CONCLUSIONS  AND  RECOMMENDATIONS 

7 . 1 Conclusions 
7.1.1  Tendon  Element  Without  Friction 

The  tendon  element  without  friction  has  been  verified 
and  proved  to  be  useful  in  modeling  the  prestressing 
reinforcement  for  both  the  pretensioned  and  the  post- 
tensioned  prestressed  concrete  structures. 

The  formation  of  the  tendon  stiffness  is  based  upon 
the  continuity  along  the  whole  tendon  which  is  the  major 
difference  with  the  formation  of  truss  element  or  frame 
element.  Consequently,  this  is  a more  reasonable  approach 
than  the  conventional  method.  Further,  considering  the  tendon 
stiffness  is  more  plausible  in  the  analysis  of  the  bonded 
prestressed  concrete  structures. 

The  initial  force  of  the  tendon  is  automatically 
transfered  to  the  structure  in  the  form  of  the  equivalent 
loading.  This  is  another  advantage  of  the  element,  it 
avoids  having  to  calculate  the  equivalent  loading, 
especially  when  many  tendons  are  used.  For  the  tendon 
element  without  friction,  the  initial  force  is  constant. 

This  matches  the  behavior  of  the  pretensioned  concrete 
beams  and  is  still  applicable  to  the  bonded  post-tensioned 
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structures  with  the  proper  estimate  for  the  immediate 
losses . 
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The  prestress  loss  due  to  the  deformation  of  the 
structure  is  evaluated  from  the  change  in  the  geometry  of 
the  tendon.  Since  the  angle  change  is  very  small,  the 
calculation  is  according  to  the  change  of  the  total  tendon 
length.  The  final  tendon  force  is  then  calculated  by 
deducting  the  loss  from  the  initial  force. 

The  applicability  of  the  element  has  been  verified  in 
Chapter  6.  The  displacement,  due  to  the  prestress,  and  the 
increment  of  the  tendon  force,  due  to  external  loading, 
proved  to  be  satisfactory  for  the  analysis  of  prestressed 
concrete  structures. 

7.1.2  Tendon  Element  With  Friction 

The  element  was  developed  to  simulate  the  external 
prestressing  reinforment  used  in  post-tensioned  prestressed 
concrete  structures.  Immediate  prestress  losses,  such  as 
friction  loss  and  anchorage  loss,  are  included. 

The  friction  loss  is  much  less  in  external  pre- 
stresssing than  that  in  internal  prestressing  because  of 
the  less  contact  area  between  steel  and  concrete.  The 
simplied  formula  presented  here  proved  to  be  adequate 
compared  to  the  more  frequently  used  formula.  The  anchorage 
loss  can  be  predicted  by  using  the  step-by-step  method 
proposed  in  Chapter  5.  The  reverse  friction,  caused  by  the 
seating  loss,  is  solved  by  the  slip  stiffness  method 
proposed  in  this  study. 
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The  formation  of  the  stiffness  matrix  is  strongly 
dependent  upon  the  structural  model  of  the  prestressing 
steel.  Because  of  the  existence  of  friction,  the  tendon 
element  may  not  be  one  continuous  element  but  a combination 
of  truss  elements  and  tendon  elements.  The  structural  model 
is  determined  by  the  sliding  condition  at  each  node.  Since 
the  element  is  nolinear,  the  stiffness  should  be  updated  if 
the  structural  model  is  changed. 

The  initial  force  is  not  constant  along  the  tendon  due 
to  the  instaneous  prestress  losses.  After  the  immediate 
losses  have  been  calculated,  the  equivalent  loading  is  then 
formed  and  transfered  to  the  concrete. 

The  slip  stiffness  was  developed  to  solve  the  tendon 
force  with  regard  to  the  friction  at  the  deviation  block. 
The  formation  is  basically  the  combination  of  the  truss 
stiffness  for  the  segmental  members  of  the  tendon  and  the 
friction  stiffness  at  the  deviation  block.  The  degrees—of— 
freedom  of  the  element  are  determined  by  the  structural 
model.  The  unbalanced  load  is  determined  by  the  equilibrium 
at  each  degree-of-freedom.  The  final  tendon  force  is  then 
solved  by  the  force  distribution  inside  the  tendon  element. 

Although  the  friction  at  the  deviation  block  has  not 
been  investigated,  the  method  used  here  is  an  available 
approach  to  verify  the  importance  of  the  friction.  The 
results  presented  in  Chapter  6 have  shown  the 


effects  of  the  friction  at  the  deviation  blocks,  especially 
for  the  segmental  bridges.  It  can  be  concluded  that  enough 
friction  at  the  deviation  block  increases  the  tendon 
stiffness  and  help  reduce  the  maximum  gap  width  and 
increases  the  strength  of  the  segmental  bridge. 

7.2  Recommendations 

The  following  recommendations  are  made  based  upon  the 
results  of  this  study: 

1.  The  tendon  element  without  friction  can  provide 
good  results  for  the  analysis  of  the  prestressed 
concrete  structure  with  bonded  tendons. 

2.  External  tendon  at  the  deviation  block  of  the 
segmental  bridges  should  be  well  bonded  to  provide 
enough  force  to  resist  gap  opening. 

3 . The  external  tendon  used  in  the  non— segmental 
structure  can  be  unbonded  at  the  deviation  block. 
But  the  concrete  block  should  be  designed  to  resist 
the  high  tension  of  the  tendon. 

4.  The  force  jump  of  the  external  tendon  at  the 
deviation  block  due  to  the  friction  should  be  taken 
into  account  in  design. 

5.  The  profile  of  the  external  tendon  is  better 
arranged  to  form  the  bigger  angle  at  the  harped 
point  to  provide  larger  friction.  But  the  high 
tension  at  the  block  should  be  considered. 

6.  The  location  of  the  deviation  block  should  be  well 


distributed  along  the  concrete  segments  in  order  to 
control  the  gap  opening  in  segmental  bridges. 

7 . 3 Future  Work 

Convergence  problems  have  been  encountered  when  the 
tendon  element  (with  friction)  and  link  elements  are  used 
at  the  same  time.  A possible  reason  is  due  to  the 
nonlinearities  of  both  elements,  the  nonlinear  analysis 
method  used  in  ANSR  is  not  able  to  find  the  right 
convergence  path.  The  improvement  of  the  convergence 
process  will  be  valuable  for  the  analysis  of  the  segmental 
bridges. 

The  long-time  prestress  losses  due  to  creep  and 
shinkage  of  concrete  and  relaxation  of  the  prestressing 
steel  are  not  included  here.  If  the  shell  element  has  the 
capability  of  handling  creep  and  shinkage,  then  the 
prestress  losses  can  be  done  without  any  change  in  this 
tendon  element.  If  the  convergence  can  be  improved,  the 
loss  due  to  relaxation  of  the  prestressing  steel  can  also 
be  added. 

The  structures  used  in  the  examples  were  modeled  by 
shell  elements  in  the  x-y  plane  only.  This  is  because  there 
was  difficulty  in  modeling  the  box  section  with  the 
existing  shell  elements.  The  displacement  from  the  model 
didn't  match  that  from  beam  theory.  Hopefully,  the 
improvement  of  shell  elements  can  help  make  box  section 
model  more  reasonable. 
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Due  to  the  lack  of  experimental  data  for  the  friction 
at  the  deviation  block,  the  method  proposed  here  can  not  be 
verified  and  compared.  If  experiments,  measuring  the 
friction  coefficient  and  the  variation  of  the  tendon  force, 
can  be  conducted,  it  will  help  understand  the  effects  of 
friction. 


APPENDIX  A 


USERS  MANUAL  FOR  THE  TENDON  ELEMENT 


TENDON  ELEMENT  USER'S  GUIDE 


A.l  CONTROL  INFORMATION 

- One  Card 

Columns 

Notes 

Name 

Data 

1-5(1) 

NGR 

Element  group  indicator. 

6-10(1) 

NLES 

Number  of  elements  in  group. 

11-15(1) 

MFST 

Element  number  of  first  element 

group . 

16-25 (F) 

DKO 

Initial  stiffness  damping 

factor,  Pq. 

26-35(F) 

DKT 

Current  tangent  stiffness 

damping  factor, 

41-80(A) 

GRHED 

Optional  group  heading. 

A.  2 MATERIAL  PROPERTY  INFORMATION 

A. 2 (a) 

Control 

Card 

Columns 

Notes 

Name 

Data 

1-5(1) 

NMAT 

Number  of  different  material 

types.  Default=l 


A.  2(b)  Properties  - NMAT  set  of  cards 


Columns  Notes 

Name 

Data 

1-5(1) 

M 

Material  Number. 

*6-15(F) 

E 

Young's  modulus  of  elasticity. 
Static  friction  coefficient 
between  duct  and  tendon. 

*16-25 (F) 

/^is 

*26-35(F) 

Mid 

Dynamic  friction  coefficient  between 
duct  and  tendon. 

*36-45(F) 

Mfs 

Static  friction  coefficient  between 
concrete  and  tendon. 

*46-55(F) 

Mfd 

Dynamic  friction  coefficient 
between  concrete  and  tendon. 

*A.2(c)  Tensioning  method 

& Anchorage  slip  length 

1-5(1) 

nt 

Index  for  tensioning  at  one  end 
or  both  two  ends. 

= 1 one  end 
= 2 two  ends 

6-10(1) 

ip 

The  load  segment  number  from 
which  the  analysis  starts 
including  the  tendon  stiffness. 

11-20 (F) 

ansi 

Anchorage  slip  length  due  to  the 

prestressing  method. 
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A. 3 ELMENT  DATA  - one  card 
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Columns  Notes 

Name 

1-5(1) 

N 

6-10(1) 

NODC 

11-15(1) 

NODC 

16-20(1) 

NODC 

21-25(1) 

NODC 

26-30(1) 

NODC 

31-35(1) 

NODC 

36-40(1) 

MT 

41-50(F) 

AO 

51-60(F) 

T 

62-66(1) 

ND 

Data 

Element  number 

Node  number  I 

Node  number  II 

Node  number  III 

Node  number  JJJ 

Node  number  JJ 

Node  number  J 

Material  type  number 

Area  of  tendon 

Initial  tension  of  tendon 

Node  number  differences 


*For  the  version  without  friction,  the  specified  items  or 
section  should  be  omitted. 


APPENDIX  B 


LISTINGS  OF  THE  TENDON  ELEMENT  WITHOUT  FRICTION 


c 


SUBROUTINE  INELll  (NJT,NDKOD, X, Y, Z , KEXEC) 


C 

C 

C 

C 

C 

C 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  TO  READ,  GENERATE  AND  PRINT  ELEMENT  DATA. 
ELEMENT  TYPE  = 11  THREE  DIMENSIONAL  TENDON  ELEMENTS. 

***************************************************** 

IMPLICIT  DOUBLE  PRECISION (A-H,0-Z) 

DIMENSION  NDKOD(NJT,3) ,X(NJT) ,Y(NJT) ,Z(NJT) ,ICOM(756) 
COMMON  /TAPES  / NIU, NOU, NTl , NT2 , NT3 , NT4 , NTS , NT6 , NT7 , NTS , 

* NT9 , NTIO , NTH , NT12 , NT13 , NT14 , NT15 , NT16 
COMMON  /INFGR  / NGR,NELS,MFST, IGRHED(IO) ,NINFC,LSTAT, 

* LSTF,LSTC,NDOF,DKO,DKT,EPROP(1000) 

COMMON  /INFEL  / IMEM, KST, LM ( 18 ) , NODE ( 6) , XX ( 6 ) , YY ( 6) , 

* ZZ(6) ,PROP(l) ,AREA,XL,B(7) ,SEL,VTOT, 

* STOT,TT(3,3) ,T,SKP(324) 

COMMON  /PRECIS/  ZERO, ONE 

COMMON  /WORK  / NODC (6) , XE ( 6) , YE ( 6) , ZE ( 6) 

EQUIVALENCE  ( IMEM , ICOM ( 1 ) ) 

DIMENSION  IAST(3) 

DATA  lAST  /2H  ,2H  *,2H**/ 


MEANINGS  OF  VARIABLES  IN  COMMON  BLOCK  /INFEL/ 


IMEM 
KST 
LM(18) 

NODE (6)  - 

KGEOM 

KTHO 

KBUCK 

KOD 

PROP(l)  - 

AREA 

XL 

XYZ(3,6)  - 
B(8) 

TT(3,3)  - 

VTOT 

STOT 

SEP 

SEL 

VENP,  VENN 


- ELEMENT  NUMBER. 

- STIFFNESS  UPDATE  CODE. 

- LOCATION  MATRIX. 

- NODE  NUMBERS  FROM  END  I TO  J. 

- LARGE  DISPLACEMENT  CODE. 

- RESPONSE  OUTPUT  CODE. 

- BUCKLING  BEHAVIOR  CODE. 

- CURRENT  YIELD  CODE. 

- MATERIAL  PROPERTIES. 

- CROSS  SECTIONAL  AREA. 

- ELEMENT  LENGTH. 

- COORDINATES  OF  NODE  I TO  J. 

- DIRECTION  PARAMETERS. 

- TRANSFORMATION  MATRIX. 

- TOTAL  AXIAL  DEFORMATION. 

- TOTAL  STRESS 

- AXIAL  FORCE  IN  ELASTO-PLASTIC  COMPONENT. 
AXIAL  FORCE  IN  ELASTIC  COMPONENT. 


- POS.  AND  NEG.  ENVELOPE  VALUES  OF  AXIAL 
DEFORMATION. 

VPACP,  VPACN 

- ACCUMULATED  POS.  AND  NEG.  PLASTIC  DEFORMATIONS. 
VBUCK  - BUCKLING  DEFORMATION. 

SENP,  SENN 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

100 

c 

c 

c 

c 

c 

c 


c 


110 

c 

c 

c 


c 

140 

c 
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- POS.  AND  NEG.  ENVELOPE  VALUES  OF  AXIAL  FORCE. 
TVENP,  TVENN 

- TIMES  AT  POS.  AND  NEG.  ENVELOPE  VALUES  OF 
DEFORMATION. 

TSENP,  TSENN 

- TIMES  AT  POS.  AND  NEG.  ENVELOPE  VALUES  OF  AXIAL 
FORCE . 

TOLY  - OVERSHOOT  TOLERANCE  FOR  YIELD  EVENTS 

TOLU  - OVERSHOOT  TOLERANCE  FOR  UNLOADING  EVENTS 

KODP  - PREVIOUS  YIELD  CODE. 

LDUM  - INTEGER  VARIABLE  NOT  USED 

SKP(324)  - CURRENT  TANGENT  OR  EFFECTIVE  STIFFNESS  MATRIX. 


ELEMENT  GROUP  PARAMETERS 

IF  (MFST.LE.O)  MFST  = 1 

NDOF  =18 
LstaT  = 108 
LSTF  = 648 
LSTC  = 0 

NINFC  = 756 
KST  = 1 

DO  100  1=3,756 
ICOM(I)  = 0 

PRINT  ELEMENT  GROUP  INFORMATION 

WRITE (NOU ,2000)  NGR, IGRHED , NELS , MFST , NDOF , NINFC , DKO , DKT 

READ  AND  PRINT  MATERIAL  PROPERTIES 

READ  (NIU,1000)  NMAT 
IF  (NMAT.LE.O)  NMAT  = 1 

WRITE (NOU, 2010)  NMAT 

DO  110  1=1, NMAT 

READ  (NIU,1005)  M,EPROP(M) 

WRITE (NOU, 2020)  (M,EPROP(M) ,M=1,NMAT) 

READ,  GENERATE  AND  PRINT  ELEMENT  DATA 

IMEM  = MFST 

NLAST  = MFST  + NELS  - 1 

WRITE (NOU, 2030) 

READ  (NIU,1010)  N,NODC,MT,AO,T,ND 

IF  (MT.LE.O)  MT  = 1 
IF  (ND.EQ.O)  ND  = 1 


O O O O to  o 
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C 

150 

C 

155 


C 

160 

170 

C 


C 


C 


C 

190 

200 


10 


220 

C 


IF  (N  - IMEM)  155,160,210 

JFAST  = IAST(3) 

KEROR  = 1 

WRITE (NOU, 2040)  JFAST, N,N0DC,MT, AO, T 
GO  TO  290 

DO  170  1=1,6 
NODE (I)  = NODC(I) 

KEROR  = 0 
NDIF  = ND 
MTYP  = MT 
AREA  = AO 
FORCE  = T 
SEL=t 

JFAST  = lAST(l) 

IF  (N.NE.NLAST)  GO  TO  140 
GO  TO  210 

DO  200  1=1,6 

NODE (I)  = NODE (I)  + NDIF 
JFAST  = IAST(2) 

WRITE ( NOU ,2050) JFAST , IMEM , NODE , MTYP , AREA , SEL 

ELEMENT  CONNECTIVITY  AND  DIRECTION 
TRANSFORMATION  OF  COORDINATES 

DO  220  1=1,6 
NOD  = NODE (I) 

II  = 3*1 

CALL  NCODLM(LM(II-2) ,NJT,NOD,l) 

CALL  NCODLM(LM(II-l) ,NJT,NOD,2) 

CALL  NCODLM(LM(II) ,NJT,NOD,3) 

XE(I)=X(NOD) 

YE (I) =Y (NOD) 

ZE(I)=Z(NOD) 

CALL  TRAN1(XE,YE,ZE,XX,YY,ZZ,TT) 

Z4=ZZ(5)-ZZ(4) 

DZ4=DABS (Z4) 

DX1=XX(2)-XX(1) 

IF(XX(1) .EQ.XX(2) ) DXl=ZERO 
DX2=XX(3)-XX(2) 

DX3=XX(4)-XX(3) 

DX4=XX(5)-XX(4) 

DX5=XX(6)-XX(5) 

IF(XX(6) .EQ.XX(5))  DX5=ZERO 
DY2=YY(2) -YY(3) 
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DY4=YY(5) -YY(4) 

XL1=DX1 

XL2=SQRT ( DX2  *DX2+DY2  *DY2 ) 

XL3=DX3 

XLH4  5=SQRT ( DX4  *DX4+DZ4  *DZ  4 ) 

XL4  5=XLH4  5 *XLH4  5+DY4  *DY4 
XL4=SQRT(XL45) 

XL5=DX5 

XL=XL1+XL2+XL3+XL4+XL5 
BX=EPROP (MT) *AO/XL 
CY1=DY2/XL2 
SY1=DX2/XL2 
CY2=DY4/XL4 
CX2=DX4/XL4 
CZ2=Z4/XL4 
SX11=1.  -SYl 
CX22=1.  -CX2 
C 

c INITIALIZE  DIRECTION  PARAMETERS 

C 

C 

B(1)=CY1 
B(2)=CY2 
B(3)=SY1 
B(4)=CX2 
B(5)=CZ2 
B(6)=SX11 
B(7)=CX22 
PROP(l)  =BX 
C 

C COMPUTE  STIFFNESS  MATRIX  PROFILE  AND  TRANSFER 
C ELEMENT  DATA  TO  TAPE 

CALL  BAND(LM,NDOF,NDKOD) 

CALL  COMPAC 
C 

C CHECK  LAST  ELEMENT 

C 

IF  (IMEM.EQ.NLAST)  GO  TO  290 
IMEM  = IMEM  + 1 
IF  (IMEM.EQ.N)  GO  TO  160 
GO  TO  190 
C 

C ERROR  AND  ELEMENT  DATA  GENERATION  MESSAGES 

290  IF  (KEROR.EQ.O)  GO  TO  300 
C 

WRITE (NOU, 2060) 

CALL  EXIT 
C 

300  RETURN 
C 

1000  FORMAT  (15) 
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1005 

1010 

2000 


2010 


2020 

2030 


2040 

2050 

2060 


C**** 

c 

c*  * * * 


100 


200 


FORMAT  (I5,F10.0) 

FORMAT  (8I5,2F10.0,1X,I5) 

FORMAT  (26H  ELEMENT  GROUP  INDICATOR  = 13, 

* 36H  (THREE  DIMENSIONAL  TENDON  ELEMENTS)//5X, 10A4// 

* 5X,43HNUMBER  OF  ELEMENTS  IN  THIS  GROUP = 15/ 

* 5X,43HNUMBER  OF  FIRST  ELEMENT  IN  THIS  GROUP  ...  = 15// 

* 5X,43HNUMBER  OF  DEGREES  OF  FREEDOM  PER  ELEMENT  = 15/ 

* 5X,43HLENGTH  OF  ELEMENT  INFORMATION  ARRAY  = 15// 

* 5X,43HDAMPING  COEFFICIENT,  BETA-0  = F11.5/ 

* 5X,43HDAMPING  COEFFICIENT,  BETA-T  = Fll!5) 

FORMAT  (//20H  MATERIAL  PROPERTIES// 

* 5X,36HNUMBER  OF  DIFFERENT  MATERIAL  TYPES  = 15// 

* 5X, 7HMAT.no. ,5X,7HMODELUS/) 

FORMAT  (I9,E12.5) 

FORMAT  (///20H  ELEMENT  INFORMATION// 

* 5X , 4HELEM , 5X , 4HNODE , 5X , 4HNODE , 5X , 4HNODE , 5X , 4HNODE , 

* 5X,4HNODE,5X,4HNODE,5X,4HMAT. , 8X, 4HAREA, 5X, 7HINITIAL/ 

* 5X,4H  NO.,5X,4H  I,5X,4H  II,5X,4H  III,5X,4H  JJJ, 

* 5X,4H  JJ,5X,4H  J, 5X, 4H  NO. , 17X, 7H  FORCE  /) 

FORMAT  (1X,A2,I6,7I9,2F12.4) 

FORMAT  (1X,A2,I6,7I9,2F12.4) 

FORMAT  (///15H  ERROR  MESSAGES// 

* 5X,  3HKEY,  5X,  29H MESSAGE , 5X,  5HERROR// 

* 5X,3H  , 5X, 29HELEMENT  DATA  INPUT  DIRECTLY  ,5X,5HNONE  / 

* 5X,3H  *,5X,29HELEMENT  DATA  GENERATED  ,5X,5HNONE  / 

* 5X,3H  **,5X,29HELEMENT  CARD  OUT  OF  SEQUENCE  , 5X, 5HFATAL/) 

END 

******************************************************** 
SUBROUTINE  TRAN 1- TRANS FORMAT I ON  of  COORDINATES  * 

********************************************^^^^^^^^^^^^ 
SUBROUTINE  TRANl (X, Y, Z , XX, YY, ZZ , T) 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  x(6) ,y(6) ,2(6) ,T(3,3) ,XX(6) ,YY(6) 

* ,22(6) ,V(4,4) 

CALL  VECTOR(V(l,l) ,X(3) ,Y(3) ,Z(3) ,X(4) ,Y(4) , 

* Z(4)  ) 

CALL  VECTOR(V(l,4) ,X(3) ,Y(3) ,Z(3) ,X(1) ,Y(1) ,Z(1)) 

CALL  CROS(V(l,l),V(l,4),V(l,3)) 

CALL  CROS(V(l,3) ,V(1,1) ,V(1,2) ) 

DO  100  K=l,3 
DO  100  L=l,3 
T(K,L)=V(L,K) 

DO  200  N=l,6 

XX(N)=T(1,1) *X(N)+T(1,2) *Y(N)+T(1,3)*Z(N) 

YY(N)=T(2,1) *X(N)+T(2,2) *Y(N)+T(2,3)*Z(N) 

ZZ (N)=T(3,1) *X(N)+T(3,2) *Y(N)+T(3, 3) *Z(N) 

CONTINUE 

RETURN 

END 

VECTOR 

SUBROUTINE  VECTOR  (V, XI , YI , ZI , XJ, YJ, ZJ) 

IMPLICIT  REAL*8  (A-H,0-Z) 
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DIMENSION  V(4) 

X=XJ-XI 

Y=yj-YI 

Z=ZJ-ZI 

V(4)=SQRT(X*X+Y*Y+Z*Z) 

V(3)=Z/V(4) 

V(2)=Y/V(4) 

V(1)=X/V(4) 

RETURN 

END 

* CROSS 

SUBROUTINE  CROS(A,B,C) 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  A(4) ,B(4) ,C(4) 

COMMON  /lOLIST/  NTM,NTR, NIN, NOT, NTl , NFL, NT2 ,NT3 , NT4 ,NT5 


X=A(2)*B(3)-A(3)*B(2) 

Y=A(3)*B(1)-A(1)*B(3) 

Z=A(1)*B(2)-A(2)*B(1) 

XYZ  = X*X  + Y*Y  + Z*Z 
IF(XYZ.GT.O.O)  GO  TO  100 
WRITE  (NTM,2000) 

WRITE  (NOT, 2000) 

STOP 

100  C(4)=SQRT(XYZ) 

C(3)=Z/C(4) 

C(2)=Y/C(4) 

C(1)=X/C(4) 

RETURN 

C 

2000  FORMAT  ( ' **  K NODE  IS  ALONG  ELEMENT  SEE  MANUAL  **'/) 
END 



SUBROUTINE  STIFll (ISTEP, NDF, CDKO, CDKT, FK, INDFK, ISTFC) 

c ******************************************************* 

C SUBROUTINE  TO  COMPUTE  ELEMENT  TANGENT  STIFFNESS  MATRIX. 

C ELEMENT  TYPE  = 11  THREE  DIMENSIONAL  TENDON  ELEMENTS. 

C ******************************************************* 

c 

IMPLICIT  DOUBLE  PRECISION (A-H, 0-Z) 

DIMENSION  FK (NDF, NDF) 

COMMON  /INFGR  / NGR, NELS , MFST, IGRHED ( 10) , NINFC, LSTAT, 

* LSTF , LSTC , NDOF , DKO , DKT , EPROP (1000) 
COMMON  /INFEL  / IMEM, KST, LM ( 18 ) , NODE ( 6) , XX ( 6) , YY ( 6) , 

* ZZ(6) ,PROP(l) , AREA, XL, B (7) ,SEL,VTOT, 

* STOT,TT(3,3) ,T,SKP(324) 

COMMON  /WORK  / S ( 12 , 12 ) , SK ( 18 , 18 ) 

COMMON  /PRECIS/  ZERO, ONE 

c set  to  full  matrix  version 

indfk=0 


C 


137 

C LINEAR  PART  OF  EFFECTIVE  STIFFNESS 

C 

C TENDON ( I, J)  OR  TENDON ( I, JJ)  OR  TENDON (I I, JJ)  OR 

C TENDON ( II, J) 

DO  10  1=1,12 
DO  10  J=I,12 
10  S(I,J)=ZERO 

S(l,l)=ONE 
S(l,2)=-B(6) 

S(l,3)=-B(l) 

S(l,4)=-S(l,2) 

S(l,5)=-S(l,3) 

S(l,6)=-B(7) 

S(1,7)=B(2) 

S(1,8)=B(5) 

S(l,9)=-S(l,6) 

S(l,10)=-S(l,7) 

S(l,ll)=-S(l,8) 

S(l,12)=-S(l,l) 

S(2,2)=B(6) *B(6) 

S(2,3)=B(6) *B(1) 

S(2,4)=-S(2,2) 

S(2,5)=-S(2,3) 

S(2,6)=B(6)*B(7) 

S(2,7)=-B{6) *B(2) 

S(2,8)=-B(6) *B(5) 

S(2,9)=-S(2,6) 

S(2,10)=-S(2,7) 

S(2,ll)=-S(2,8) 

S(2,12)=B(6) 

S(3,3)=B(1)*B(1) 

S(3,4)=-B(l) *B(6) 

S(3,5)=-S(3,3) 

S(3,6)=B(1)*B(7) 

S(3,7)=-B(l) *B(2) 

S(3,8)=-B(l) *B(5) 

S(3,9)=-S(3,6) 

S(3,10)=-S(3,7) 

S(3,ll)=-S(3,8) 

S(3,12)=B(1) 

DO  100  1=4,12 
100  S(4,I)=-S(2,I) 

DO  200  1=5,12 
200  S(5,I)=-S(3,I) 

S(6,6)=B(7) *B{7) 

S(6,7)=-B(7) *B(2) 

S(6,8)=-B(7) *B(5) 

S(6,9)=-S(6,6) 

S(6,10)=-S(6,7) 

S(6,ll)=-S(6,8) 

S(6,12)=B(7) 

S(7,7)=B(2)*B(2) 
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S(7,8)=B(2)*B(5) 

S(7,9)=-S(6,7) 

S(7,10)=-S{7,7) 

S(7,ll)=-S(7,8) 

S(7,12)=-B(2) 

S(8,8)=B(5) *B(5) 

S(8,9)=-S(6,8) 

S(8,10)=-S(7,8) 

S(8,ll)=-S(8,8) 

S(8,12)=-B(5) 

DO  300  1=9,12 
300  S(9,I)=-S(6,I) 

DO  400  1=10,12 
400  S(10,I)=-S(7,I) 

S(ll,ll)=-S(8,ll) 

S(ll,12)=-S(8,12) 

S(12,12)=ONE 

IF(XX(1) .EQ.XX(2)  ) THEN 

C FOR  TENDON ( II -JJ,J) 

DO  500  1=1,12 
500  S(l,I)=ZERO 

S(2,2)=B(3) *B(3) 

S(2,3)=-B(3)*B(1) 

S(2,4)=B(3)*B(6) 

S(2,5)=-S(2,3) 

S(2,6)=-B(3)*B(7) 

S(2,7)=B(3)*B(2) 

S(2,8)=B(3) *B(5) 

S(2,9)=-S(2,6) 

S(2,10)=-S(2,7) 

S(2,ll)=-S(2,8) 

S(2,12)=-B(3) 

ENDIF 

C FOR  TENDON ( I, II-JJ) 

550  IF(XX(5) .EQ.XX(6) ) THEN 
DO  600  1=1,12 
600  S(I,12)=ZERO 
S(l,9)=-B(4) 

S(2,9)=B(4)*B(6) 

S(3,9)=B(4)*B(1) 

S(4,9)=-S(2,9) 

S(5,9)=-S(3,9) 

S(6,9)=B(4) *B(7) 

S(7,9)=-B(4)*B(2) 

S(8,9)=-B(4)*B(5) 

S(9,9)=B(4) *B(4) 

S(9,10)=-S(7,9) 

S(9,ll)=-S(8,9) 

ENDIF 

C for  tendon (II-JJ) 

650  IF(XX(1) ,EQ.XX(2) ,AND.XX(5) .EQ.XX(6) ) THEN 
S(l,9)=ZERO 
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S(2,9)=-B(4) *B(3) 

ENDIF 

DO  700  K=l,12 

DO  700  L=l,12 

S(L,K)=S(K,L) 

700  CONTINUE 

DO  750  K=l,12 
DO  750  L=l,12 

S(L,K)=S(L,K)*PROP(l) 

750  CONTINUE 

CALL  TRAN2(TT,S,SK,18,1) 

IJ  =0 

DO  800  I=l,NDOF 
do  800  j=l,ndof 
IJ  = IJ  + 1 

STIF  = SK(I,J) 

FK  (I,J)  = STIF 

c FK  (J,I)  = STIF 

800  SKP(IJ)  = STIF 

C 

RETURN 
END 

C TRAN2 

SUBROUTINE  TRAN2 (T, SL, S,ND,NL) 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  T(3,3) ,TG(12,18) , SL( 12 , 12 ) , S (18 , 18 ) ,ST(12,18) 
COMMON/ PRECIS/  ZERO, ONE 
DO  50  1=1,12 
DO  50  J=l,18 
ST(I,J)=ZERO 
50  TG(I,J)=ZERO 

DO  60  1=1,18 
DO  60  J=l,18 
60  S(I,J)=ZERO 

c form  the  transformation  matrix  TG 

DO  100  J=l,3 
100  TG(1,J)=T(1,J) 

DO  150  1=2,3 
DO  150  J=4,6 
11=1-1 
JJ=J-3 

150  TG(I,J)=T(II,JJ) 

DO  200  1=4,5 
DO  200  J=7,9 
11=1-3 
JJ=J-6 

200  TG(I, J)=T(II, JJ) 

DO  250  L=6,9,3 
DO  250  I=L,L+2 
DO  250  J=L+4,L+6 
II=I-L+1 
JJ=J-L-3 


o n o o o o 
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250  TG(I,J)=T(II, JJ) 

DO  300  J=16,18 
JJ=J-15 

300  TG(12, J)=T(1, JJ) 

C EVALUATION  ST=SL*TG 

DO  400  1=1,12 
DO  400  J=l,18 
DO  400  K=l,12 

400  ST(I,J)=ST(I, J)+SL(I,K) *TG(K, J) 

C S=TTG*SL*TG 

DO  500  1=1,18 
DO  500  J=I,18 
DO  500  K=l,12 

500  S(I,J)=S(I,J)+TG(K,I) *ST(K,J) 

DO  600  1=1,18 
DO  600  J=I,18 
600  S(J,I)=S(I, J) 

RETURN 

END 



SUBROUTINE  STATU  (NDF, Q, TIME, FACAL, FACIL, ALFA) 



C **************************************************** 

C SUBROUTINE  FOR  STATE  DETERMINATION  CALCULATIONS. 

C ELEMENT  TYPE  = 11.  THREE  DIMENSIONAL  TENDON  ELEMENTS. 

C **************************************************** 

C 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

DIMENSION  Q(NDF) ,FACAL(1) ,FACIL(1) 

COMMON  /INFGR  / NGR, NELS ,MFST, IGRHED (10) , NINFC, LSTAT, 

* LSTF , LSTC , NDOF , DKO , DKT , EPROP (1000) 

COMMON  /INFEL  / IMEM, KST, LM(18) ,NODE ( 6) , XX(6) , YY (6) , 

* ZZ(6) ,PROP(l) ,AREA,XL,B(7) ,SEL,VTOT, 

* STOT,TT(3,3) ,T,SKP(324) 

COMMON  /PRECIS/  ZERO, ONE 

C 

KST  = 0 

RETURN 

END 



SUBROUTINE  RINTll  (NDF, Q, VEL, FE, FD, TIME, FACAL, FACIL, ALFA) 

**************************************************** 
SUBROUTINE  FOR  INTERNAL  FORCE  CALCULATIONS, 

ELEMENT  TYPE  = 11  THREE  DIMENSIONAL  TENDON  ELEMENTS. 
*************************************************** 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

DIMENSION  Q(NDF) , VEL (NDF) ,FE(NDF) ,FD(NDF) ,FT(12) , 

* FACAL ( 1) , FACIL (1) 

COMMON  /INFGR  / NGR, NELS , MFST, IGRHED ( 10) , NINFC, LSTAT, 

* LSTF , LSTC , NDOF , DKO , DKT , EPROP (1000) 


non  non  oooo 
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COMMON  /INFEL  / IMEM, KST, LM(18 ) , NODE (6) , XX (6) , YY (6) , 

* ZZ (6) ,PROP(l) , AREA, XL, B (7) , SEL, VTOT, STOT, 

* TT(3,3) ,T,SKP(324) 

COMMON  /STLDPT/  NPNF,NPFF,NPTP,NPP4 ,NPDP,NETP,NESW,NITSR, 

* NSPATT,NSEG 
COMMON  /PRECIS/  ZERO, ONE 

CALCULATE  THE  INTERNAL  RESISTING  FORCE 

IF(LM(1) .EQ. 0.and.LM(2) . EQ. 0 . and. LM(3 ) .EQ.O)  THEN 
IF(XX(1) .EQ.XX(2) ) THEN 
NA=4 
NB=5 
NC=6 
ELSE 
NA=1 
NB=2 
NC=3 
ENDIF 
SA=ZERO 
SB=ZERO 
SC=ZERO 

CALCUALTE  THE  TENDON  LOSS  DUE  TO  DEFORMATION 

DO  400  K=l,NDOF 

NAK=(NA-1) *NDOF+K 
NBK=(NB-1) *NDOF+K 
NCK=(NC-1) *NDOF+K 
sa=sa+skp (nak) *q(k) 
sb=sb+skp(nbk) *q(k) 
sc=sc+skp(nck) *q(k) 

400  CONTINUE 

LOADING  FACTOR 

IF(NITSR.EQ.l)  THEN 

compute  factor  of  tendon  force  to  include  

npt=npnf+npff+nptp+npp4+npdp+netp+nesw+nitsr 

nptl=npt-npdp 

fac=facal (nptl) +alfa*facil (npt) 
else 

fac=ZERO 

endif 

pSL=SQRT ( SA*SA+SB*SB+SC*SC) 
if  (sa.lt. zero)  dsl=-dsl 
SEL=fac*SEL-dsl 
STOT=SEL/AREA 
VTOT=STOT/EPROP ( 1 ) 

CALL  FORVER(sel,B,FT,XX,YY,ZZ) 


C 
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C TRANSFORMATION  FROM  LOCAL  TO  GLOBAL 

C 

CALL  TRAN3 (TT,FT,FE, 18) 

C 

C CALCULATE  THE  DAMPING  FORCE 

C 


DO  500  I=l,NDOF 
500  FD(I)=ZERO 
1000  RETURN 
END 


C*********************** ********** ********** 

C SUBROUTINE  TRAN3 

C******************************************* 
SUBROUTINE  TRAN3 (T, FL, F, NDF) 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  T ( 3 , 3 ) , F (NDF) , FL ( 12 ) 
COMMON/PRECIS/ZERO , ONE 
DO  100  1=1, NDF 

100  F(I)=ZERO 

DO  200  1=1,3 

200  F(I)=T(1,I) *FL(1) 

DO  300  1=16,18 
11=1-15 

300  F(I)=T(1,II)*FL(12) 

DO  400  1=4,6 
11=1-3 


**************** 

* 

**************** 


DO  400  K=2,3 
KK=K-1 

400  F(I)=F(I)+T(KK,II) *FL(K) 
DO  500  1=7,9 
11=1-6 

DO  500  K=4,5 
KK=K-3 

500  F(I)=F(I)+T(KK,II) *FL(K) 
DO  600  L=6,9,3 
DO  600  J=L+4,L+6 
DO  600  K=L,L+2 
JJ=J-L-3 
KK=K-L+1 

600  F(J)=F(J)+T(KK, JJ) *FL(K) 

DO  700  1=16,18 
11=1-15 

700  F(I)=T(1,II)*FL(12) 

RETURN 
END 


C**************** ************** ******** 

C SUBROUTINE  FORVER  

C* ************************************* 


******************* 

******************* 


SUBROUTINE  FORVER (T1 , B , FT , XX, YY , ZZ ) 
IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  B(7) ,FT(12) ,XX(6) ,YY(6) ,ZZ(6) 
FT(1)=-T1  ' 


* 

* 

* 


onooooooo  nooooo 


143 


FT(2)=T1*B(6) 

FT(3)=T1*B(1) 

FT(4)=-FT(2) 

FT(5)=-FT(3) 

FT(6)=T1*B(7) 

FT(7)=-T1*B(2) 

FT(8)=-T1*B(5) 

FT{9)=-FT(6) 

FT(10)=-FT{7) 

FT(11)=-FT(8) 

FT(12)=T1 

IF  (XX(1) .EQ.XX(2) ) THEN 
FT(1)=0.0 
FT(2)=-T1*B(3) 

ENDIF 

IF  (XX(6) .EQ.XX(5) ) THEN 
FT(9)=T1*B(4) 
FT(12)=0.0 
ENDIF 
RETURN 
END 


SUBROUTINE  INITll (NJT,NDF,RF, FACIL) 


***************************************************** 

SUBROUTINE  TO  COMPUTE  THE  INITIAL  FORCE 

ELEMENT  TYPE  = 11  THREE  DIMENSIONAL  TENDON  ELEMENTS. 

**********************************1fk*it**ifk******1i1t1t 

IMPLICIT  DOUBLE  PRECISION (A-H, 0-Z) 

COMMON  /INFGR  / NGR,NELS,MFST, IGRHED(IO) ,NINFC,LSTAT, 

* LSTF , LSTC , NDOF , DKO , DKT , EPROP (1000) 

COMMON  /INFEL  / IMEM, KST, LM ( 18) ,NODE (6) , XX ( 6) , YY (6) , 

* ZZ(6) ,PROP(l) ,AREA,XL,B(7) ,SEL,VTOT, 

* STOT,TT(3,3) ,T,SKP(324) 

COMMON  /PRECIS/  ZERO, ONE 

COMMON  /STLDPT/  NPNF,NPFF,NPTP,NPP4 ,NPDP,NETP,NESW,NITSR, 

* NSPATT,NSEG 
DIMENSION  RF(NDF) ,fi (12) , FACIL (1) 

-CHECK  IF  INITIAL  FORCE  IS  TO  BE  INCLUDED  

IF(NITSR.EQ. 1)  THEN 

-compute  factor  of  tendon  force  to  include  

npt=npnf+npf f+nptp+npp4+npdp+netp+nesw+nitsr 

-for  facal,  nptl  = npt  - npdp 
-fac  = facal (nptl)  + alfa*facil (npt) 
fac=facil (npt) 

COMPUTE  INITIAL  FORCE 

P=FAC*T 
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CALL  FORVER(P,B,FI,XX,YY,ZZ) 

C 

C TRANSFORMATION  FROM  LOCAL  TO  GLOBAL 
C 

CALL  TRAN3 (TT,FI,RF,18) 

ENDIF 

RETURN 

END 



SUBROUTINE  OUTS 11  (KPR,TIME) 



C *************************************************** 

C SUBROUTINE  TO  PRINT  TIME  HISTORY  OF  CURRENT  STATE, 

C ELEMENT  TYPE  = 11  THREE  DIMENSIONAL  TENDON  ELEMENTS. 

C *************************************************** 

C 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

COMMON  /TAPES  / NIU, NOU, NTl , NT2 , NT3 , NT4 , NTS , NT6 , NT7 , NTS  , 

* NT9 , NTIO , NTH , NT12  , NT13  , NT14  , NT15 , NT16 
COMMON  /INFGR  / NGR,NELS,MFST, IGRHED(IO) ,NINFC,LSTAT, 

* LSTF , LSTC , NDOF , DKO , DKT , EPROP (1000) 

COMMON  /INFEL  / IMEM, KST, LM ( 18 ) , NODE (6) , XX ( 6) , YY (6) , 

* ZZ (6) ,PROP(l) , AREA, XL, B (7) , SEL, VTOT, STOT, 

* TT(3,3) ,T,SKP(324) 

C 

IF  (IMEM.EQ.MFST)  KHED  = 0 
FTOT=SEL 

IF  (KPR.EQ.O)  RETURN 
C 

IF  (KHED.NE.O)  GO  TO  100 
C 

KHED  = 1 

KKPR  = lABS(KPR) 

WRITE (NOU, 2000)  KKPR, TIME, IGRHED 
C 

100  WRITE (NOU, 2010)  IMEM,NODE ( 1) , NODE (6) , FTOT, STOT, VTOT 

C 

2000  FORMAT  (///18H  RESULTS  FOR  GROUP, 13, 

* 44H  (THREE  DIMENSIONAL  TENDON  ELEMENTS),  TIME  = Ell. 4// 

* 5X,10A4// 

* 5X, 4HELEM, 5X, 4HNODE, 5X, 4HNODE, 8X, 6HTENDON, 8X, 6HTENDON, 

* 5X,9H  TENDON  / 

* 5X, 4H  NO. , 5X, 4H  I,5X,4H  J, 4X, 8X, 6HFORCE  , 8X, 6HSTRESS , 

* 3X,11HDEF0RMATI0N  /) 

2010  FORMAT  (3I9,3E14.4) 

C 

RETURN 

END 


APPENDIX  C 


LISTINGS  OF  THE  TENDON  ELEMENT  WITH  FRICTION 


c 


SUBROUTINE  INELll  (NJT,NDKOD, X, Y, Z , KEXEC) 


C 

C 

C 

C 

C 

C 


C 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 


**********-k******  ******  *******************  ********** 
SUBROUTINE  TO  READ,  GENERATE  AND  PRINT  ELEMENT  DATA. 
ELEMENT  TYPE  =11  THREE  DIMENSIONAL  TENDON  ELEMENTS. 

***********************************************it***i, 

IMPLICIT  DOUBLE  PRECISION (A-H,0-Z) 

DIMENSION  NDKOD(NJT,3) ,X(NJT) ,Y(NJT) ,Z(NJT) ,ICOM(882) 
COMMON  /TAPES  / NIU, NOU, NTl , NT2 ,NT3 ,NT4 ,NT5 , NT6 , NT7 , NT8 , 

* NT9  , NTIO , NTH , NT12  , NT13  , NT14  , NT15 , NT16 
COMMON  /INFGR  / NGR, NELS ,MFST, IGRHED (10) ,NINFC, LSTAT, 

* LSTF , LSTC , NDOF , DKO , DKT , EPROP (5,200) 
COMMON  /INFEL  / IMEM, KST, LM( 18 ) , NODE (6) , INS (6) , INSO (6) , 

* XX(6),  YY(6) ,ZZ(6) ,PROP(5) ,AREA,XL(5) ,B(7) ,SEL(5) , 

* STOT(5) ,TT(3,3) ,T,TI(5) ,QL(12) ,SHA1,SHA2, 

* BK(4,4) ,SK(5) ,IP,IDUM,ANSL,SKP(324) 

COMMON  /PRECIS/  ZERO, ONE 

COMMON  /FRICT  /F (5) , AK(4 , 4) , SF (4 , 4 ) , R (4) , SS (4 , 4 ) , RR(4 , 4) 
COMMON  /WORK  / NODC (6) , XE ( 6) , YE ( 6) , ZE (6) 

EQUIVALENCE  ( IMEM, I COM ( 1 ) ) 

DIMENSION  I AST (3) 

DATA  lAST  /2H  ,2H  *,2H**/ 


MEANINGS  OF  VARIABLES  IN  COMMON  BLOCK  /INFEL/ 


IMEM 

KST 

LM(18) 

NODE (6) 

INS (6) 

INSO (6) 

PROP (5) 

AREA 

XL(5) 

XX(6) 

B(8) 

TT(3,3) 

STOT(5) 

SEL(5) 

SKP(324 


- ELEMENT  NUMBER. 

- STIFFNESS  UPDATE  CODE. 

- LOCATION  MATRIX. 

- NODE  NUMBERS  FROM  END  I TO  J. 

- CODE  OF  NODE'S  SLIDING. 

- PREVIOUS  CODE  OF  NODE'S  SLIDING. 

- MATERIAL  PROPERTIES,EA, fs, fd,ds,dd. 

- CROSS  SECTIONAL  AREA. 

- ELEMENT  LENGTH  OF  EACH  SEGMENT. 

- LOCAL  COORDINATES  OF  NODE  I TO  J. 

- DIRECTION  PARAMETERS. 

- TRANSFORMATION  MATRIX. 

- STRESS  OF  EACH  SEGMENT. 

- AXIAL  FORCE  OF  EACH  SEGMENT. 

- CURRENT  TANGENT  OR  EFFECTIVE  STIFFNESS  MATRIX. 


ELEMENT  GROUP  PARAMETERS 

IF  (MFST.LE.O)  MFST  = 1 

NDOF  =18 
LSTAT  = 234 
LSTF  =648 
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147 


C 

100 

C 

C 

C 

C 

C 

C 


C 


110 


c 

c 

c 


c 

140 

C 


C 

150 

C 

155 


C 

160 

170 

C 


LSTC  = 0 
NINFC  = 882 
KST  = 1 

DO  100  1=3,882 
ICOM(I)  = 0 

PRINT  ELEMENT  GROUP  INFORMATION 

WRITE (NOU ,2000)  NGR , IGRHED , NELS , MFST , NDOF , NINFC , DKO , DKT 

READ  AND  PRINT  MATERIAL  PROPERTIES  & ANCHORAGE  SLIP 
LENGTH 

READ  (NIU,1000)  NMAT 
IF  (NMAT.LE.O)  NMAT  = 1 

WRITE (NOU, 2010)  NMAT 
DO  110  1=1, NMAT 

READ  (NIU,1005)  M, (EPROP( J,M) , J=l, 5) 

WRITE (NOU, 2020)  (M, (EPROP(J,M) ,J=1,5) ,M=1,NMAT) 

READ  (NIU,1008)  NT,IP,ANSL 
WRITE (NOU, 2025)  NT,IP,ANSL 
WRITE (NOU, 2026) 

READ,  GENERATE  AND  PRINT  ELEMENT  DATA 

IMEM  = MFST 

NLAST  = MFST  + NELS  - 1 

READ  (NIU,1010)  N,NODC,MT,AO,FORC,ND 

IF  (MT.LE.O)  MT  = 1 
IF  (ND.EQ.O)  ND  = 1 

IF  (N  - IMEM)  155,160,210 

JFAST  = IAST(3) 

KEROR  = 1 

WRITE (NOU ,2040)  JFAST , N , NODC , MT , AO , FORC 
GO  TO  290 

DO  170  1=1,6 
NODE (I)  = NODC (I) 

KEROR  = 0 
NDIF  = ND 
MTYP  = MT 
AREA  = AO 
T = FORC 
JFAST  = lAST(l) 

IF  (N.NE. NLAST)  GO  TO  140 
GO  TO  210 


C 
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190  DO  200  1=1,6 
200  NODE (I)  = NODE (I)  + NDIF 
JFAST  = IAST(2) 

C 

210  WRITE ( NOU ,2040) JFAST , IMEM , NODE , MTYP , AREA , T 

C 

C ELEMENT  CONNECTIVITY  AND  DIRECTION 

C TRANSFORMATION  OF  COORDINATES 

C 

DO  220  1=1,6 
NOD  = NODE (I) 

II  = 3*1 

CALL  NCODLM(LM(II-2) ,NJT,NOD, 1) 

CALL  NCODLM(LM(II-l) ,NJT,NOD,2) 

CALL  NCODLM(LM(II) ,NJT,NOD,3) 

XE(I)=X(NOD) 

YE(I)=Y(NOD) 

220  ZE(I)=Z(NOD) 

CALL  TRAN1(XE,YE,ZE,XX,YY,ZZ,TT) 

C 

Z4=ZZ (5) -ZZ (4) 

DZ4=DABS(Z4) 

DX1=XX(2) -XX(1) 

IF(XX(1) .EQ.XX(2) ) DXl=ZERO 
DX2=XX(3) -XX(2) 

DX3=XX(4)-XX(3) 

DX4=XX(5) -XX (4) 

DX5=XX(6)-XX(5) 

IF(XX(6) .EQ.XX(5) ) DX5=ZERO 
DY2=YY(2)-YY(3) 

DY4=YY(5) -YY(4) 

XL(1)=DX1 

XL ( 2 ) =SQRT ( DX2  *DX2+DY2  *DY2 ) 

XL(3)=DX3 

XL(4)=SQRT(DZ4*DZ4+DX4*DX4+DY4*DY4) 

XL(5)=DX5 

CY1=DY2/XL(2) 

CY2=DY4/XL(4) 

SY1=DX2/XL(2) 

CX2=DX4/XL(4) 

CZ2=Z4/XL(4) 

SX11=1.-SY1 

CX22=1.-CX2 

C 

C INITIALIZE  DIRECTION  PARAMETERS 

C 

C 

B(1)=CY1 

B(2)=CY2 

B(3)=SY1 

B(4)=CX2 

B(5)=CZ2 


149 


B(6)=SX11 

B(7)=CX22 

PROP ( 1 ) =EPROP ( 1 , MT ) * AREA 
PROP ( 2 ) =EPROP ( 2 , MT ) 

PROP ( 3 ) =EPROP ( 3 , MT) 

PROP { 4 ) =EPROP ( 4 , MT ) 

PROP ( 5 ) =EPROP ( 5 , MT ) 

C 

C INITIALIZE  INS,INSO,IST 
C 

DO  250  1=1,6 
INS(I)=0 
INSO(I)=0 
250  CONTINUE 
C 

C INITIALIZE  QL 

C 

DO  260  1=1,12 
260  QL(I)=ZERO 
C 

C COMPUTE  THE  FRICTION  LOSS  DURING  TENSIONING 

C 

ANG1=DAC0S(B(3) ) 

ANG2=DACOS(B(4) ) 

SHA1=DSIN(ANG1*0.5D0) 

SHA2=DSIN (ANG2  *0 . 5D0) 
write(*,*) 'angl,2 ' ,angl,ang2 
C if  the  tendon  is  not  symmetrical  then  stop 
IF(NT.NE.2)  go  to  265 
DIV=DABS (ANG1-ANG2 ) /DMINl (ANGl , ANG2 ) 
IF(DIV.GT. l.E-5)  THEN 
KEROR=2 
GO  TO  290 
ENDIF 
C 

265  IF(XL(1) .EQ.ZERO)  THEN 
TI(l)=ZERO 
TI(2)=T 
IB=3 
ELSE 

TI(1)=T 

IB=2 

ENDIF 

C 

IF (XL (5) .EQ.ZERO)  THEN 
TI(5)=ZERO 
IE=4 
ELSE 
IE=5 
ENDIF 

C for  tensioning  at  both  ends 
IF(NT.EQ.2)  IE=3 
DO  310  I=IB,IE 
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IF(I.LT.4)  THEN 
COFF=PROP ( 3 ) *SHA1 
ELSE 

COFF=PROP ( 3 ) *SHA2 
ENDIF 

F(I)=2.*TI(I-1) *COFF/(l.+COFF) 
TI(I)=TI(I-1)-F(I) 

310  CONTINUE 

C for  tensioning  at  two  ends 

IF(NT.EQ.l)  go  to  311 
TI(4)=TI(2) 

TI(5)=TI(1) 

C 

C PRINT  THE  INITIAL  FORCE  AFTER  THE  FRICTION  LOSS 

C 

311  WRITE (NOU, 3000) 

WRITE (NOU, 3050) IMEM,NODE(l) ,NODE(6) 

DO  320  1=1,5 

320  WRITE (NOU, 3100)  NODE (I) , NODE (I+l) , TI (I) 

C 

C EVALUATE  THE  LOSS  DUE  TO  ANCHORAGE  SLIP 
C CALCULATE  THE  STIFFNESS  OF  EACH  SEGMENT 
C 

DO  350  1=1,5 
IF(XL(I) .EQ.ZERO)  THEN 
SK(I)=ZERO 
else 

SK(I)=PROP(l)/XL(I) 

ENDIF 

350  CONTINUE 

c 

IF (PROP (2) .EQ.ZERO)  GO  TO  410 
C 

C INITIALIZE  THE  TENSION  LOSS 

C 

IF(XL(1) .NE.ZERO)  THEN 
TLOSS=ANSL*SK(l) 

TI(l)=TI(l)-TLOSS 

ELSE 

TLOSS=ANSL*SK(2) 

TI(2)=TI(2)-TLOSS 

ENDIF 

C 

C 

C FORM  THE  SLIP  STIFFNESS  MATIX  (W/0  FRICTION) 

C 

CALL  BSTIFF(SK,BK) 

C 

C FORM  THE  FRICTION  TERMS 

C ASSUME  THE  TENDON  IS  TENSIONED  AT  NODE  I 

C for  one  tensioning  end 

C 


COFFl=PROP ( 3 ) *SHA1 


non 
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C0FF2=PR0P ( 3 ) *SHA2 
AK(l,l)=coffl*(SK(l)-SK(2) ) 

AK(l,2)=coffl*SK(2) 

AK(2,l)=-coffl*SK(2) 

AK(2,2)=coffl*(SK(2)-SK(3)  ) 

AK(2,3)=coffl*SK(3) 

AK(3,2)=-coff2*SK(3) 

AK(3,3)=coff2*(SK(3)-SK(4) ) 

AK(3,4)=coff2*SK(4) 

AK(4,3)=“COff2*SK(4) 

AK(4,4)=coff2*(SK(4) -SK(5) ) 

C for  tensioning  at  two  ends 
If(nt.eq.l)  go  to  360 
ak(3,2)=-ak(3,2) 
ak(3,3)=-ak(3,3) 
ak(3,4)=-ak(3,4) 
ak(4,3)=-ak(4,3) 
ak(4,4)=-ak(4,4) 

FORM  THE  TOTAL  SLIP  STIFFNESS  MATRIX 

360  CALL  TSTIFF(AK,BK,SF) 
c 

c check  node  II  if  the  loss  of  anchorage  slip  affects 
c ti(2) 
c 

ns=0 

pro=PROP ( 2 ) /PROP ( 3 ) 
if (xl(l) .eq.ZERO)  go  to  500 
call  zeror(r) 

ffl=(ti(l)+ti(2) ) *coffl*pro 
dfl=ti(2)-ti(l) 
if (dfl.le.ffl)  go  to  700 
ns=l 

R(l) =dfl-ff 1/pro 
call  sol(r,sf,ns,2,ss,rr) 
call  force(ti,sk,r,2) 
c 

c check  node  III  if  the  loss  of  anchorage  slip  affects 
c ti(3) 
c 

500  call  zeror(r) 

ffl=(ti(2)+ti(3) ) *cof f l*pro 
dfl=ti(3) -ti(2) 
if (dfl. le. f fl)  go  to  700 
for  tensioning  at  two  ends 
if(nt.eq.2)  then 
if (xl (1) .eg. zero)  then 
ns=2 
idd=4 
else 
ns=4 
idd=5 


C 
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endif 

ti(5)=ti(l) 
ti(4)=ti(2) 
r ( 2 ) =df 1-f f i/pro 
r(3)=-r(2) 

call  sol (r, sf , ns, idd, ss , rr) 
call  force(ti,sk,r,5) 
go  to  710 
endif 
ns=ns+l 

R(2) =dfl-ffl/pro 
call  sol(r,sf ,ns,3,ss,rr) 
call  force(ti,sk,r,3) 
c 

c check  node  JJJ  if  the  loss  of  anchorage  slip  affects 
c ti(4)  y ^ 

c 

call  zeror(r) 

ffl= (ti (3) +ti (4) ) *coffl*pro 

dfl=ti(4)-ti(3) 

if (dfl. le. ffl)  go  to  700 

ns=ns+l 

R(3) =dfl-ffi/pro 
call  sol (r , sf ,ns, 4 , ss, rr) 
call  force(ti,sk,r,4) 
c 

c check  node  JJ  if  the  loss  of  anchorage  slip  affects 
c ti(5) 
c 

if (xl(5) .eq.ZERO)  go  to  700 
call  zeror(r) 

ffl=(ti(4)+ti(5) ) *coffl*pro 

dfl=ti(5)-ti(4) 

if (df 1 . le. f fl)  go  to  500 

ns=ns+l 

R(4) =dfl-ffl/pro 

call  sol (r, sf , ns, 5, ss, rr) 

call  force(ti,sk,r,5) 

C 

C PRINT  THE  INITIAL  FORCE  AFTER  ANCHORAGE  SLIP 
C 

700  if(nt.eq.2)  then 
ti(4)=ti(2) 
ti(5)=ti(l) 
endif 

710  WRITE (NOU, 3200) 

WRITE(NOU, 3050) IMEM,NODE(l) ,NODE(6) 

DO  400  1=1,5 

400  WRITE (NOU, 3100)  NODE (I) , NODE (I+l) , TI ( I) 

C COMPUTE  STIFFNESS  MATRIX  PROFILE  AND  TRANSFER  ELEMENT 
C DATA  TO  TAPE 

C 
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410 

C 

C 

C 


C 

C 

C 

290 

C 


300 

1000 

1005 

1008 

1010 

2000 


2010 

2020 


2025 


2026 

2040 


2050 


CALL  BAND(LM,NDOF,NDKOD) 

CALL  COMPAC 

CHECK  LAST  ELEMENT 

IF  (IMEM.EQ.NLAST)  GO  TO  290 
IMEM  = IMEM  + 1 
IF  (IMEM.EQ.N)  GO  TO  160 
GO  TO  190 


ERROR  AND  ELEMENT  DATA  GENERATION  MESSAGES 


IF  (KEROR.EQ.O)  GO  TO  300 


IF(KEROR.EQ.2)  THEN 
WRITE (NOU, 2061) 

CALL  EXIT 
ENDIF 

WRITE (NOU, 2060) 

CALL  EXIT 
RETURN 
FORMAT  (15) 

FORMAT  (I5,5F10.0) 

FORMAT  (2I5,F10.0) 

FORMAT  (8I5,2F10.0,1X,I5) 

FORMAT  (26H  ELEMENT  GROUP  INDICATOR  = 13, 


* 36H  (THREE  DIMENSIONAL  TENDON  ELEMENTS) //5X, 10A4// 

* 5X,43HNUMBER  OF  ELEMENTS  IN  THIS  GROUP = 15/ 

* 5X,43HNUMBER  OF  FIRST  ELEMENT  IN  THIS  GROUP  ...  = 15// 

* 5X,43HNUMBER  OF  DEGREES  OF  FREEDOM  PER  ELEMENT  = 15/ 

* 5X,43HLENGTH  OF  ELEMENT  INFORMATION  ARRAY  = 15// 

* 5X,43HDAMPING  COEFFICIENT,  BETA-0  = FI 1.5/ 

* 5X,43HDAMPING  COEFFICIENT,  BETA-T  = F11.5) 

FORMAT  (//20H  MATERIAL  PROPERTIES// 

* 5X,36HNUMBER  OF  DIFFERENT  MATERIAL  TYPES  = 15//) 

FORMAT  (5X,38HMATERIAL  NUMBER  = 19/ 

* 5X,43HMODULUS  OF  THE  TENDON  = E12.5/ 

* 5X,43HSTA.  COFF.  OF  FRIC.  BETWEEN  DUCT  & TENDON  = E12!5/ 

* 5X,43HDYN.  COFF.  OF  FRIC.  BETWEEN  DUCT  & TENDON  = E12.5/ 

* 5X,43HSTA.  COFF.  OF  FRIC.  BETWEEN  CON.  & TENDON  = E12.5/ 

* 5X,43HDYN.  COFF.  OF  FRIC.  BETWEEN  CON.  & TENDON  = E12.5) 

FORMAT  (//40H  TENSIONING  METHOD  & ANCHORAGE  PROPERTY  // 

* 5X,43HNUMBER  OF  TENSIONING  ENDS  = 15/ 

* 5X, 43HPRESTRESSING  METHODS  = 15/ 

* 5X,43HANCHORAGE  SLIP  LENGTH  FOR  EACH  TENDON  ...  = F12.5) 

FORMAT  (//20H  ELEMENT  INFORMATION//) 

FORMAT  (/A2/ 

* 5X,43HELEMENT  NUMBER  = 15// 

* 5X,10HNODE  I =18, 5X, lOHNODE  II  =18 , 5X, lOHNODE  III  =18/ 

* 5X, lOHNODE  JJJ  =18, 5X, lOHNODE  JJ  =18 , 5X, lOHNODE  J =18/ 

* 5X, 13HMATERIAL  NO  =I7,5X,8H  AREA  =F12.4/ 

* 5X,15HINITIAL  FORCE  =F12.4/) 

FORMAT  (1X,A2,I6,7I9,2F12.4) 
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2060  FORMAT  (///15H  ERROR  MESSAGES// 

* 5X,  3HKEY,  5X,  29H MESSAGE , 5X,  5HERROR// 

* 5X,3H  ,5X,29HELEMENT  DATA  INPUT  DIRECTLY  ,5X,5HNONE  / 

* 5X;3H  *,5X,29HELEMENT  DATA  GENERATED  ,5X,5HNONE  / 

* 5X,3H  **,5X,29HELEMENT  CARD  OUT  OF  SEQUENCE  , 5X , 5HFATAL/ ) 

2061  FORMAT  (///15H  ERROR  MESSAGES// 

* 5X,43H TENDON  IS  NOT  SYMMETRICAL /) 

3000  FORMAT(//33HINITIAL  FORCE  AFTER  FRICTION  LOSS//) 

3050  FORMAT (5X,8HELEM  NO. , 15, 5X, 5HNODE1, 15, 5X, 5HNODE6, 15// 

* 5X,4HN0DE,1X,3H— ,1X,4HNODE,5X,10HTEN.  FORCE/) 

3100  FORMAT (2 19, E15. 4) 

3200  FORMAT (// 3 4HINITIAL  FORCE  AFTER  ANCHORAGE  SLIP//) 

END 

C* ********************************************************* 

C SUBROUTINE  TRANl-TRANSFORMATION  of  COORDINATES  * 

C*  ********************************************************* 

SUBROUTINE  TRANl (X, Y, Z , XX, YY, ZZ , T) 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  X(6),Y(6),Z(6),T(3,3),XX(6),YY(6) 

* ,ZZ(6) ,V(4,4) 

CALL  VECTOR (V( 1,1) ,X(3),Y(3),Z(3),X(4),Y(4), 

* Z (4)  ) 

CALL  VECTOR (V(l, 4) ,X(3),Y(3),Z(3),X(1),Y(1),Z(1)) 

CALL  CROS(V(l,l) ,V(1,4) ,V(1,3)) 

CALL  CROS(V(l,3) ,V(1,1) ,V(1,2) ) 

DO  100  K=l,3 
DO  100  L=l,3 
100  T(K,L)=V(L,K) 

DO  200  N=l,6 

XX(N)=T(1, 1) *X(N)+T(1,2) *Y(N)+T(1,3) *Z(N) 

YY(N)=T(2,1) *X(N)+T(2,2)*Y(N)+T(2,3) *Z(N) 

ZZ(N)=T(3,1) *X(N)+T(3,2) *Y(N)+T(3,3) *Z (N) 

200  CONTINUE 
RETURN 
END 

: VECTOR 

SUBROUTINE  VECTOR  (V, XI , YI , ZI , XJ , YJ , Z J) 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  V(4) 

X=XJ-XI 

Y=YJ-YI 

Z=ZJ-ZI 

V(4)=SQRT(X*X+Y*Y+Z*Z) 

V(3)=Z/V(4) 

V(2)=Y/V(4) 

V(1)=X/V(4) 

RETURN 

END 

' CROSS 

SUBROUTINE  CROS(A,B,C) 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  A(4) ,B(4) ,C(4) 

COMMON  /lOLIST/  NTM,NTR,NIN,NOT,NTl,NFL,NT2 ,NT3 ,NT4 ,NT5 


o o o 
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X=A(2) *B(3)-A(3) *B(2) 
Y=A(3) *B(1)-A(1) *B(3) 
Z=A(1)*B(2)-A(2)*B(1) 

XYZ  = X*X  + Y*Y  + Z*Z 
IF(XYZ.GT.O.O)  GO  TO  100 
WRITE  (NTM,2000) 

WRITE  (NOT, 2000) 

STOP 

C(4)=SQRT(XYZ) 

C(3)=Z/C(4) 

C(2)=Y/C(4) 

C(1)=X/C(4) 

RETURN 


2000 


FORMAT 

END 


('  **  K NODE  IS  ALONG  ELEMENT  SEE  MANUAL  **'/) 


SUBROUTINE  BSTIFF ( SK , BK) 

IMPLICIT  DOUBLE  PRECISION (A-H, 0-Z ) 
COMMON  /PRECIS/  ZERO, ONE 
DIMENSION  SK(5) ,BK(4,4) 

DO  150  J=l,4 
DO  150  1=1,4 

150  BK(I,J)=ZERO 
DO  170  1=1,4 

170  BK(I,I)=SK(I)+SK(I+1) 

BK(1,2)=-SK(2) 

BK(2,1)=BK(1,2) 

BK(2,3)=-SK(3) 

BK(3,2)=BK(2,3) 

BK(3,4)=-SK(4) 

BK(4,3)=BK(3,4) 

RETURN 

END 



SUBROUTINE  TSTIFF (AK, BK, SF) 

IMPLICIT  DOUBLE  PRECISION (A-H , 0-Z ) 
DIMENSION  AK(4,4) ,BK(4,4) ,SF(4,4) 

TOTAL  STIFFNESS 


•BSTIFF 


■TSTIFF 


5000 


DO  5000  J=l,4 

DO  5000  1=1,4 

SF(I,J)=BK(I,J)+AK(I,J) 

RETURN 

END 


SUBROUTINE  ZEROR(R) 

IMPLICIT  DOUBLE  PRECISION (A-H, 0-Z ) 
COMMON  /PRECIS/  ZERO, ONE 
dimension  r(4) 
do  100  i=l,4 


■ZEROR 


156 


100  r(i)=ZERO 
return 
end 

C FORCE 

C 

C CALCULATE  THE  TENDON  FORCES 

C 

SUBROUTINE  FORCE (sel , sk,R, ind) 

IMPLICIT  DOUBLE  PRECISION (A-H,0-Z) 
dimension  sel ( 5 ) , sk ( 5 ) , r ( 4 ) 
in=ind 

SEL(1)=SEL(1)+SK(1)  *R(1) 
if(ind.eq.5)  in=4 
DO  300  1=2, in 
11=1-1 

300  SEL(I)=SEL(I)+SK(I) *(R(I)-R(I1) ) 
if(ind.le.4)  return 
SEL(5)=SEL(5) -SK(5) *R(4) 

RETURN 

END 



SUBROUTINE  STIFll ( ISTEP , NDF , CDKO , CDKT , FK , INDFK , ISTFC) 



C ******************************************************* 

C SUBROUTINE  TO  COMPUTE  ELEMENT  TANGENT  STIFFNESS  MATRIX. 

C ELEMENT  TYPE  =11  THREE  DIMENSIONAL  TENDON  ELEMENTS. 

C ******************************************************* 

C 

IMPLICIT  DOUBLE  PRECISION (A-H, 0-Z ) 

DIMENSION  FK (NDF, NDF) 

COMMON  /INFGR  / NGR,NELS ,MFST, IGRHED (10) , NINFC, LSTAT, 

* LSTF,LSTC,NDOF,DKO,DKT,EPROP(5,200) 
COMMON  /INFEL  / IMEM, KST, LM(18) , NODE ( 6) , INS (6) , INSO ( 6) , 

* XX(6) ,YY(6) ,ZZ(6) , PROP (5) ,AREA,XL(5) ,B(7) ,SEL(5) , 

* STOT(5) ,TT(3,3) ,T,TI(5) , QL(12 ) , SHAl , SHA2 , 

* BK(4,4) ,SK(5) ,IP,IDUM,ANSL,SKP(324) 

COMMON  /STLDPT/  NPNF, NPFF, NPTP, NPP4 , NPDP,NETP, NESW, NITSR, 

* NSPATT  NSEG 

COMMON  /WORK  / ST ( 12 , 12 ) , SSK ( 18 , 18 ) , TS ( 12 , 12 ) 

COMMON  /PRECIS/  ZERO, ONE 

C set  to  full  matrix  version 

INDFK=0 

c 

c LINEAR  PART  OF  EFFECTIVE  STIFFNESS 

C 

C COMPUTE  STIFFNESS 

C 

DO  50  1=1,12 
DO  50  J=l,12 
50  ST(I,J)=ZERO 

C 
C 
C 


TURN  OFF  STIFFNESS  AT  TRANSFER  FOR  POST-TENSIONED 


CJCJUU  uuuu 


c 
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IF(IP.EQ.O.AND.NSEG.EQ.l)  GO  TO  510 
IF(NSEG.LT. IP)  GO  TO  510 

CHECK  NODE (II) 

IF(INS(2) .EQ.O)  THEN 
TRUSS (I, II) 

IF(XX(1) .EQ.XX(2) ) GO  TO  100 
TPR=PROP(l)/XL(l) 

ST(1,1)=TPR 
ST(1,2)=-ST(1,1) 

ST(2,2)=ST(1,1) 

GO  TO  100 
ENDIF 

INS(2)=1  CHECK  NODE (III) 

GO  TO  150 

CHECK  NODE (III) 

100  IF(INS(3) .EQ.O)  THEN 
C TRUSS (II, III) 

TPR=PROP(l)/XL(2) 
ST(2,2)=ST(2,2)+TPR*B(3) *B(3) 
ST(2,3)=-TPR*B(1) *B(3) 
ST(2,4)=-TPR*B(3) *B(3) 
ST(2,5)=-ST(2,3) 

ST(3,3)=TPR*B(1)  *B(1) 

ST(3,4)=ST(2,5) 

ST(3,5)=-ST(3,3) 

ST(4,4)=-ST(2,4) 

ST(4,5)=ST(2,3) 

ST(5,5)=ST(3,3) 

GO  TO  200 
ENDIF 

GO  TO  210 

150  IF(INS (3) .EQ. 0)  THEN 
C TENDON (I-III) 

TPR=PROP(l)/(XL(l)+XL(2) ) 

ST(1, 1)=TPR 
ST(1,2)=-TPR*B(6) 

ST(1,3)=-TPR*B(1) 

ST(1,4)=-TPR*B(3) 

ST(1,5)=-ST(1,3) 

ST(2,2)=TPR*B(6) *B(6) 

ST(2,3)=TPR*B(6) *B(1) 

ST(2,4)=TPR*B(6) *B(3) 

ST(2,5)=-ST(2,3) 

ST(3,3)=TPR*B(1) *B(1) 

ST(3,4)=TPR*B(1) *B(3) 

ST(3,5)=-ST(3,3) 

ST(4,4)=TPR*B(3) *B(3) 

ST(4,5)=-ST(3,4) 

ST(5,5)=-ST(3,5) 
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GO  TO  200 
ENDIF 

GO  TO  230 

C CHECK  NODE(JJJ) 

200  IF(INS(4) .EQ.O)  THEN 

C TRUSS  (III-JJJ) 

TPR=PROP(l)/XL(3) 
ST(4,4)=ST(4,4)+TPR 
ST(4,6)=-TPR 
ST(6,6)=TPR 
GO  TO  300 
ENDIF 

GO  TO  310 

210  IF(INS(4) .EQ.O)  THEN 

C TENDON  ( II- JJJ) 

TPR=PROP(l)/(XL(2)+XL(3) ) 
ST(2,2)=ST(2,2)+TPR*B(3) *B(3) 
ST(2,3)=-TPR*B(1) *B{3) 

ST ( 2 , 4 ) =TPR*B ( 3 ) *B ( 6 ) 

ST(2,5)=-ST(2,3) 

ST(2,6)=-TPR*B(3) 

ST(3,3)=TPR*B(1) *B(1) 

ST(3,4)=-TPR*B(1) *B(6) 

ST(3,5)=-ST(3,3) 

ST(3,6)=TPR*B(1) 

ST(4,4)=TPR*B(6) *B(6) 

ST(4,5)=-ST(3,4) 

ST(4,6)=-TPR*B(6) 

ST(5,5)=ST(3,3) 

ST(5,6)=-ST(3,6) 

ST(6,6)=TPR 
GO  TO  300 
ENDIF 

GO  TO  330 

230  IF(INS(4) .EQ.O)  THEN 

C TENDON  (I-JJJ) 

TPR=PROP(l)/(XL(l)+XL(2)+XL(3) ) 

ST(1,1)=TPR 

ST(1,2)=-TPR*B(6) 

ST(1,3)=-TPR*B(1) 

ST(1,4)=-ST(1,2) 

ST(1,5)=-ST(1,3) 

ST(1,6)=-ST(1,1) 

ST ( 2 , 2 ) =TPR*B ( 6 ) *B ( 6 ) 
ST(2,3)=TPR*B(6) *B(1) 
ST(2,4)=-ST(2,2) 
ST(2,5)=-ST(2,3) 

ST(2,6)=-ST(1,2) 

ST(3,3)=TPR*B(1) *B(1) 
ST(3,4)=-ST(2,3) 

ST(3,5)=-ST(3,3) 

ST(3,6)=-ST(1,3) 

ST(4,4)=-ST(2,4) 
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ST(4,5)=-ST(2,5) 

ST(4,6)=-ST(2,6) 

ST(5,5)=-ST(3,5) 

ST(5,6)=-ST(3,6) 

ST(6,6)=-ST(1,6) 

GO  TO  300 
ENDIF 

GO  TO  370 

C CHECK  NODE(JJ) 

300  IF(INS(5) .EQ.O)  THEN 

C TRUSS (JJJ,JJ)  TRUSS (JJ,J) 

CALL  TRU46(ST,PROP(l) ,XL,B,XX) 

GO  TO  400 
ENDIF 

C TENDON  (JJJ,J) 

CALL  TEN46(ST,PROP(l) ,XL,B) 

GO  TO  400 

310  IF(INS(5) .EQ.O)  THEN 

C TENDON (III, JJ)  + TRUSS (JJ,J) 

CALL  TT36(ST,PROP(l) ,XL,B,XX) 

GO  TO  400 
ENDIF 

C TENDON  ( III,  J) 

CALL  TEN36(ST,PROP(l) ,XL,B) 

GO  TO  400 

330  IF(INS(5) .EQ.O)  THEN 

C TENDON(II, JJ)  + TRUSS (JJ,J) 

TPR=PROP(l)/ (XL(2)+XL(3)+XL(4) ) 

CALL  TEN16(INS,TS,TPR,B) 

DO  335  1=1,12 
DO  335  J=I,12 

335  ST(I, J)=ST(I, J)+TS(I, J) 

TPR=PROP ( 1 ) /XL ( 5 ) 

IF(XX(5) .EQ.XX(6) ) GO  TO  400 
ST(9,9)=ST(9,9)+TPR 
ST(9,12)=-TPR 
ST(12,12)=TPR 
GO  TO  400 
ENDIF 

C TENDON  ( II- J) 

TPR=PROP(l)/(XL(2)+XL(3)+XL(4)+XL(5) ) 
CALL  TEN16(INS,TS,TPR,B) 

DO  338  1=1,12 
DO  338  J=I,12 

338  ST(I,J)=ST(I,J)+TS(I,J) 

GO  TO  400 

370  IF(INS(5) .EQ.O)  THEN 

C TENDON  ( I,  JJ) 

TPR=PROP(l)/ (XL(1)+XL(2)+XL(3)+XL(4) ) 
CALL  TEN16(INS,ST,TPR,B) 

TPR=PROP ( 1 ) /XL ( 5 ) 

IF(XX(5) .EQ.XX(6) ) GO  TO  400 
ST(9,9)=ST(9,9)+TPR 


o o o 
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ST(9,12)=-TPR 
ST(12,12)=TPR 
GO  TO  400 
ENDIF 

C TENDON  ( I,  J) 

TPR=PROP(l)/(XL(l)+XL(2)+XL(3)+XL(4)+XL(5) ) 

CALL  TEN16(INS,ST,TPR,B) 

400  DO  500  1=1,12 
DO  500  J=I,12 
ST(J,I)=ST(I, J) 

500  CONTINUE 

C TRANSFORM  THE  LOCAL  STIFFNESS  MATRIX  INTO  GLOBAL  MATRIX 

510  CALL  TRAN2 (TT,ST,SSK, 18, 1) 

IF(ISTFC.EQ.O)  GO  TO  550 
IJ  =0 

DO  140  I=l,NDOF 
c DO  140  J=I,NDOF 

do  140  j=l,ndof 
IJ  = IJ  + 1 

STIF  = SSK(I,J) 

FK  (I,J)  = STIF 

c FK  (J,I)  = STIF 

140  SKP(IJ)  = STIF 

RETURN 

CHANGE  STIFFNESS 

550  IJ=0 

DO  600  I=l,NDOF 
DO  600  J=l,NDOF 
IJ=IJ+1 
STIF=SSK(I, J) 

STIFD=STIF-SKP (IJ) 

FK(I, J)=STIFD 
SKP(IJ)=STIF 
600  CONTINUE 

RETURN 
END 


SUBROUTINE  TRU46 (S, PRO,XL, B,XX) 
IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  3(12,12) ,XL(5) ,B(7) ,XX(6) 
COMMON/PRECIS/ZERO , ONE 

C TRUSS  (JJJ-J) 

TMP= PRO/XL (4) 

S(6,6)=S(6,6)+TMP*B(4)*B(4) 
S(6,7)=TMP*B(4) *B(2) 
S(6,8)=TMP*B(4)*B(5) 
S(6,9)=-TMP*B(4) *B(4) 
S(6,10)=-S(6,7) 

S(6,ll)=-S(6,8) 

S(7,7)=TMP*B(2) *B(2) 
S(7,8)=TMP*B(2)*B(5) 
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S(7,9)=-S(6,7) 

S(7,10)=-S(7,7) 

S(7,ll)=-S(7,8) 

S(8,8)=TMP*B(5) *B(5) 

S(8,9)=S(6,11) 

S(8,10)=-S(7,8) 

S(8,ll)=-S(8,8) 

TMP1=PR0/XL(5) 

IF(XX(5) .EQ.XX(6) ) TMP1=ZER0 
S(9,9)=TMP*B(4) *B(4)+TMP1 
S(9,10)=-S(7,9) 

S(9,ll)=-S(6,ll) 

S(9,12)=-TMP1 

S(10,10)=S(7,7) 

S(10,11)=S(7,8) 

S(11,11)=S(8,8) 

S(12,12)=TMP1 

RETURN 

END 

TEN46 

SUBROUTINE  TEN46 (S , PRO, XL, B) 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  S(12,12) ,XL(5) ,B(7) 

C TENDON  (JJJ,J) 

TMP=PRO/(XL(4)+XL(5)  ) 

S(6,6)=S(6,6)+TMP*B(4) *B(4) 

S(6,7)=TMP*B(4) *B(2) 

S(6,8)=TMP*B(4) *B(5) 

S(6,9)=TMP*B(4) *B(7) 

S(6,10)=-S(6,7) 

S(6,ll)=-S(6,8) 

S(6,12)=-TMP*B(4) 

S(7,7)=TMP*B(2) *B(2) 

S(7,8)=TMP*B(2) *B(5) 

S(7,9)=TMP*B(2) *B(7) 

S(7,10)=-S(7,7) 

S(7,ll)=-S(7,8) 

S(7,12)=-TMP*B(2) 

S(8,8)=TMP*B(5) *B(5) 

S(8,9)=TMP*B(5) *B(7) 

S(8,10)=-S(7,8) 

S(8,ll)=-S(8,8) 

S(8,12)=-TMP*B(5) 

S(9,9)=TMP*B(7) *B(7) 

S(9,10)=-S(7,9) 

S(9,ll)=-S(8,9) 

S(9,12)=-TMP*B(7) 

S(10,10)=S(7,7) 

S(10,11)=S(7,8) 

S(10,12)=-S(7,12) 

S(ll,ll)=-S(8,ll) 

S(ll,12)=-S(8,12) 

S(12,12)=TMP 


c 


RETURN 

END 
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C 


C 


C 


SUBROUTINE  TT36 (S , PRO , XL, B, XX) 
IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  S(12,12) ,XL(5) ,B(7) ,XX(6) 
COMMON/PRECIS/ZERO, ONE 
-TENDON ( III, JJ)  + TRUSS (JJ,J) 
TMP=PRO/ (XL ( 3 ) +XL ( 4 ) ) 
S(4,4)=S(4,4)+TMP 
S(4,6)=-TMP*B(7) 

S(4,7)=TMP*B(2) 

S(4,8)=TMP*B(5) 

S(4,9)=-TMP*B(4) 

S(4,10)=-S(4,7) 

S(4,ll)=-S(4,8) 

S(6,6)=TMP*B(7) *B(7) 
S(6,7)=-TMP*B(7) *B(2) 
S(6,8)=-TMP*B(7) *B(5) 
S(6,9)=TMP*B(7) *B(4) 

S(6,10)=-S(6,7) 

S(6,ll)=-S(6,8) 

S(7,7)=TMP*B(2) *B(2) 

S(7,8)=TMP*B(2) *B(5) 
S(7,9)=-TMP*B(2) *B(4) 
S(7,10)=-S(7,7) 

S(7,ll)=-S(7,8) 

S(8,8)=TMP*B(5) *B(5) 
S(8,9)=-TMP*B(5)*B(4) 

S(8,10)=-S(7,8) 

S(8,ll)=-S(8,8) 

TMPl=PRO/XL(5) 

IF(XX(5) .EQ.XX(6) ) TMPl=ZERO 
S ( 9 , 9 ) =TMP*B ( 4 ) *B ( 4 ) +TMP1 
S(9,10)=-S(7,9) 

S(9,ll)=-S{8,9) 

S(9,12)=-TMP1 

S(10,10)=-S(7,10) 

S(10,ll)=-S(7,ll) 

S(ll,ll)=-S(8,ll) 

S(12,12)=-S(9,12) 

RETURN 

END 


SUBROUTINE  TEN36 (S, PRO,XL, B) 
IMPLICIT  REAL*8  (A-H,0-Z) 
DIMENSION  S(12,12) ,XL(5) ,B(7) 
■TENDON  ( III,  J) 


TMP=PRO/ ( XL ( 3 ) +XL ( 4 ) +XL ( 5 ) ) 
S(4,4)=S(4,4) +TMP 
S(4,6)=-TMP*B(7) 
S(4,7)=TMP*B(2) 
S(4,8)=TMP*B(5) 


TT3  6 


TEN3  6 


o o 
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S(4,9)=-S(4,6) 

S(4,10)=-S(4,7) 

S(4,ll)=-S(4,8) 

S(4,12)=-TMP 

S(6,6)=TMP*B(7) *B(7) 

S(6,7)=-TMP*B(7) *B(2) 

S(6,8)=-TMP*B(7) *B(5) 

S(6,9)=-S(6,6) 

S(6,10)=-S(6,7) 

S(6,ll)=-S(6,8) 

S(6,12)=-S(4,6) 

S(7,7)=TMP*B(2)*B(2) 

S(7,8)=TMP*B(2) *B(5) 

S(7,9)=-S(6,7) 

S(7,10)=-S(7,7) 

S{7,ll)=-S(7,8) 

S(7,12)=-S(4,7) 

S(8,8)=TMP*B(5) *B(5) 

S(8,9)=-S(6,8) 

S(8,10)=-S(7,8) 

S(8,ll)=-S(8,8) 

S(8,12)=-S(4,8) 

S(9,9)=-S(6,9) 

S(9,10)=-S(6,10) 

S(9,ll)=-S(6,ll) 

S(9,12)=-S(6,12) 

S(10,10)=-S(7,10) 

S(10,ll)=-S(7,ll) 

S(10,12)=-S(7,12) 

S(ll,ll)=-S(8,ll) 

S(ll,12)=-S(8,12) 

S(12,12)=-S(4,12) 

RETURN 

END 


SUBROUTINE  TEN16 (ID, S , PRO , B) 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  ID(6) ,S(12,12) ,B(7) 

COMMON/PRECIS/ZERO , ONE 

-TENDON ( I, J)  OR  TENDON ( I, JJ)  OR  TENDON (I I, JJ)  OR 
TENDON ( II, J)  ' 

DO  10  1=1,12 
DO  10  J=I,12 
S(I, J)=ZERO 
S(l,l)=ONE 
S(l,2)=-B(6) 

S(l,3)=-B(l) 

S(l,4)=-S(l,2) 

S(l,5)=-S(l,3) 

S(l,6)=-B(7) 

S(1,7)=B(2) 

S(1,8)=B(5) 

S(l,9)=-S(l,6) 


TEN16 


10 
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S(l,10)=-S(l,7) 

S(l,ll)=-S(l,8) 

S(l,12)=-S(l,l) 

S(2,2)=B(6) *B(6) 

S(2,3)=B(6)*B(1) 

S(2,4)=-S(2,2) 

S(2,5)=-S(2,3) 

S(2,6)=B(6)  *B(7) 

S(2,7)=-B(6) *B(2) 

S(2,8)=-B(6) *B(5) 

S(2,9)=-S(2,6) 

S(2,10)=-S(2,7) 

S(2,ll)=-S(2,8) 

S(2,12)=B(6) 

S(3,3)=B(1) *B(1) 

S(3,4)=-B(1)*B(6) 

S(3,5)=-S(3,3) 

S(3,6)=B(1) *B(7) 

S(3,7)=-B(l) *B(2) 

S(3,8)=-B(l) *B(5) 

S(3,9)=-S(3,6) 

S(3,10)=-S(3,7) 

S(3,ll)=-S(3,8) 

S(3,12)=B(1) 

DO  100  1=4,12 
100  S(4,I)=-S(2,I) 

DO  200  1=5,12 
200  S(5,I)=-S(3,I) 

S(6,6)=B(7)*B(7) 

S(6,7)=-B(7) *B(2) 

S(6,8)=-B(7)*B(5) 

S(6,9)=-S(6,6) 

S(6,10)=-S(6,7) 

S(6,ll)=-S(6,8) 

S(6,12)=B(7) 

S(7,7)=B(2)*B(2) 

S(7,8)=B(2) *B(5) 

S(7,9)=-S(6,7) 

S(7,10)=-S(7,7) 

S(7,ll)=-S(7,8) 

S(7,12)=-B(2) 

S(8,8)=B(5) *B(5) 

S(8,9)=-S(6,8) 

S(8,10)=-S(7,8) 

S(8,ll)=-S(8,8) 

S(8,12)=-B(5) 

DO  300  1=9,12 
300  S(9,I)=-S(6,I) 

DO  400  1=10,12 
400  S(10,I)=-S(7,I) 

S(ll,ll)=-S(8,ll) 

S(ll,12)=-S(8,12) 

S(12,12)=ONE 
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C FOR  TENDON ( II- JJ,J) 

IF(ID(2) .EQ.O)  THEN 
DO  500  1=1,12 
500  S(l,I)=ZERO 

S(2,2)=B(3) *B(3) 

S(2,3)=-B(3)*B(1) 

S(2,4)=B(3) *B(6) 

S(2,5)=-S(2,3) 

S(2,6)=-B(3) *B(7) 

S(2,7)=B(3)  *B(2) 

S(2,8)=B(3) *B(5) 

S(2,9)=-S(2,6) 

S(2,10)=-S(2,7) 

S(2,ll)=-S(2,8) 

S(2,12)=-B(3) 

ENDIF 

c for  tendon ( I, II-JJ) 

IF(ID(5) .EQ. 0)  THEN 
DO  600  1=1,12 
600  S(I,12)=ZERO 
S(l,9)=-B(4) 

S(2,9)=B(4)*B(6) 

S(3,9)=B(4)*B(1) 

S(4,9)=-S(2,9) 

S(5,9)=-S(3,9) 

S(6,9)=B(4)*B(7) 

S(7,9)=-B(4)*B(2) 

S(8,9)=-B(4) *B(5) 

S(9,9)=B(4) *B(4) 

S(9,10)=-S(7,9) 

S(9,ll)=-S(8,9) 

ENDIF 

C FOR  TENDON (II-JJ) 

IF ( ID ( 2 ) .EQ.O. AND . ID ( 5 ) .EQ.O)  THEN 

S(l,9)=ZERO 

S(2,9)=-B(4) *B(3) 

ENDIF 

DO  700  1=1,12 
DO  700  J=I,12 
700  S(I,J)=S(I, J) *PRO 

RETURN 
END 



SUBROUTINE  TRAN2 (T, SL, S , ND, NL) 
IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  T(3,3) ,TG(12,18) ,SL(12,12) 
COMMON/PRECIS/  ZERO, ONE 
DO  50  1=1,12 
DO  50  J=l,18 
ST(I, J)=ZERO 
50  TG(I,J)=ZERO 
DO  60  1=1,18 
DO  60  J=l,18 


TRAN2 

,S(18,18) ,ST(12,18) 
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60  S(I,J)=ZERO 

C FORM  THE  TRANSFORMATION  MATRIX  TG 

DO  100  J=l,3 
100  TG(1, J)=T(1, J) 

DO  150  1=2,3 
DO  150  J=4,6 
11=1-1 
JJ=J-3 

150  TG(I,J)=T(II,JJ) 

DO  200  1=4,5 
DO  200  J=7,9 
11=1-3 
JJ=J-6 

200  TG(I, J)=T(II,  JJ) 

DO  250  L=6,9,3 
DO  250  I=L,L+2 
DO  250  J=L+4,L+6 
II=I-L+1 
JJ=J-L-3 

250  TG(I, J)=T(II,JJ) 

DO  300  J=16,18 
JJ=J-15 

300  TG(12, J)=T(1, JJ) 

C EVALUATION  ST=SL*TG 

DO  400  1=1,12 
DO  400  J=l,18 
DO  400  K=l,12 

400  ST(I, J)=ST(I, J)+SL(I,K) *TG(K, J) 

C S=TTG*SL*TG 

DO  500  1=1,18 
DO  500  J=I,18 
DO  500  K=l,12 

500  S(I, J)=S(I, J)+TG(K,I) *ST(K, J) 

DO  600  1=1,18 
DO  600  J=I,18 
600  S(J,I)=S(I, J) 

RETURN 

END 



SUBROUTINE  STATU  (NDF, Q, TIME, FACAL, FACIL, ALFA) 



C **************************************************** 

C SUBROUTINE  FOR  STATE  DETERMINATION  CALCULATIONS. 

C ELEMENT  TYPE  = 11.  THREE  DIMENSIONAL  TENDON  ELEMENTS. 

C ***************************************************** 

C 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

DIMENSION  Q(NDF) , FACAL (1) , FACIL (1) ,IN(5) 

COMMON  /INFGR  / NGR, NELS , MFST, IGRHED ( 10) , NINFC, LSTAT, 

* LSTF , LSTC , NDOF , DKO , DKT , EPROP (5,200) 
COMMON  /INFEL  / IMEM, KST , LM ( 18 ) , NODE ( 6) , INS ( 6) , INSO ( 6 ) , 

* XX(6) ,YY(6) ,ZZ(6) ,PROP(5) , AREA,XL(5) , B(7) ,SEL(5) , 

* STOT(5) ,TT(3,3) ,T,TI(5) ,QL(12) ,SHA1,SHA2, 
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C 

C 

C 

C 

10 

C 

C 

C 

C 


C 

C 

C 


* BK(4,4) ,SK(5) ,IP,IDUM,ANSL,SKP(324) 

COMMON  /PRECIS/  ZERO, ONE 

COMMON  /PRITC  / TOL, WTFAC , RNORM, DENOM, RNUM, WFAC, CONV, 

* GRACE,NIT,MAXIT,ITKOD,ITPR,RNMN 

COMMON  /STLDPT/  NPNF, NPFF, NPTP, NPP4 , NPDP, NETP, NESW, NITSR, 

* NSPATT,NSEG 
COMMON  /TRUSS/  TD ( 5 ) , TRUSF ( 5 ) 

COMMON  /DIST/  R (4 ) , SF (4 , 4 ) , SS (4 , 4 ) , RR(4 ) , AK(4 , 4 ) 


CALCULATE  THE  INTERNAL  RESISTING  FORCE 


SET  INSO(I)=INS 
DO  10  1=1,6 
INSO(I)=INS(I) 

TRANSFORM  Q TO  QL 
CALL  TRAN4(TT,Q,QL,NDOF) 

CALCULATE  THE  AXIAL  FORCE  AND  DEFORMATION  OF  EACH 
SEGMENT (TRUSS) 

TD(1)=QL(2)-QL(1) 

TD(2)  = (QL(4)-QL(2)  ) *B(3)  + (QL(3) -QL(5) ) *B(1) 
TD(3)=QL(6)-QL(4) 

TD(4)=(QL(9) -QL(6) ) *B ( 4 ) + (QL ( 10) -QL(7) ) *B (2 ) + (QL ( 11) - 
* QL(8))*B(5) 

TD(5)=QL(12)-QL(9) 

LOADING  FACTOR 


IF (NITSR. EQ.l)  THEN 

c compute  factor  of  tendon  force  to  include  — 

npt=npnf+npf f+nptp+npp4+npdp+netp+nesw+nitsr 
nptl=npt-npdp 

fac=facal (nptl) +alfa*facil (npt) 
else 

fac=ZERO 

endif 

C 

C INITIALIZE  SEL 

C 


DO  50  1=1,5 

IF(XL(I) .EQ.ZERO)  THEN 
TRUSF(I)=ZERO 
SEL(I)=ZERO 
else 

TRUSF(I)=fac*TI(I)+SK(I) *TD(I) 
SEL(I)=TRUSF(I) 

ENDIF 

50  CONTINUE 
IT=1 


INITIALIZE  INS 


C 

C 

C 


DO  60  1=1,6 
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60  INS(I)=0 

C 

C CALCULATE  THE  DRIVING  FORCE  & DETERMINE  THE  NODE'S  SLIDING 

C 

70  DO  65  1=2,5 

IN(I)=INS(I) 

65  CONTINUE 

NID=0 
NS=0 

DO  100  1=1,4 
11=1+1 

IF(XL(I) .EQ.ZERO.OR.XL(Il) .EQ. ZERO. OR. INS (II) .EQ.l)  THEN 
R(I)=ZERO 
ELSE 

IF(I.GE.3)  THEN 
SHA=SHA2 
ELSE 

SHA=SHA1 

ENDIF 

DF=SEL(I1)-SEL(I) 

FF=(SEL(I)+SEL(I1) ) *SHA*PROP(4) 

IF(PROP(4) .EQ.ZERO)  THEN 
FF1=ZER0 
ELSE 

FF1=FF*PR0P ( 5 ) /PROP ( 4 ) 

ENDIF 

R(I)=DABS (DF) -FF 
IF(R(I) .GT.ZERO)  THEN 
INS (II) =1 
R(I)=R(I)-FF1+FF 
NID=1 

IF(DF.LT.ZERO)  R(I)=-R(I) 

ELSE 

R(I)=ZERO 

ENDIF 

ENDIF 

100  CONTINUE 

IF(NID.EQ.O)  GO  TO  350 
C 

C FOMULATION  OF  THE  STIFFNESS  MATRIX  OF  FORCE  DISTRIBUTION 
C 

IF(IT.GT.l)  GO  TO  180 
CALL  BSTIFF(SK,BK) 

DO  150  J=l,4 
DO  150  1=1,4 
150  AK(I,J)=ZERO 
C 

C FRICTION  TERMS 

C 

180  IF(R(1) .NE.ZERO)  THEN 

AK(l,l)=PROP(5)*(SK(l)-SK(2) ) *SHA1 
AK(l,2)=PROP(5) *SK(2) *SHA1 
IF(R(1) .LT.ZERO)  THEN 
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AK(1,1)=-AK(1,1) 

AK(1,2)=-AK(1,2) 

ENDIF 

ENDIF 

IF(R(2) .NE.ZERO)  THEN 

AK(2,l)=-PROP(5) *SK(2) *SHA1 
AK(2,2)=PROP(5) *(SK(2)-SK(3) )*SHA1 
AK(2,3)=PROP(5) *SK(3) *SHA1 
IF(R(2) .LT.ZERO)  THEN 
AK(2,1)=-AK(2,1) 

AK(2,2)=-AK(2,2) 

AK(2,3)=-AK(2,3) 

ENDIF 

ENDIF 

IF(R(3) .NE.ZERO)  THEN 

AK(3,2)=-PROP(5) *SK(3) *SHA2 
AK(3,3)=PROP(5) *(SK(3) -SK(4) ) *SHA2 
AK(3,4)=PROP(5) *SK(4) *SHA2 
IF(R(3) .LT.ZERO)  THEN 
AK(3,2)=-AK(3,2) 

AK(3,3)=-AK(3,3) 

AK{3,4)=-AK(3,4) 

ENDIF 

ENDIF 

IF(R(4) .NE.ZERO)  THEN 

AK(4,3)=-PROP(5) *SK(4) *SHA2 
AK(4,4)=PROP(5) *(SK(4) -SK(5) ) *SHA2 
IF(R(4) .LT.ZERO)  THEN 
AK(4,3)=-AK(4,3) 

AK(4,4)=-AK(4,4) 

ENDIF 

ENDIF 

C 

C TOTAL  STIFFNESS 

C 

CALL  TSTIFF(AK,BK,SF) 

C 

C CALCULATE  DISPLACEMENTS  OF  FORCE  DISTRIBUTION 

C 

190  DO  200  1=2,5 

IF(INS(I) .EQ.O)  THEN 
NS=0 
else 
NS=NS+1 

IF(INS(I+1) .EQ.O)  THEN 
C IF(NS.EQ.4)  THEN 

C CALL  CROUT(SF,R,4,1,0) 

C GO  TO  200 

C ENDIF 

CALL  SOL(R,SF,NS, I,SS,RR) 

ENDIF 

ENDIF 

200  CONTINUE 
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C 

C CALCULATE  THE  TENDON  FORCES 
C 

SEL(1)=SEL(1)+SK(1) *R(1) 

DO  300  1=2,4 
11=1-1 

300  SEL(I)=SEL(I)+SK(I) *(R(I)-R(I1) ) 

SEL(5)=SEL(5)-SK(5) *R(4) 

C 

C CHECK  IF  THE  RECHECKING  IS  NECESSARY 

C 

DO  305  1=2,5 

IF(INS(I) .NE.IN(I) ) GO  TO  310 
305  CONTINUE 
GO  TO  350 
C 

C RECHECK  THE  UNSLIDING  NODES 

C 

310  DO  320  1=2,5 

IF(INS(I) .EQ.O)  THEN 
IT=IT+1 
GO  TO  70 
ENDIF 

320  CONTINUE 
C 

C COMPARE  INS  AND  INSO 

C 

350  KST=0 

DO  600  1=2,5 

IF(INS(I) .NE.INSO(I) ) KST=1 
600  CONTINUE 
C 

C TURN  ON  STIFFNESS 

C 

C IF(IP.EQ.0.AND.NSEG.EQ.2)  KST=1 

IF(IP.EQ.NSEG.AND.NIT.EQ. 1)  KST=1 
RETURN 
END 

TRAN4 

SUBROUTINE  TRAN4 (T, Q, QL, NDF) 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  T ( 3 , 3 ) , Q (NDF) , QL ( 12 ) 

C DO  100  1=1,14 

ClOO  QL(I)=0.0 

DO  200  1=1,3 

200  QL(1)=QL(1)+T(1,I) *Q(I) 

DO  300  1=2,3 
DO  300  K=4,6 
11=1-1 
KK=K-3 

300  QL(I)=QL(I)+T(II,KK)*Q(K) 

DO  400  1=4,5 
DO  400  K=7,9 
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11=1-3 
KK=K-6 

400  QL(I)=QL(I)+T(II,KK) *Q(K) 

DO  500  1=16,18 
11=1-15 

500  QL(12)=QL(12)+T(1,II) *Q(I) 

DO  1500  L=6,9,3 
DO  1500  J=L,L+2 
DO  1500  K=L+4,L+6 
JJ=J-L+1 
KK=K-L-3 

1500  QL(J)=QL(J)+T(JJ,KK) *Q(K) 

RETURN 
END 

C 

SUBROUTINE  SOL(R,SF,NS,N,SS,RR) 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  R(4) ,SF(4,4) ,RR(NS) ,SS(NS,NS) 

DO  100  1=1, NS 
DO  100  J=1,NS 
II=N-NS-1+I 
JJ=N-NS-1+J 
SS(I, J)=SF(II, JJ) 

100  CONTINUE 

DO  200  1=1, NS 
II=N-NS-1+I 
RR(I)=R(II) 

200  CONTINUE 

CALL  CROUT(SS,RR,NS, 1, 0) 

DO  300  1=1, NS 
II=N-NS-1+I 
R(II)=RR(I) 

300  CONTINUE 
RETURN 
END 



SUBROUTINE  CROUT (A, B, N, LD, KEY) 



C GENERAL  EQUATION  SOLVER  - M.  I.  HOIT  1/17/83 

C SYMMETRIC  AND  NON  SYMMETRIC 
C 

C KEY  = 0 TRIANGULARIZE  AND  SOLVE 

C KEY  = 1 TRIANGULARIZE  ONLY 

C KEY  = 2 FOWARD  AND  BACK  SUBSTITUTE  ONLY 



IMPLICIT  REAL* 8 (A-H,0-Z) 

DIMENSION  A(N,N) ,B(N,LD) 

C 

IM  = 1 

IF (KEY. EQ. 2)  GO  TO  200 
IF(N.EQ.l)  GO  TO  200 
DO  140  J=2,N 
JM  = J - 1 


SOL 


CROUT 
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A(J,1)  = A(J,1)/A(1,1) 

IF(J.EQ.2)  GO  TO  130 
DO  120  1=2, JM 
IM  = I - 1 
DO  100  K=1,IM 

C FORM  L(I,J)  U(I,J)  

A(J,I)  = A(J,I)  - A{J,K)*A(K,I) 

100  A(I,J)  = A(I,J)  - A(I,K) *A(K, J) 

D = A(I,I) 

C CHECK  IF  ZERO  ON  DIAGONAL 

IF(D.NE.O.O)  GO  TO  120 
WRITE (NTM, 1000)  I 
GO  TO  140 

120  A(J,I)  = A(J,I)/D 
130  DO  135  K=1,JM 

135  A(J,J)  = A(J,J)  - A(J,K) *A(K, J) 

140  CONTINUE 

C FOWARD  AND  BACK  SUBSTITUTE 

200  IF(KEY.EQ.l)  GO  TO  900 
DO  400  L=1,LD 

C FORM  Y(I,L) 

IF(N.EQ.l)  GO  TO  300 
DO  210  1=2, N 
IM  = I - 1 
DO  210  K=1,IM 

210  B(I,L)  = B(I,L)  - A(I,K)  *B(K,L) 

C form  X(I,L) 

300  B(N,L)  = B(N,L)/A(N,N) 

IF(N.EQ.l)  GO  TO  400 
IM  = N - 1 
DO  330  1=2, N 
JM  = IM  + 1 
DO  310  J=JM,N 

310  B(IM,L)  = B(IM,L)  - A(IM, J) *B(J,L) 

320  B(IM,L)  = B(IM,L)/A(IM,IM) 

330  IM  = IM  - 1 

C 

400  CONTINUE 
900  RETURN 

1000  FORMATC  ****  ZERO  ON  THE  DIAGONAL  OF  EQUATION' 
* ' NUMBER  ',15,'  ****•/) 

END 


^ SUBROUTINE  RINTll  (NDF, Q,VEL, FE, FD, TIME, FACAL, FACIL, ALFA) 

C SUBROUTINE  FOR  INTERNAL  FORCE  CALCULATIONS, 

C ELEMENT  TYPE  = 11  THREE  DIMENSIONAL  TENDON  ELEMENTS. 

C ********************************************************* 

C 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

DIMENSION  Q(NDF) ,VEL(NDF) ,FE(NDF) ,FD(NDF)  ,FT(12) , FACAL (1)  . 
* FACIL(l) ,TTI(5) 
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COMMON  /INFGR  / NCR, NELS , MFST, IGRHED ( 10 ) , NINFC, LSTAT, 

* LSTF,LSTC,NDOF,DKO,DKT,EPROP(5,200) 
COMMON  /INFEL  / IMEM, KST, LM(18) ,NODE ( 6) , INS (6) , INSO (6)  , 

* XX(6) ,YY(6) ,ZZ(6) , PROP (5) ,AREA,XL(5) ,B(7) ,SEL(5)  , 

* STOT(5) ,TT(3,3) ,T,TI(5) ,QL(12) ,SHA1,SHA2, 

* BK(4,4) ,SK(5) ,IP,IDUM,ANSL,SKP(324) 

COMMON  /STLDPT/  NPNF,NPFF,NPTP,NPP4 ,NPDP,NETP,NESW,NITSR, 

* NSPATT,NSEG 
COMMON  /PRECIS/  ZERO, ONE 

C 

C CALCULATE  STRESS 
C 

if (nitsr.eq. 1)  then 

c compute  factor  of  tendon  force  to  include  

npt=npnf+npff+nptp+npp4+npdp+netp+nesw+nitsr 

nptl=npt-npdp 

fac=facal (nptl) +alfa*facil (npt) 
else 

fac=ZERO 

endif 

do  100  i=l,5 
100  tti(i)=fac*ti(i) 

350  DO  400  1=1,5 
400  STOT(I)=SEL(I)/AREA 
if (nseg.ge. ip)  then 
CALL  FORVER(SEL,B,FT,XX,YY,ZZ) 
else 

call  forver (tti,b, ft,xx,yy, zz) 
endif 
C 

C TRANSFORMATION  FROM  LOCAL  TO  GLOBAL 
C 

CALL  TRAN3 (TT,FT,FE, 18) 

C 

C CALCULATE  THE  DAMPING  FORCE 
C 

DO  500  I=l,NDOF 
500  FD(I)=ZERO 
RETURN 
END 

C****************** ****** *************************** 

C SUBROUTINE  TRAN3  * 

C*************************************************^^ 

SUBROUTINE  TRAN3 (T, FL, F, NDF) 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  T(3,3) ,F(NDF) ,FL(12) 

COMMON/ PRECIS/ ZERO , ONE 
DO  100  1=1, NDF 
100  F(I)=ZERO 

DO  200  1=1,3 
F(I)=T(1,I) *FL(1) 

DO  300  1=16,18 
11=1-15 


200 


300  F(I)=T(1,II) *FL(12) 

DO  400  1=4,6 
11=1-3 

DO  400  K=2,3 
KK=K-1 

400  F(I)=F(I)+T(KK,II) *FL(K) 
DO  500  1=7,9 
11=1-6 

DO  500  K=4,5 
KK=K-3 

500  F(I)=F(I)+T(KK,II) *FL(K) 
DO  600  L=6,9,3 
DO  600  J=L+4,L+6 
DO  600  K=L,L+2 
JJ=J-L-3 
KK=K-L+1 

600  F(J)=F(J)+T(KK,JJ)*FL(K) 
DO  700  1=16,18 
11=1-15 

700  F(1)=T(1,11) *FL(12) 

RETURN 
END 


C SUBROUTINE  FORVER  


SUBROUTINE  FORVER (SEL,B,FT,XX,YY,ZZ) 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  SEL(5) , B (7) , FT (12 ) , XX (6) , YY ( 6) , ZZ (6) 
FT(1)=-SEL(1) 

FT(2)=SEL(1) -SEL(2) *B(3) 

FT(3)=SEL(2)  *B(1) 

FT(4)=SEL(2) *B(3)-SEL(3) 

FT(5)=-FT(3) 

FT(6)=SEL(3) -SEL(4) *B(4) 

FT(7)=-SEL(4) *B(2) 

FT(8)=-SEL(4) *B(5) 

FT(9)=SEL(4) *B(4) -SEL(5) 

FT(10)=-FT(7) 

FT(11)=-FT(8) 

FT(12)=SEL(5) 

IF  (XX(1) .EQ.XX(2) ) THEN 
FT(1)=0.0 

FT(2)=-SEL(2) *B(3) 

ENDIF 


IF  (XX(6) .EQ.XX(5) ) THEN 
FT(9)=SEL(4) *B(4) 
FT(12)=0.0 
ENDIF 
RETURN 
END 


C 
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C 


SUBROUTINE  INITI 1 ( NJT , NDF , RF , FACIL) 
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C 

C 

C 

C 

C 


C 

C— 

C 

c — 

c — 
c — 

C 

C 

C 

600 

C 

C 

C 


C 

C 

C 

C 

c 

c 

c 


***************************************************** 

SUBROUTINE  TO  COMPUTE  THE  INITIAL  FORCE 

ELEMENT  TYPE  = 11  THREE  DIMENSIONAL  TENDON  ELEMENTS. 

***************************************iiit********** 

IMPLICIT  DOUBLE  PRECISION (A-H,0-Z) 

COMMON  /INFGR  / NCR, NELS ,MFST, IGRHED (10) , NINFC, LSTAT, 

* LSTF,LSTC,NDOF,DKO,DKT,EPROP(5,200) 
COMMON  /INFEL  / IMEM, KST, LM( 18) , NODE (6) , INS (6) , INSO (6) , 

* XX(6) ,YY(6) ,ZZ(6) ,PROP(5) , AREA, XL(5) , B (7) ,SEL(5)  , 

* STOT(5) ,TT(3,3) ,T,TI(5) ,QL(12) ,SHA1,SHA2, 

* BK(4,4) ,SK(5) ,IP,IDUM,ANSL,SKP(324) 

COMMON  /STLDPT/  NPNF,NPFF,NPTP,NPP4 ,NPDP,NETP,NESW,NITSR, 

* NSPATT,NSEG 
COMMON  /PRECIS/  ZERO, ONE 

DIMENSION  RF(NDF) ,FI(12) ,P(5) ,FACIL(1) 

-CHECK  IF  INITIAL  FORCE  IS  TO  BE  INCLUDED  

IF(NITSR.EQ.l)  THEN 

-compute  factor  of  tendon  force  to  include  

npt=npnf+npff+nptp+npp4+npdp+netp+nesw+nitsr 
-for  facal,  nptl  = npt  - npdp 
-fac  = facal (nptl)  + alfa*facil (npt) 
fac=facil (npt) 

CALCULATE  THE  EQUVALENT  INITIAL  FORCE 

DO  600  1=1,5 
P(I)=FAC*TI(I) 

CALL  FORVER(P,B,FI,XX,YY,ZZ) 

TRANSFORMATION  FROM  LOCAL  TO  GLOBAL 

CALL  TRAN3(TT,FI,RF,18) 

ENDIF 

RETURN 

END 


SUBROUTINE  OUTSll  (KPR,TIME) 


**************************************************** 
SUBROUTINE  TO  PRINT  TIME  HISTORY  OF  CURRENT  STATE, 
ELEMENT  TYPE  = 11  THREE  DIMENSIONAL  TENDON  ELEMENTS. 
**************************************************^^ 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

COMMON  /TAPES  / NIU, NOU, NTl , NT2 , NT3 , NT4 , NT5 , NT6 , NT7 , NT8 , 

* NT9,NT10,NT11,NT12,NT13,NT14,NT15,NT16 
COMMON  /INFGR  / NGR, NELS , MFST, IGRHED ( 10) , NINFC, LSTAT, 

* LSTF , LSTC , NDOF , DKO , DKT , EPROP (5,200) 
COMMON  /INFEL  / IMEM, KST, LM ( 18 ) , NODE ( 6) , INS ( 6 ) , INSO ( 6) , 

* XX(6) ,YY(6) ,ZZ(6) , PROP (5) ,AREA,XL(5) ,B(7) ,SEL(5) , 
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* ST0T(5) ,TT(3,3) ,T,TI(5) ,QL(12) ,SHA1,SHA2, 

* BK(4,4) ,SK(5) ,IP,IDUM,ANSL,SKP(324) 

C 

IF  (IMEM.EQ.MFST)  KHED  = 0 
IF  (KPR.EQ.O)  RETURN 
C 

IF  (KHED.NE.O)  GO  TO  100 
C 

KHED  = 1 

KKPR  = lABS(KPR) 

WRITE (NOU, 2000)  KKPR,TIME, IGRHED 
100  WRITE (NOU, 2100)  IMEM,NODE ( 1) , NODE (6) 

DO  200  1=1,5 

200  WRITE (NOU, 2010)  NODE (I) ,NODE (I+l) , SEL(I) , STOT (I) 

WRITE (NOU, 2200)  INS 

2200  FORMAT (5X, 13HSLIDING  INDEX, 5X, 612 ) 

C 

2000  FORMAT  (///18H  RESULTS  FOR  GROUP, 13, 

* 44H  (THREE  DIMENSIONAL  TENDON  ELEMENTS),  TIME  = Ell. 4// 

* 5X,10A4//) 

2100  FORMAT  (5X, 8HELEM  NO. , 15, 5X, 5HNODE1, 15, 5X, 5HNODE6, 15// 

* 5X,4HN0DE,1X,3H , IX, 4HNODE , 5X, lOHTEN.  FORCE, 

* 5X, lOHTEN. STRESS/) 

2010  FORMAT  (2I9,2E15.4) 

RETURN 

END 
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